Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	41.8%	93.6 / 34.3	67	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (34/5);  new base A (1/27);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 2.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTG  >  NC_000913/1755722‑1756016                                                                                                                                                      |                                                                                                                                                  aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGt                                                                                                                                                   >  1:798772/1‑150 (MQ=255)       ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccg                                                                                                                                            <  1:894622/151‑1 (MQ=255)         aCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgcc                                                                                                                                          <  1:870677/151‑1 (MQ=255)             aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGAAAAACcgccgc                                                                                                                                        <  1:322754/149‑1 (MQ=255)                  cGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTg                                                                                                                                 >  1:25338/1‑151 (MQ=255)                     aTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCtt                                                                                                                              >  1:273517/1‑151 (MQ=255)                        tGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGAt                                                                                                                           >  1:700373/1‑151 (MQ=255)                          aaGAGATGTTACAAAAAATGCGGGACGCGGGCAGGAAGGTTGTGCGTCTGAACTTCTCTCGGGGTGATTATGCAAAACAGGGTCAGCCCATCCAGAATCTGCGAATGGTGATGAAAAAAAGGGAAAAAACCGCCGCTATCCTGCTTGAt                                                                                                                           <  1:816608/149‑1 (MQ=255)                                ttttAGCTAAAATGCTGGGCGCTGGCATAAACGTTATGCGTCTAAACTTCTCCCATGGTGACTATGCAAAACACGGTCGGCGCATTCAGAATCTGCGCAACGTGATGACAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                                        <  2:1015659/144‑1 (MQ=255)                                tGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                   >  2:977041/1‑151 (MQ=255)                                tGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                   >  1:206897/1‑151 (MQ=255)                                  ttagtaaaaaTGCTGGGCGTTGGCATGACCGTTATGCGTATGACCTTTTCTCATGTTGACTATCAAAACCCCGGTCAGCGCATGAAGATTCTGCGCAAGGTAATAACCAAAGTTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGt                                                                                                                  <  1:896276/144‑1 (MQ=255)                                  ttagaaaaatTGCGGGACGCTGGCATAAACTTTATGCGTCTGAACTTCTCTCATGGTTACTAGGCAGAACACGGTCAGCGCATTCAAAATCTGCCCAACGTGAAGAAAAAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                                      <  2:916712/140‑1 (MQ=255)                                  ttAGCAAAAATGCGGGACGCGGCAATGACCGTTATGCTTCTAATCTTCCCTCATGGTGACTATGAAAAACACGGTCAGCGCATTAAAAATCTGCGCAAGGTGATGAAAAAAATTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                   <  2:679267/149‑1 (MQ=255)                                    aacaaaaaTGCTGGCCGCTGGCATGACCTTTAGGCGTCTAAACTTCTCTCATGGTGATTATGCAAAACACGGCCAGCGCATCCAGAATCGGCGCAACGTAATGAGAAAACCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTcc                                                                                                                <  1:527507/146‑1 (MQ=255)                                      cTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCgg                                                                                                              >  2:634785/1‑150 (MQ=255)                                       tAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                             >  2:36991/1‑150 (MQ=255)                                        aaaaTGCTGGCCGCTGGCATGAACGTTATGCGTCTGACCTTCTCTCAGGGTGACTATGCAGAACCCGGTCAGCGCATTCAGAATCTGCGCAACGTGATGACAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCgg                                                                                                              <  2:906714/148‑1 (MQ=255)                                        aaaaTGCTGGACGCTGGCATGAACGTTATGCGTCTGACCTTCTCTCAGGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCGGCGCAACGTGATGAGAAAACCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGaa                                                                                                            <  2:1019320/150‑1 (MQ=255)                                        aaaaGGCGGGCCCCTGGCATAACCTTTAGGCGTCTGACCTTCTCTCGTGTTAATTAGCCAAAACTCGCCCGGCAAATCAAGACCCGGCGAACGGTGATGATAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                             <  2:460857/149‑1 (MQ=255)                                          aaTGCGGGACGCGGGCATGAACTTTCTGCGTCTAAACTTCTCTCATGGTGACTATGCAAAACAGGGCCAGCGCATTCAAAATCTGCGCAACGTGATAAAAAAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                                         <  2:374077/151‑1 (MQ=255)                                               ggggcgcgGGCATGACGGTTATGCGTCTAAACTTCTCTCATGGTGATTATGAAGACCACGGTCGGCGCATACAAAATCTCCGCAACGTGATGAGAAAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                    <  1:766938/147‑1 (MQ=255)                                                ggACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAcc                                                                                                   >  2:932052/1‑151 (MQ=255)                                                       ggCATAAACTTTATGCGTATGACCTTCTCTCAGGGTGACTATGAAAAACAGGGTCAGCCCTTCAAAAATCTGCGAAGGGTGATGAGAAAATCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAAc                                                                                            <  2:196242/151‑1 (MQ=255)                                                             aaCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                        >  2:714654/1‑149 (MQ=255)                                                                 ttATGCGTCTGACCTTCTCTATTGGTATCTATGCAAAACACGGGCGGCGCATCAAGACTCTGGGCACCGTGATGAACAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                    <  2:1045483/149‑1 (MQ=255)                                                                    tGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTa                                                                                >  2:886116/1‑150 (MQ=255)                                                                      cGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                             >  1:467540/1‑151 (MQ=255)                                                                       gtttGAACTTCTCCGATGGTGACTAGGAAAAACGCGGACAGCCCATTCAGAATCGGCGCAACGTGATGAGAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAac                                                                              <  1:230979/146‑1 (MQ=255)                                                                        tCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                           >  1:709223/1‑151 (MQ=255)                                                                        tCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                           >  2:702000/1‑151 (MQ=255)                                                                             aCTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtc                                                                      >  2:456984/1‑151 (MQ=255)                                                                                 ctctcATGGTGACTATGCAAAACACGGCCAGCGCTTTAAAAATCGGCGCAACGTGATGAGAAAAAGTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                    <  2:420561/149‑1 (MQ=255)                                                                                      aTGGAGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                             >  1:26673/1‑151 (MQ=255)                                                                                       tGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                            >  2:434286/1‑151 (MQ=255)                                                                                        gtgaattATGCAAAACTCGGTCGGCCAATTCAAAACCCGCGAAAGTTGATAGGCAAAAGGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                                <  1:434086/140‑1 (MQ=255)                                                                                        ggTGACTATGCAGACCACGGCCAGCGCATTCAAAATCGGCGAAACGTGAAGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAg                                                              <  2:987349/148‑1 (MQ=255)                                                                                                tGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAcc                                                   >  2:541999/1‑151 (MQ=255)                                                                                                 gCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                  >  1:955845/1‑151 (MQ=255)                                                                                                  cAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                                 >  2:301267/1‑151 (MQ=255)                                                                                                     aaCGCGGTCAGCCCATCCAAAATCGGCGCAACGTGAAGAGAAAACCGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                                <  1:91609/149‑1 (MQ=255)                                                                                                         cGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                          >  1:681623/1‑151 (MQ=255)                                                                                                              aGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                     >  1:102404/1‑151 (MQ=255)                                                                                                                cgcATTCAAAATCTGCGCAACGTGATGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                     <  2:789781/149‑1 (MQ=255)                                                                                                                  cATTCAGAGTCTGCGCAACGGGATGAGCAAAACTGGAAAAACCGCCGCGATCCTGCTTGATACCGAAGGGGAGGAAATCCGCACCATGAAACTGGAATGCGGTAACGACGTATCTCTGAACGCTGGTCAGACCTTTACTTTCACCACTGAt                                 >  1:505685/1‑151 (MQ=255)                                                                                                                  cATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                 >  2:91351/1‑151 (MQ=255)                                                                                                                    ttCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATa                                >  2:41980/1‑150 (MQ=255)                                                                                                                        gAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                           >  1:892426/1‑151 (MQ=255)                                                                                                                         actCTGCGCAAGGTGAAGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                           <  1:853321/148‑1 (MQ=255)                                                                                                                           tCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                        >  2:202983/1‑151 (MQ=255)                                                                                                                              gcgcAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATc                     >  2:49939/1‑151 (MQ=255)                                                                                                                                  agcGTGATGAAAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                  <  1:476701/148‑1 (MQ=255)                                                                                                                                  aaCGTGATGAGAAAAACTGATAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCa                 <  2:838708/151‑1 (MQ=255)                                                                                                                                   aCGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCa                 >  1:49178/1‑150 (MQ=255)                                                                                                                                   aCGTGATAAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                        <  2:671852/143‑1 (MQ=255)                                                                                                                                    cgagaaaaaaaaacCTGGAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                      <  2:25339/134‑1 (MQ=255)                                                                                                                                    cGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc               >  2:292033/1‑151 (MQ=255)                                                                                                                                     gTGATGAAAAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc               <  2:758807/150‑1 (MQ=255)                                                                                                                                        attaGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg           >  2:388966/4‑151 (MQ=255)                                                                                                                                        atgaGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg           >  1:55957/1‑151 (MQ=255)                                                                                                                                         tgacaaaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc            <  2:750544/144‑1 (MQ=255)                                                                                                                                         tgaGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg           >  1:471577/1‑150 (MQ=255)                                                                                                                                           agaaaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa          <  1:647233/146‑1 (MQ=255)                                                                                                                                              aaaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg     >  1:177391/1‑151 (MQ=255)                                                                                                                                              aaaaCGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg      <  2:771404/150‑1 (MQ=255)                                                                                                                                               aaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGt    >  1:507404/1‑151 (MQ=255)                                                                                                                                               aaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGt    <  2:206933/151‑1 (MQ=255)                                                                                                                                                 aCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt   >  2:316497/1‑150 (MQ=255)                                                                                                                                                     ttAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTg  <  2:519919/146‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTG  >  NC_000913/1755722‑1756016

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