Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	598,486	T→C	7.5%	intergenic (+35/+247)	tufA → / ← yhaA	elongation factor EF1A/M20D subfamily peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	598,486	0	T	C	7.5%	190.6 / 6.4	67	intergenic (+35/+247)	tufA/yhaA	elongation factor EF1A/M20D subfamily peptidase
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

GTAGAATTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGTTGAGTCTTGGA  >  CP000730/598347‑598624                                                                                                                                            |                                                                                                                                            gTAGAATTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATc                                                                                                                                            >  1:260005/1‑141 (MQ=255)     aaTTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTc                                                                                                                                        <  1:512507/141‑1 (MQ=255)          aTCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTcgattg                                                                                                                                   >  1:181954/1‑136 (MQ=255)           tCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGa                                                                                                                                  >  2:580973/1‑141 (MQ=255)             gCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgg                                                                                                                                <  2:571955/141‑1 (MQ=255)              ccccAATAGCGATTGAAGACGGGACACGTTGATCAATACGGGGAGGGCAACGAACTGTAGAATCAGGGGATGTATCTGAAATCATTAAATAATTTCTAAGATATTAGATTTTATATAAAAAGACAATCGCTCAATCGAggg                                                                                                                               >  2:512179/3‑141 (MQ=255)               gattggagcgATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                                                              >  2:85550/8‑141 (MQ=255)                  aaTCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCt                                                                                                                           <  2:403015/141‑1 (MQ=255)                   aTCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtt                                                                                                                          >  2:172191/1‑141 (MQ=255)                   aTCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCtt                                                                                                                          <  1:237313/141‑1 (MQ=255)                     cgcgATTGAAGAAGGTACTCGTTTATCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttct                                                                                                                        <  2:137122/141‑3 (MQ=255)                       cgATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtt‑tttt                                                                                                                     >  1:86984/1‑141 (MQ=255)                           tGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAAt                                                                                                                 <  2:514785/141‑1 (MQ=255)                             aaGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAAtgt                                                                                                               <  2:506244/141‑1 (MQ=255)                             aaGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatat                                                                                                                <  2:181954/141‑4 (MQ=37)                                    tACTCGTTTCTCAATCCGTGAAGGTTGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAAtt                                                                                                        <  2:335850/141‑1 (MQ=255)                                    tACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAAtt                                                                                                        <  1:152796/141‑1 (MQ=255)                                    tACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAAtt                                                                                                        <  1:69948/141‑1 (MQ=255)                                    tACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAAtt                                                                                                        <  1:152813/141‑1 (MQ=255)                                          tgtctcAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTGATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATTTGTAAATTTTGTaa                                                                                                  <  2:244452/139‑1 (MQ=255)                                           ttCTCAATCAGAGAACGAGGACTATCTTTATCATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGTGTCTT‑TTTTTAATGTGTAAATTTTGTAAt                                                                                                 >  2:174313/1‑141 (MQ=255)                                              tcAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGc                                                                                              >  1:240347/1‑141 (MQ=255)                                                aaTCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTa                                                                                            <  1:688093/141‑1 (MQ=255)                                                 aTCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTAt                                                                                           <  2:131527/141‑1 (MQ=255)                                                  tCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTAtt                                                                                          <  2:317161/141‑1 (MQ=255)                                                  tCCGTGAAGGTGGACGAACTGTAGGATCAGGCGTTGTTACTGAAAACATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCATCAATCAAGTGGTCTT‑TTTTTAATGCGTAAATTTTGTAATGGCTAtt                                                                                          >  2:210829/1‑141 (MQ=255)                                                   ccGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTc                                                                                         >  1:642165/1‑141 (MQ=255)                                                    cgcgAAGGTGGACGTACTGTAGGAACAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCg                                                                                        <  1:348068/138‑1 (MQ=255)                                                       gAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGAtt                                                                                     >  1:516572/1‑141 (MQ=255)                                                       gAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGAtt                                                                                     <  2:204008/141‑1 (MQ=255)                                                         aGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAAtgt                                                                                                               >  1:249359/1‑113 (MQ=255)                                                         aGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAAtgt                                                                                                               <  2:249359/113‑1 (MQ=255)                                                         aGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTa                                                                                   <  2:136431/141‑1 (MQ=255)                                                          ggTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGTTTTTATATAAAAAGAAGATTCATCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTTTAATGGCTAttaattttat                                                                                  >  2:209826/1‑133 (MQ=255)                                                                     tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                                                       >  1:320947/1‑141 (MQ=255)                                                                      gTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataat                                                                      >  1:164270/1‑141 (MQ=255)                                                                         ggATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtga                                                                   >  1:204506/1‑141 (MQ=255)                                                                            tCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                                                                >  1:448499/1‑141 (MQ=255)                                                                            tCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                                                                <  2:478309/141‑1 (MQ=255)                                                                              aGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaaa                                                              <  1:680191/141‑1 (MQ=255)                                                                               ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTNTTCGATTTAGAAGAAc                                                                            <  1:461534/126‑1 (MQ=255)                                                                               ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAAc                                                                            >  2:461534/1‑126 (MQ=255)                                                                                    tgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga                                                        <  2:9321/141‑1 (MQ=255)                                                                                     gttACTGAAATCCTTATATATTTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac                                                       <  2:571310/141‑1 (MQ=255)                                                                                     gttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac                                                       <  1:69616/141‑1 (MQ=255)                                                                                      ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTc                                                                                                                            >  1:502564/1‑72 (MQ=255)                                                                                      ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTc                                                                                                                            <  2:502564/72‑1 (MQ=255)                                                                                       tACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                                                     <  1:114590/141‑1 (MQ=255)                                                                                       tACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                                                     <  1:343932/141‑1 (MQ=255)                                                                                        acttaaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa                                                    <  1:237553/137‑1 (MQ=255)                                                                                          agAAATTATTTAATGATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGAtt                                                                                     <  2:269222/105‑1 (MQ=255)                                                                                          agAAATTATTTAATGATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGAtt                                                                                     >  1:269222/2‑106 (MQ=255)                                                                                                  ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAACCTtat                                          <  2:320947/141‑1 (MQ=255)                                                                                                    aaaTAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTAGTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATaa                                        <  1:544438/141‑1 (MQ=255)                                                                                                           ttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGataa                                 >  2:572274/1‑141 (MQ=255)                                                                                                                   tctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGttt                         >  1:486635/1‑141 (MQ=255)                                                                                                                    ctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTta                        >  2:149198/1‑141 (MQ=255)                                                                                                                    ctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTta                        >  2:361589/1‑141 (MQ=255)                                                                                                                     ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTtaa                       <  2:240347/141‑1 (MQ=255)                                                                                                                        gATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTtaaata                    >  2:575241/1‑141 (MQ=255)                                                                                                                         aTTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTtaaataa                   >  2:652957/1‑141 (MQ=255)                                                                                                                               tataAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAAttgt             >  2:252224/1‑141 (MQ=255)                                                                                                                               cattaaaaaaAGACCCCTCGATTGAGGGATCTTCTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAAttg              >  1:583667/5‑141 (MQ=37)                                                                                                                                 taAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAAttgttg           <  2:86984/141‑1 (MQ=255)                                                                                                                                    aaaGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGTTGAGt        >  2:220800/1‑141 (MQ=255)                                                                                                                                     aagaagATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGTTGAGTc       >  1:183534/1‑141 (MQ=255)                                                                                                                                          gaTCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGTTGAGTCTTGGa  >  1:514341/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                            GTAGAATTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGTTGAGTCTTGGA  >  CP000730/598347‑598624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑