Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,985	0	T	C	25.5%	201.8 / 68.5	98	intergenic (‑195/+113)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (24/1);  total (61/37)
Fisher's exact test for biased strand distribution p-value = 2.48e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084                                                                                                   |                                                                                                           gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCGCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  1:1420635/1‑100 (MQ=255)      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      <  1:935223/100‑1 (MQ=255)         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  2:265270/100‑1 (MQ=255)         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  1:186953/100‑1 (MQ=255)          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                  <  1:1350381/100‑1 (MQ=255)          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                  <  1:1346543/100‑1 (MQ=255)           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                 <  2:6838/100‑1 (MQ=255)            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                                <  2:1038347/100‑3 (MQ=25)                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                           >  1:47946/1‑100 (MQ=255)                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                          >  2:994385/1‑100 (MQ=255)                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                         <  2:705535/100‑1 (MQ=255)                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                        >  2:715645/1‑100 (MQ=255)                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       <  2:677211/100‑1 (MQ=255)                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                      <  2:1460629/100‑1 (MQ=255)                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                      <  1:984862/100‑1 (MQ=255)                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                      <  1:1199226/100‑1 (MQ=255)                       agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                                     >  1:242053/1‑100 (MQ=255)                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                 <  2:61338/100‑1 (MQ=255)                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCATACAAAATCGTGCAAATTc                                                                                 <  1:120279/100‑1 (MQ=255)                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                                <  1:624266/100‑1 (MQ=255)                             acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                               >  1:1069831/1‑100 (MQ=255)                                atctcACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                            <  1:1918157/96‑1 (MQ=255)                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                          >  1:1246908/1‑100 (MQ=255)                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACagat                                                                        >  2:1918145/1‑97 (MQ=255)                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                        >  1:225974/1‑100 (MQ=255)                                      ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑a                                                                    >  1:1612474/1‑100 (MQ=255)                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                       <  2:598050/98‑1 (MQ=255)                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                  <  1:651302/100‑1 (MQ=255)                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                            >  2:538928/1‑100 (MQ=255)                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                            >  2:544868/1‑100 (MQ=255)                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                            >  2:1920274/1‑100 (MQ=255)                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑tt                                                           >  2:1043804/1‑100 (MQ=255)                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑ttt                                                          <  2:1680332/100‑1 (MQ=255)                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTc                                                         >  2:527322/1‑100 (MQ=255)                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                                      >  2:1877862/1‑100 (MQ=255)                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑Gt                                                     <  2:483381/100‑1 (MQ=255)                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTa                                                    <  2:1069825/100‑1 (MQ=255)                                                  ttccgatctaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAg                                                   <  2:981688/91‑1 (MQ=255)                                                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGc                                                 >  1:1866184/1‑100 (MQ=255)                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGcc                                                <  1:544867/100‑1 (MQ=255)                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCt                                               <  2:1572388/100‑1 (MQ=255)                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTg                                              >  2:1460997/1‑100 (MQ=255)                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                             >  2:1624261/1‑100 (MQ=255)                                                        gACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                             >  1:1325610/1‑100 (MQ=255)                                                        gACAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                             >  1:1325647/1‑100 (MQ=255)                                                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAgagatcggaa                                          >  1:981690/1‑91 (MQ=255)                                                            aaGTGCGCTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAg                                         >  1:1603951/1‑100 (MQ=255)                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                        <  1:538927/100‑1 (MQ=255)                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                    <  1:1814239/100‑1 (MQ=255)                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  2:1478234/1‑100 (MQ=255)                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  1:1944092/1‑100 (MQ=255)                                                                   ccgatcttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:2018964/94‑1 (MQ=255)                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  2:376388/1‑100 (MQ=255)                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  1:343608/1‑100 (MQ=255)                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  2:1671658/1‑100 (MQ=255)                                                                    attttgtaggCCGGATAAGGCGTTAACGCCGCATCCTGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                                 >  2:471679/10‑100 (MQ=17)                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  1:1040961/7‑100 (MQ=17)                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            <  1:1252198/100‑1 (MQ=255)                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGC‑‑A‑GA‑GA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:560033/5‑100 (MQ=17)                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgagatcgg                           >  2:2018953/1‑93 (MQ=255)                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                           >  1:236987/1‑100 (MQ=255)                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:730057/3‑100 (MQ=18)                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          <  1:1624273/100‑1 (MQ=255)                                                                            ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:73401/11‑100 (MQ=21)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                        >  1:563293/2‑98 (MQ=21)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1132936/2‑100 (MQ=21)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:217134/2‑100 (MQ=17)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  1:1329195/2‑100 (MQ=21)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  1:643340/2‑100 (MQ=21)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  1:572536/2‑98 (MQ=21)                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:1906487/1‑100 (MQ=255)                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:1165565/1‑100 (MQ=255)                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  1:335593/1‑98 (MQ=18)                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                       >  1:1336869/1‑97 (MQ=21)                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAACGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  2:1077409/9‑98 (MQ=17)                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:1714715/9‑98 (MQ=21)                                                                               cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  2:618536/8‑97 (MQ=18)                                                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:1485094/8‑97 (MQ=18)                                                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:192989/1‑98 (MQ=255)                                                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  1:1143073/7‑98 (MQ=21)                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaggtt                     >  2:1204936/7‑96 (MQ=17)                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:1723011/7‑100 (MQ=17)                                                                                 gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                    <  1:1671670/100‑1 (MQ=255)                                                                                   taaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctcgc                  >  2:1025950/4‑93 (MQ=17)                                                                                   taaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctcgc                  >  2:1088188/4‑93 (MQ=18)                                                                                    aaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCATTtatgcgt                >  2:279898/3‑98 (MQ=25)                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                <  1:527321/100‑1 (MQ=255)                                                                                       gcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgttt              >  2:1336865/1‑100 (MQ=18)                                                                                       gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              >  2:940242/1‑100 (MQ=255)                                                                                       gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              >  2:940427/1‑100 (MQ=255)                                                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             >  1:919810/1‑100 (MQ=255)                                                                                         gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            >  2:1272430/1‑100 (MQ=255)                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:1944080/100‑1 (MQ=255)                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:919808/100‑1 (MQ=255)                                                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:1043809/100‑1 (MQ=255)                                                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  2:751675/100‑1 (MQ=255)                                                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  1:783479/100‑1 (MQ=255)                                                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  1:1732902/100‑1 (MQ=255)                                                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:255911/100‑1 (MQ=255)                                                                                                   tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  <  2:1603939/100‑1 (MQ=255)                                                                                                   |                                                                                                           GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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