Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,153,598	A→G	5.7%	Y216H (TAT→CAT)	scrR ←	sucrose operon repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,153,598	0	A	G	5.7%	212.3 / 3.1	70	Y216H (TAT→CAT)	scrR	sucrose operon repressor
Reads supporting (aligned to +/- strand):  ref base A (25/41);  new base G (2/2);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGATGTGGTTTCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACTTGTTGATATCCCTGTTGA  >  CP000730/2153459‑2153737                                                                                                                                            |                                                                                                                                            tGATGTGGTTTCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATaa                                                                                                                                            <  1:605817/141‑1 (MQ=255)     gtgGTTTCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTaa                                                                                                                                        <  1:369033/141‑1 (MQ=255)       ggTTTCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTaaaa                                                                                                                                      <  1:171278/141‑1 (MQ=255)       ggTTTCATAACATCGTTTTTATCAGAATAATATTTATAGGCAGCTAATGCACTCAAATCGTTTGCTGCGACAACCGCATCTACTGGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTaaaa                                                                                                                                      <  1:493546/141‑1 (MQ=255)        gTTTCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAAt                                                                                                                                     >  1:643645/1‑141 (MQ=255)         tttCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAAtt                                                                                                                                    <  1:161865/141‑1 (MQ=255)         tttCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAAtt                                                                                                                                    <  1:394138/141‑1 (MQ=255)            cATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGt                                                                                                                                 <  2:689726/141‑1 (MQ=255)                       ttttATCAGAATAATATTTATCGGCAGCTAATGCGATCGTATCAGTTGCTCCAACAACCGCATATACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAGTTTGTTTCATGTatat                                                                                                                      >  2:496819/1‑141 (MQ=255)                           aTCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAgg                                                                                                                  <  2:603509/141‑1 (MQ=255)                             cAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGtt                                                                                                                <  2:663510/141‑1 (MQ=255)                              aGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATGAGTAAAATTTGTTTCATGTATATTAGGttt                                                                                                               >  1:341089/1‑141 (MQ=255)                              aGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATGAGTAAAATTTGTTTCATGTATATTAGGttt                                                                                                               >  1:557749/1‑141 (MQ=255)                               gAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTACTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTa                                                                                                              >  2:454740/1‑141 (MQ=255)                                    atatTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTg                                                                                                         <  2:194997/141‑1 (MQ=255)                                         tataGGCAGCTAATGCAATCGTTTCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTAtt                                                                                                    <  1:590526/141‑1 (MQ=255)                                             ggCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGc                                                                                                <  1:546120/141‑1 (MQ=255)                                               cAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTa                                                                                              <  2:45786/141‑1 (MQ=255)                                                aGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTaa                                                                                             >  1:47400/1‑141 (MQ=255)                                                aGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTaa                                                                                             >  1:361864/1‑141 (MQ=255)                                                aGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTaa                                                                                             <  2:453749/141‑1 (MQ=255)                                                    aaTGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTgg                                                                                         <  1:683478/141‑1 (MQ=255)                                                       gCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCa                                                                                      >  1:202697/1‑141 (MQ=255)                                                       gCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCa                                                                                      >  2:364945/1‑141 (MQ=255)                                                            ggTATCAGTTGCTCCAACAACCGCATCTGCTTGCTCCACAGTTTCTATATCATTTGCAACATCGTTTTTTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTGGCTAACTGGTCAAGTaa                                                                                 <  2:167440/140‑1 (MQ=255)                                                            cGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTaa                                                                                 <  2:705162/141‑1 (MQ=255)                                                             gTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTaaa                                                                                <  1:271122/141‑1 (MQ=255)                                                               aTCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAAcc                                                                              <  2:647971/141‑1 (MQ=255)                                                                tCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAGGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCa                                                                             <  2:60166/141‑1 (MQ=255)                                                                  aGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACg                                                                           >  1:702355/1‑141 (MQ=255)                                                                    ttGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGtt                                                                         <  1:554467/141‑1 (MQ=255)                                                                     tGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGttt                                                                        <  2:361864/141‑1 (MQ=255)                                                                       cTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTc                                                                      >  1:462433/1‑141 (MQ=255)                                                                            acaacCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGt                                                                 >  1:305783/1‑141 (MQ=255)                                                                              aacCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGtat                                                               >  1:270901/1‑141 (MQ=255)                                                                                 cGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATAcc                                                            <  1:382172/141‑1 (MQ=255)                                                                                 cGCATCTACTTGCTCCACATTTGCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTCGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATAcc                                                            <  1:588892/141‑1 (MQ=255)                                                                                   cATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCaa                                                          >  2:390352/1‑141 (MQ=255)                                                                                     tCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTAt                                                                                                     >  1:593976/1‑96 (MQ=255)                                                                                     tCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTAt                                                                                                     <  2:593976/96‑1 (MQ=255)                                                                                       tACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGc                                                      <  1:74586/141‑1 (MQ=255)                                                                                       tACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGc                                                      >  1:694561/1‑141 (MQ=255)                                                                                        aCTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCa                                                     <  2:220283/141‑1 (MQ=255)                                                                                        aCTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCa                                                     >  1:401692/1‑141 (MQ=255)                                                                                         cTTGCTTCACATCTTCCAAAACATTTGCAACTTCTTTTTGTGCTTCCACATACGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTCTGTATACCAACTGCaa                                                    <  2:405275/141‑1 (MQ=255)                                                                                            gCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAatat                                                 <  1:647304/141‑1 (MQ=255)                                                                                                cacaTTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATCTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCttt                                             <  1:108726/141‑1 (MQ=255)                                                                                                 acaTTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCtttt                                            <  2:202697/141‑1 (MQ=255)                                                                                                     tttCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCAc                                        >  1:99648/1‑141 (MQ=255)                                                                                                      ttCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACt                                       >  1:645162/1‑141 (MQ=255)                                                                                                      ttCCAAAACATTTGCAACATCTTTTTGTACTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACt                                       <  1:454740/141‑1 (MQ=255)                                                                                                         cAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTac                                    <  2:694602/141‑1 (MQ=255)                                                                                                          aaaaCATTGGCAACATCTTTTTGAGCTTGCTGATAAGTAAAATTAGTTAAATTTATAGTAGATTTAATTTGATATTTAGCTAACTAGTCAAGTAAACCACGTTTTCTATGTATACCAAGTGCAATTTCTTTTTCCCTTaca                                   >  2:74954/1‑141 (MQ=255)                                                                                                            aaCATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACt                                 <  2:664287/141‑1 (MQ=255)                                                                                                               aTTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAGCTGCAATATCTTTTTCACTTACACTaaa                              <  1:95272/141‑1 (MQ=255)                                                                                                               aTTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTaaa                              <  2:345499/141‑1 (MQ=255)                                                                                                                  tGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAacac                           >  2:550258/1‑141 (MQ=255)                                                                                                                   gCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACt                          >  1:623642/1‑141 (MQ=255)                                                                                                                   gCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACt                          >  1:460709/1‑141 (MQ=255)                                                                                                                    tgcagtgtctgtatgCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACtt                         <  2:36826/128‑1 (MQ=255)                                                                                                                    cAACATCTTTTTGTGCTTCCACATGAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACtt                         <  2:341089/141‑1 (MQ=255)                                                                                                                    cAACATCTTTTTGTGCTTCCACATGAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACtt                         <  2:557749/141‑1 (MQ=255)                                                                                                                    cAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACtt                         <  1:489555/141‑1 (MQ=255)                                                                                                                    cAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACtt                         <  2:702355/141‑1 (MQ=255)                                                                                                                    cAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATATATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACATTAAAGACtt                         >  2:201457/1‑141 (MQ=255)                                                                                                                             tttGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACttgtt                <  1:507575/141‑1 (MQ=255)                                                                                                                             tttGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACttgtt                >  2:679239/1‑141 (MQ=255)                                                                                                                               tgtgCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTCTGTATACCAACTGCAATATCTTTTTCACTCACACTAAACACTTCAACTTGTTGa              >  1:200871/1‑141 (MQ=255)                                                                                                                                gtgCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACTTGTTGat             <  2:474387/141‑1 (MQ=255)                                                                                                                                     tCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACTTGTTGatatgct        >  2:588763/1‑138 (MQ=255)                                                                                                                                           tAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACTTGTTGATATCCCTGTTGa  <  1:550258/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            TGATGTGGTTTCATAACATCTTTTTTATCAGAATAATATTTATAGGCAGCTAATGCAATCGTATCAGTTGCTCCAACAACCGCATCTACTTGCTCCACATTTTCCAAAACATTTGCAACATCTTTTTGTGCTTCCACATAAGTAAAATTTGTTTCATGTATATTAGGTTTAATTTGGTATTTAGCTAACTGGTCAAGTAAACCACGTTTTCTATGTATACCAACTGCAATATCTTTTTCACTTACACTAAACACTTCAACTTGTTGATATCCCTGTTGA  >  CP000730/2153459‑2153737

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