Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,162,120	0	T	G	30.2%	92.2 / 10.7	54	G219G (GGT→GGG)	mdtD	putative multidrug efflux pump MdtD
Reads supporting (aligned to +/- strand):  ref base T (10/27);  new base G (16/0);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 3.58e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTCATGCCGGACTGATG  >  NC_000913/2161884‑2162369                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            ttACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCCCTGG‑GCTGGGCGGCCGTAAAGGGACAGGttttttgcc                                                                                                                                                                                                                                               >  1:142348/1‑245 (MQ=255)   aCCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATcgccg                                                                                                                                                                                                                                              <  1:286351/250‑1 (MQ=255)     cctgTCAGGTCGGTCCGCTGCCCGGTCCGGCGCACGGCGGTCTGCTGGTGGAGTACTCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACAGTGCTGTTAATTCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTTGCATGGCGGTATTAC‑CCC‑CGGCGCTGGACGGCAGTTAAGTTACAGGTTTATCGCCGCTg                                                                                                                                                                                                                                           <  2:260047/248‑1 (MQ=255)      cGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCCTGGCGGTATTAA‑ACC‑TGGCGCTGGACGGCAGGAAAGGGACAGGGTTTTCGCCGCTgg                                                                                                                                                                                                                                          >  1:79733/1‑250 (MQ=255)               cggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCCTGGCGGGATTTA‑CCC‑CGGCGCCGGGCGGCGGTAAAAGTTAAGGTTTATTGCCGCttgcggttgccg                                                                                                                                                                                                                                 >  2:45337/1‑240 (MQ=255)                ggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTA‑CCC‑TGGGGCTGGACGGCAGTAAAAGGACAGGGTTTTCCGCGCTGGAGATTGCAgg                                                                                                                                                                                                                                >  1:231328/1‑251 (MQ=255)                  tcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTA‑CCC‑TGGCGCTGGGCGGCAGTAAAGGTACAGGGTTATCGCCCCTGACGATTGCCGGcc                                                                                                                                                                                                                              >  1:346106/1‑251 (MQ=255)                            ggTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGt                                                                                                                                                                                                                    <  1:247222/251‑1 (MQ=255)                                       cggcggTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtggt                                                                                                                                                                                                           <  1:174384/249‑1 (MQ=255)                                       cggcggTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtggt                                                                                                                                                                                                           <  2:145274/249‑1 (MQ=255)                                                   ggtggAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑ACC‑TGGGGCTGGACGGCAGTAAAGGGACAGGGTTTTCGCCGCTGACGATTTCAAGGCTGGGCGCCGTTGGGGGGGTGGGACTGGTGCttt                                                                                                                                                                                             >  1:451896/1‑251 (MQ=255)                                                   ggtggAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATCAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGGACAGGGTTTTCGCCCCCGGCGCTTGCAGGCCTGGTCGCCGGTGGCGTGGTGGcgctcgtggttt                                                                                                                                                                                             >  2:113138/1‑240 (MQ=255)                                                      ggAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGTAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATc                                                                                                                                                                                          <  2:228057/251‑1 (MQ=255)                                                        aGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATctg                                                                                                                                                                                        <  2:79733/251‑1 (MQ=255)                                                          tACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCCGTAAAGGGACAGGGTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGGGGCACTGGGGCTTTTTctgct                                                                                                                                                                                      >  1:53197/1‑251 (MQ=255)                                                                  gTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCa                                                                                                                                                                              <  1:103532/251‑1 (MQ=255)                                                                              cTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCCCATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGTTTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGGACAGGTTTTTACCCGCTGACGATTGCAGGCCTGGGCTCAGTTGGGGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCCGAAAAAACAAcc                                                                                                                                                                  >  1:397702/1‑251 (MQ=255)                                                                                       cAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTTGCATGGCGGGATTAA‑ACC‑CGGCGCCGGGCGGCAGTAAAGGGAAAGGTTTATCGCCGCTGACGGTTGCAGGGCTGGTCGCAGTTGGGGGGGTGGGACTGGTGCTTTATCTGCTGTACGGCCGAAAAAACAAACGTGGCCTGt                                                                                                                                                         >  1:254914/1‑251 (MQ=255)                                                                                              ccGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCt                                                                                                                                                  <  2:337475/251‑1 (MQ=255)                                                                                               cGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTg                                                                                                                                                 <  1:362991/251‑1 (MQ=255)                                                                                                   ggggATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑CGGGGGTGGCCGGCAGAAAAGGAACCGGTTTACCGCCGCTGAACATTGGAGGCCTGGGTGCCGTTGGTGTGGGGGCCCCGGTGCTTTTTCTGCTGCACGCCAGAAATAAAAACCGGGCCCTGTTCAGGCTGaaa                                                                                                                                              >  1:451297/1‑250 (MQ=255)                                                                                                    gggATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACt                                                                                                                                            <  1:349165/251‑1 (MQ=255)                                                                                                            cGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGt                                                                                                                                     <  1:326781/250‑1 (MQ=255)                                                                                                             ggTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAAACCC‑GGGCG‑TGGACGGCAGGAAAGGGACAGGGTTATTGCCCCCGACGCTTGCAGGCCCGGGCGCAGTTGGGGGGGTGGCACTGGTGCGTTTTTTTCTTCCCGCCAGAAATAACAAACCTGCCCCGTTTACTCTGAAACTGTTCCGtcc                                                                                                                                   >  2:293325/1‑249 (MQ=255)                                                                                                                gCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTc                                                                                                                                 <  1:148044/250‑1 (MQ=255)                                                                                                                     cGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTGTTAGTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑GGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGTGCTGGTGCGTTATCTGGTGCGCGCCAGAAATAGCAACCGGGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTAcc                                                                                                                            <  1:45826/250‑1 (MQ=255)                                                                                                                             tGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAA‑ACC‑TGGGGCTGGACGGCAGTAAAAGGACAGGGTTTTCGCCGCTGACGGTTGCAGGCCTGGGCGCCGGTTGCGGGGGGGGACTGGGGCTTTTTTTGCTGCACGCCAGAAATAAAAACCGGGCCCTGGTCAGTCTTAAAATGGTCCGTAACCGTACCTTTTCGCTg                                                                                                                   >  1:259715/1‑251 (MQ=255)                                                                                                                                      tGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAA‑CCC‑TGGCGCTGGACGGCAGGAAAGGGACAGGGTTATTGCCGCTGACGATTGCAGGCCTGGGCGCAGTTGGCGTGGTGGCACTGGGGCTTTATCTGCTGCCCGCCAGAAATAAAAACCGGGCCCTGTTCAGTCTGAAACTGTTTCGTACTCGTACCTTTTCGCTGGGCCTGGCg                                                                                                          >  1:393724/1‑251 (MQ=255)                                                                                                                                        ccGAACTACACCCTGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑CGGCGCTGGACGGCAGTAAAGGGACAGGGTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGGGTGGGGGCACTGGTGCTTTATCTTCTGCACGCCAGAAATAAAAACCGTGCCCTGTTCAGTCTGAAAATGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCggg                                                                                                        >  1:97804/1‑251 (MQ=255)                                                                                                                                                       cAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACg                                                                                         <  2:62158/251‑1 (MQ=255)                                                                                                                                                       cAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGGACAGGGTTATCGCCGCCGACGATTGCAGGCCTGGGCGCAGGTGGCGGGGTGGCACTGGGGCTTTATCTGCTGCACGCCAGAAATAACAAACCTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCgggggggttttccgggcg                                                                                         >  1:119505/1‑237 (MQ=255)                                                                                                                                                            gcggcgCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGGGGTATTAA‑CCC‑TGGGGCTGGACGGCAGTAAAGGGACAGGGTTATCGCCGCTGGCGAATGGAGGGCTGGTCGCAGGTGGCGTGGGGGCCCTGGTGGTTTATTTGGTGGACGCCCGAAATAAAAACCGTGCCCTGTTTAGTCTGAAACTGTTCCGTAATCGGACCTTTTCGCTGGGCCCGgggggggggcttggcggggggattt                                                                                    >  1:59045/1‑232 (MQ=255)                                                                                                                                                            gcggcgCTTTGATCTCTCCGGATGTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCGGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTg                                                                                    <  1:70411/251‑1 (MQ=255)                                                                                                                                                                  cTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGGACAGGGTTATCGCCGCTGACGATTGGCGGCCTGGTCGCAGTTGGCGTGGGGGCACTGGGGCTTTTTCTGCTGCACGCCAGGAATAACAACCGGGCCCTGTTCAGTCTGGAACTGTTCCGTAATCGTACCCTTTCGCTGGGCCTGGCGGGGGGCTTTGGCgggaggattgggcggg                                                                              >  1:138325/1‑237 (MQ=255)                                                                                                                                                                     tGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGGACAGGTTTATTGCCGCTGACGATTGCAGGCCTGGTCGCAATTGGCGCGGTGGCACTGGTGCTTTATCTGCTGCACCCCAGAAATAACCACCGTGCCCCGTTCAGTCTGAAAATGTTCCGCACTCGTACCTTTTCGCTTGGCCTTGAGGGGGGCTTTGCCCGGCGTTTTGGccggggcc                                                                           >  1:42116/1‑244 (MQ=255)                                                                                                                                                                               ggATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCtt                                                                <  2:359806/251‑1 (MQ=255)                                                                                                                                                                                gATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCttt                                                               <  1:347192/251‑1 (MQ=255)                                                                                                                                                                                gATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGAGTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑GTGCCCttt                                                               <  1:312886/251‑1 (MQ=255)                                                                                                                                                                                     ttATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGa                                                           <  1:426090/250‑1 (MQ=255)                                                                                                                                                                                         tGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACAcc                                                       <  1:404204/250‑1 (MQ=255)                                                                                                                                                                                              gCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCTTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGtttt                                                 <  1:465729/251‑1 (MQ=255)                                                                                                                                                                                               cGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTc                                                <  1:48433/251‑1 (MQ=255)                                                                                                                                                                                                ggTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGCCGGCAGTAAAGGTAAAGGTTTATCGCCGCTGACCATTGCAGGCCTGGTCGCCGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCCGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCCGGCGGGGGGGTTTGCCGGACGGATTGGCCGTGGGATGGTGGCCCTTTTGGACCCCGGtttac                                                >  2:439083/1‑249 (MQ=255)                                                                                                                                                                                                               aTTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTCCAGGTTTATCGCCTCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTTGCAGGGGCATG‑GTTCCCTTTATGACACCGGTTTTCCCGCAAAATTGGctcg                                >  2:364544/1‑248 (MQ=255)                                                                                                                                                                                                                 tAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCgg                              <  1:353904/250‑1 (MQ=255)                                                                                                                                                                                                                 tAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGGATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCggg                             >  1:137154/1‑250 (MQ=255)                                                                                                                                                                                                                  aa‑CCC‑TGGCGCTGGACGGCAGTAAAGGTCCAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGATGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCCGGGGCATG‑TTGCCCTTTTTGAACCCGGGTTTCCTGCAAA‑ATGGCCTCGGtt                            >  2:286762/1‑251 (MQ=255)                                                                                                                                                                                                                             gcTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGTCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGt                   >  1:181511/1‑251 (MQ=255)                                                                                                                                                                                                                                ggACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGGTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTc                >  1:89724/1‑251 (MQ=255)                                                                                                                                                                                                                                 gACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTGCCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCGTTATGACACCGTGTTGCCTGCAAA‑TTGGCCTCGGTTGCTCTCCGTTTCa               <  2:31742/251‑1 (MQ=255)                                                                                                                                                                                                                                  aCGGCAGTAAAGGTAACGGTTTATCGCCGCTGACGATTGCAGGCCGGGTCGCAGTTGGCGGGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTCAt              <  2:123884/251‑1 (MQ=255)                                                                                                                                                                                                                                   cGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTGCCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACCCCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTCATg             <  2:36456/251‑1 (MQ=255)                                                                                                                                                                                                                                            aGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCCGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCCTG‑GTGCCCTTTATTACACCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTCCTGCCGGACtg     >  1:265350/1‑250 (MQ=255)                                                                                                                                                                                                                                              gTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGCCGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTCATGCCGGACtgatg  >  1:244695/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            TTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAA‑CCC‑TGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGGGAGCTTTGCCGGACGTATTGGCAGTGGCATG‑TTGCCCTTTATGACACCGGTTTTCCTGCAAA‑TTGGCCTCGGTTTCTCGCCGTTTCATGCCGGACTGATG  >  NC_000913/2161884‑2162369

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