Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,136,467	0	A	C	35.1%	88.3 / 13.3	60	E120D (GAA→GAC)	yijE	cystine exporter
Reads supporting (aligned to +/- strand):  ref base A (15/22);  new base C (20/0);  total (38/22)
Fisher's exact test for biased strand distribution p-value = 2.83e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAGTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGGC  >  NC_000913/4136172‑4136763                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                         ggAGTTATAGCGGTATTGTCATGAAGCAAGTCACCATTTACATCGGGACCTTCGACTTTACCGCCTTACGCTGCATTTTTGTCGCTCTCGTTTTATTCATCGTCCCTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACNTG‑CGGGATGGTTGGTCTGGCGCAGTGGACGTTGGTCAGTGTAGGTGCGGGGGAGGTGGCGATCCTGAGCTATACTATGCCGTTCTGGGTGGGNATTGTCGCCGCGTTGTTTCTCGGTGAACGcc                                                                                                                                                                                                                                                                                                      <  2:599607/300‑1 (MQ=255)  gAGTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTAGCCGCGTTGTTTCTCGGTGAACGCCt                                                                                                                                                                                                                                                                                                     >  2:547054/1‑300 (MQ=255)    gTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGc                                                                                                                                                                                                                                                                                                   >  2:802730/1‑300 (MQ=255)                                     gTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCg                                                                                                                                                                                                                                                                  >  2:87039/1‑300 (MQ=255)                                        aCATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCACTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGTCGCtt                                                                                                                                                                                                                                                               >  2:123014/1‑300 (MQ=255)                                                gCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTa                                                                                                                                                                                                                                                       >  1:727739/1‑300 (MQ=255)                                                         tttACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGGGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATGTTTCGCGATTCTGATTGCCGCTTTCGGCTTATTGTTGgtg                                                                                                                                                                                                                                              >  2:717018/1‑300 (MQ=255)                                                                           aTTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTc                                                                                                                                                                                                                            <  1:730318/300‑1 (MQ=255)                                                                                         cGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCNCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGAt                                                                                                                                                                                                              >  1:601394/1‑300 (MQ=255)                                                                                                       tCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCGGTTACAAACCTG‑CGGGATGGTTAGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGTAACTCGATTTCTCTTCGATGAAAAGTGCCATGc                                                                                                                                                                                                <  1:204823/300‑1 (MQ=255)                                                                                                           tttATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGc                                                                                                                                                                                            >  2:521003/1‑300 (MQ=255)                                                                                                              attaCGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAAt                                                                                                                                                                                         <  1:654178/300‑1 (MQ=255)                                                                                                                taCGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATcc                                                                                                                                                                                       <  1:799824/300‑1 (MQ=255)                                                                                                                      gTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCTACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCCTGCTGGCAATCctctacg                                                                                                                                                                                 >  2:107416/1‑297 (MQ=255)                                                                                                                       cctctgactacacccaacaccgggtAAATACACCTCATCCATTGACCGGGTACAC‑CCTGGCGGGATGGGTAGTCGGGCGCAGTGGGTGTGGGTCAGCGTAGGTGCGGGGAGGGTGGCGATCCTGATCTGTACCATGTTGTTTTGGGTGNTGATTTTCGCCGCGTTGCTTCTCGGTNAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGTTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCgg                                                                                                                                                                                <  2:562191/276‑1 (MQ=255)                                                                                                                                  gcCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTGCGAAGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCAGTGGCAACTAGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTggg                                                                                                                                                                     >  2:148536/1‑300 (MQ=255)                                                                                                                                   cccGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTgggg                                                                                                                                                                    <  1:649696/300‑1 (MQ=255)                                                                                                                                            gTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCAGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAgcgc                                                                                                                                                           >  2:749124/1‑300 (MQ=255)                                                                                                                                                aaaTACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGAtt                                                                                                                                                       >  2:472824/1‑300 (MQ=255)                                                                                                                                                     cacCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGGTTGTTGc                                                                                                                                                  >  2:405501/1‑300 (MQ=255)                                                                                                                                                                 tGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTa                                                                                                                                      <  1:451963/300‑1 (MQ=255)                                                                                                                                                                     cTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGcc                                                                                                                                  <  1:171304/300‑1 (MQ=255)                                                                                                                                                                       gTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTGCGAAGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCAATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGccag                                                                                                                                >  2:148511/1‑298 (MQ=255)                                                                                                                                                                                                      gcAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGANGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAAc                                                                                                  <  1:589770/300‑1 (MQ=255)                                                                                                                                                                                                               gTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAANATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCa                                                                                         <  1:616706/300‑1 (MQ=255)                                                                                                                                                                                                                     ccgCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCAGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCt                                                                                   >  1:523308/3‑300 (MQ=255)                                                                                                                                                                                                                            ggTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTTCTCATTGAGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGTGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATgcg                                                                            >  2:505972/1‑300 (MQ=255)                                                                                                                                                                                                                             gTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCGTGGGCAATGTTTCGCGATTCTGATTGCAGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGCTGAAAAGTGCCATGCTGGCCATCCTCTCCGGCGTCGGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCNGTCATCCGAGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATgcgg                                                                           >  2:627502/1‑300 (MQ=255)                                                                                                                                                                                                                              tGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGTTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATgcggc                                                                          <  1:308847/300‑1 (MQ=255)                                                                                                                                                                                                                              tGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATgcggc                                                                          <  1:348400/300‑1 (MQ=255)                                                                                                                                                                                                                                      aGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTGCGCCTTAGGCAATATTTCGCGGTTCTGATTGCTGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTNCCATGCTGGCAATTCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGtga                                                                  >  2:202898/1‑300 (MQ=255)                                                                                                                                                                                                                                          ggCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCNGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAg                                                              <  1:580788/300‑1 (MQ=255)                                                                                                                                                                                                                                              aTCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCTTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGtgtg                                                          >  2:712997/1‑300 (MQ=255)                                                                                                                                                                                                                                              aTCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGtgtg                                                          >  1:258992/1‑300 (MQ=255)                                                                                                                                                                                                                                              aTCCTGAGCTATACCATGCAGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGATTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCAGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGtgtg                                                          >  1:658252/1‑300 (MQ=255)                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTTCTCCTTGGGCAGTATTTCGCGATTCTGATTGCCGCTTTAGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTAGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTAAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGTAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  2:457695/1‑300 (MQ=255)                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGATCGCCTTCTCCTTAAACAATATTTAGCGATTCTGATTGCCGCTTTCGGATTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTATTAGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTACGCGGCGCTGGTGATGAGTGTGGtagc                                                     >  2:313123/1‑297 (MQ=255)                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTTCTCATTAGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGACGTGGAAACTCGATTTCTCTTCGATGCAGAGTGCAATGCTGGCAATCCTCTAAGGCGTAAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTATATGCCCGTCATCCGCGCGTGGATTTATTGTCGCTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  2:222612/1‑300 (MQ=255)                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  1:75803/1‑300 (MQ=255)                                                                                                                                                                                                                                                     gCTATACCATGCCGTTCTGGGTGGTGATTTTNGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCGTGGGCAATATTTCGAGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGNCAACTCGATTTCTCTTCGATGAAAAGGGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCACGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCtt                                                   >  2:587930/1‑300 (MQ=255)                                                                                                                                                                                                                                                     gCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGATCGCCTTCTCCTTGAACAATGTTTAGAGATGATGAGTGCAGATTTCGGATTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGTGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTAGTTAACATCCCGGCAGATGCTGTATGCGGCGCTGGTGATGAGGGTGGTCGCtt                                                   >  2:135650/1‑300 (MQ=255)                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCNGCGTTGTTTCTCGGTGAGCGCCTGCGAAGTGGGAAAGGTTTAGAGATTATGATTGCAGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCANCTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCGTCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTAc                                                >  2:585451/1‑300 (MQ=255)                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTTCGCAGTGAGCAATATGTCGAGAGTATGATTGCCGATTTCGGCTTACTTTTGGTGTTGCAGTCGTGGCAACTAGATTTCTATTAGATGAAAAGTGCAATGCTGGCAATTCTATACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCACGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATACGGCGCTGGTGATGAGTATGGTCGCTTTAc                                                >  2:86131/1‑300 (MQ=255)                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATGTTTAGCGATTTTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGCTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTAc                                                >  2:791324/1‑300 (MQ=255)                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTAc                                                >  1:410406/1‑300 (MQ=255)                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTAGGTGACCGCCTGCTCCTTAAGCAAGAGGTAGGGAGGCTGAGTGCTGTTTTCGGTTTATTTTTGGTGTTGCAGAAGTGGCAACTCGATTTCTATTAGCTGAAGAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGTGGTTGCTGCTAAACATCTGTATGCACGTAATGAGCGTGGGGATTTCGTGTAGTTACGATCGTGGGAGGTGCTGTATGAGGCGCTGATGATGAGTGTGGTCGCGTTAc                                                >  2:505503/1‑300 (MQ=255)                                                                                                                                                                                                                                                           ccATGCCGTTCTGGTTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTNCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTgg                                             <  1:584227/300‑1 (MQ=255)                                                                                                                                                                                                                                                           ccATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTGACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTgg                                             <  1:335043/300‑1 (MQ=255)                                                                                                                                                                                                                                                                  gTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTTCGCAGTGGGCAATATTTCGAGGTTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGAGTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCa                                      >  2:179677/1‑300 (MQ=255)                                                                                                                                                                                                                                                                  gTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCTCCGTGGGCAATGTTTAGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTATTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCa                                      >  2:14507/1‑300 (MQ=255)                                                                                                                                                                                                                                                                  gTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTCCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCa                                      >  1:132781/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           gttgATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAAt                             <  1:465332/298‑1 (MQ=255)                                                                                                                                                                                                                                                                           ggtgATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAAt                             <  1:847555/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                     gccgcGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAg                   <  1:627696/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                            tgttTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTGTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCg            <  2:531520/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             gttTCTCGGTGAACGCCTGCGACGTGGGCCATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCgt           >  1:28961/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                     gTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTggg   >  2:383582/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                     gTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCAATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGNATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTggg   >  2:190521/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                      tGAACGCCTGCGACGTGTGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGgc  <  1:197888/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      tGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGgc  <  1:392317/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                         GGAGTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTG‑CGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGGC  >  NC_000913/4136172‑4136763

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