Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,839,359	0	T	G	32.4%	127.8 / 12.4	72	intergenic (‑95/‑165)	ascG/ascF	DNA‑binding transcriptional repressor AscG/beta‑glucoside specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base T (12/36);  new base G (23/0);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 3.15e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTTGGACGGTAACCGCTCTCTTCTACCGCCTGAAACACGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGT  >  NC_000913/2839120‑2839607                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                              aGTTTTGACGTTAACCGCTCTCTTCTTCCGCNTGAATCACGCGAGCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAGCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATTTTCAGGTGACCGGTTTCACAAata‑t                                                                                                                                                                                                                                                   <  2:321367/251‑1 (MQ=255)      ggACGGTAACCGCTCTCTTCTACCGCCTGAAACACGCGATCTTTAGTCTCCTGGCTGTCGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTGTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑Taaaaa                                                                                                                                                                                                                                              <  1:501910/251‑1 (MQ=255)              aCCGCTCTCTTCTACCGCCTGAAACACGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTGGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTGAGGGGGCCGGATTCACAAATA‑TAAAAAATGAACa                                                                                                                                                                                                                                       >  1:62358/1‑250 (MQ=255)               ccGCTCTCTTCTACCGCCTGAAACACGCGAGCTTTAGTCTCCTGGCTGTCGTAGCCATTACCTGAAAGCACCCTGGAAACGGTCGCTTTTGAAACCTCTTCGCGCTTCGCCACTTCCAGCAGCGTCGTCATCATTTTCATCCCTTTTCACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAAt                                                                                                                                                                                                                                     <  2:474627/251‑1 (MQ=255)                           aCCGCCTTAAACACACGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCtt                                                                                                                                                                                                                           <  1:302564/249‑1 (MQ=255)                           aCCGACTGAAACACGCGATCTTTAGTCTCCTCGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGAACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATAATCAAGAGACCGGGTACACAAATAATAAAAAA‑GAACAAtttaatctccttc                                                                                                                                                                                                                         >  2:158103/1‑240 (MQ=255)                                  gAAACACGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCACCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATA‑TAAAAAATGAAAAAATACCTCACTTGCTTAtttt                                                                                                                                                                                                                  >  2:230180/1‑250 (MQ=255)                                    aaCACGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGt                                                                                                                                                                                                                <  1:205040/251‑1 (MQ=255)                                          cgATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGGACGGTTTCACAAAAA‑TAAAAAATGAACAATTCACTCTTTTGCTTTTTTAGTG‑ACAAc                                                                                                                                                                                                         >  1:15661/1‑251 (MQ=255)                                           gATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGGCCGGGTTCACCAATA‑TAAAAAATGAAAAATTCACTCTATTGCTTATTTAGTG‑ACAAcc                                                                                                                                                                                                        >  1:46123/1‑250 (MQ=255)                                                   gTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑a                                                                                                                                                                                              <  2:15661/251‑1 (MQ=255)                                                    tctcCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGACCGGGTTCACAAAAA‑TAAAAAATGAACAATTAACTCTCTTGCTTATTTAGTG‑ACAAATA‑TTCATG‑At                                                                                                                                                                                             >  1:240297/1‑251 (MQ=255)                                                       cctggcgtaCGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑Atttt                                                                                                                                                                                          <  2:714547/243‑1 (MQ=255)                                                       ccTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑Att                                                                                                                                                                                            <  1:149774/249‑1 (MQ=255)                                                          ggCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTtgtg                                                                                                                                                                                       <  1:627860/251‑1 (MQ=255)                                                              gACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAAGGGACCCGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCCTATTTAGAG‑ACAACTA‑TTCATG‑ATTTTGTGaaa                                                                                                                                                                                    >  1:472511/1‑250 (MQ=255)                                                                   aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGtt                                                                                                                                                                              <  2:495020/251‑1 (MQ=255)                                                                    gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGttt                                                                                                                                                                             <  1:187278/251‑1 (MQ=255)                                                                      cATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTCCTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGACAAATATTCAGGTGCCCGGTTTTAACAATA‑TAAAAAATGAAAAAATCACTCTCTTTATTTTTTAGAG‑ACAACTA‑TTCCTG‑ATTTTGTTAAACCGGTTTCt                                                                                                                                                                           >  2:106602/1‑251 (MQ=255)                                                                         tACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGTCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTTACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTa                                                                                                                                                                         <  2:121977/250‑1 (MQ=255)                                                                           ccTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCTTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAAtt                                                                                                                                                                      <  2:440067/251‑1 (MQ=255)                                                                            cTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGGCCGGGTTCAAAAAAA‑AAAAAAATGAAAAATATACTCTCTGTCTTATTAAG‑G‑AAAAATACTTAATA‑ATTTTTTGAAACCGGGTTTTTAATTc                                                                                                                                                                     >  2:165809/1‑251 (MQ=255)                                                                              gAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTcc                                                                                                                                                                    <  2:543361/250‑1 (MQ=255)                                                                              gAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTc                                                                                                                                                                     <  1:407835/249‑1 (MQ=255)                                                                                 aGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCCATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGGGCGGGTTCCCCAAAA‑TAAAAAATGAACAATTTACTCTCTTGCTTATTTTGAG‑AAAAATA‑TTCATGCATTTT‑TGAAACCGGTTTCTTTATTCCGttt                                                                                                                                                                >  2:47496/1‑251 (MQ=255)                                                                                  gCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCTCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACTCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTc                                                                                                                                                               <  1:478424/251‑1 (MQ=255)                                                                                   cACGCGGGAAACGGTCGCGTTTGAAACCCCGGCGCGCTTCGCCACTTCCCTCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTTCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCa                                                                                                                                                              <  2:489447/251‑1 (MQ=255)                                                                                      gcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCa                                                                                                                                                           <  1:256708/251‑1 (MQ=255)                                                                                        gggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATc                                                                                                                                                         <  2:564808/251‑1 (MQ=255)                                                                                        gggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCAGGGTCAAAAAAA‑TAAAAAATGAACAATTTAATTCCTTTCTTATTTAATG‑GAAAAAA‑TTCAAG‑ATTTTGTGGAACCGGTTTCTTAATTCCgctttcgcagt                                                                                                                                                         >  1:381704/1‑241 (MQ=255)                                                                                         ggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGGCCGGGTTCACAAATA‑TAAAAAATGAACAATTCCCTCTCTTGCTTATTTTGTG‑ACAACTA‑TTCATG‑ATTTTTTGAAACCGGTTTCTTTATTCCGTTTCAGCATCg                                                                                                                                                        >  1:114786/1‑251 (MQ=255)                                                                                          gAAACGGTCGCTTTTGAAGCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAAATAAGGGGAACGGTTTTCCAAATA‑TAAAAAATGAAAAAtt                                                                                                                                                                                                                                    >  2:596334/1‑177 (MQ=255)                                                                                          gAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAAtt                                                                                                                                                                                                                                    <  1:596334/177‑1 (MQ=255)                                                                                                 tCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTCAGGGGGCCGGGTTCCCAAATA‑TAAAAAATGAAAAATTCACTCTCTTGCTTATTTTGAG‑ACAACTA‑TTCATG‑ATTTTTTAAAACCGGNTTCTTTATTTCGTTTCAGCCTCGGCCttttt                                                                                                                                                >  2:333463/1‑251 (MQ=255)                                                                                                  cGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCATCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTc                                                                                                                                               <  1:627236/251‑1 (MQ=255)                                                                                                    cTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGGGACCGGGTTCACAAATA‑TAAAAAATAAACAAATCACTCTCTTGCTTATTTAGTG‑ACAAATA‑TTCATG‑ATTTTGTGGAAACGGTTTCTTAATTCCGTTTTAGCATCGGCATTTTTCCg                                                                                                                                             >  1:136216/1‑251 (MQ=255)                                                                                                         gaatccccGGCGCGCTGCGCCACTTCCAGCAGCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCAc                                                                                                                                         <  2:71227/246‑1 (MQ=255)                                                                                                         gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACg                                                                                                                                        <  1:64498/251‑1 (MQ=255)                                                                                                         gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTGCCCCAAAA‑TAAAAAATGAACAATTTACTCCCTTGCTTTTTTAGTG‑AAAACTA‑TTCATG‑ATTTTGTGAAACAGGGTTTTTAATTCCGTTTCAGCATCGGCCTTTTTCCGTCACg                                                                                                                                        >  2:209031/1‑251 (MQ=255)                                                                                                          aaaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCCGGTGACCGGGATCCCAAATA‑TAAAAAAAGAAAAAATCAATCTCTTTCTTATTTAATG‑ACAACTA‑TTCATG‑ATTTTGTTAAACCGGTTTTTTAATTCCGGTTCAACAACGGCATTTTTTCGTCACgg                                                                                                                                       >  1:157597/1‑250 (MQ=255)                                                                                                                 ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGGCCGGTTGCAAAAATA‑TAAAAAATGAACAATTCCCTCTCTTGCCTATTTTGGG‑ACACCTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTTCGATTCCACCTCGGCATTTTTCCGTCACCTCGCCTg                                                                                                                                 >  1:413335/1‑250 (MQ=255)                                                                                                                  gcgcgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTTTTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCCTCGGCATTTTTCCGTCACGTCGACTGAt                                                                                                                               >  1:497299/1‑251 (MQ=255)                                                                                                                           gCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGGCCGGTATCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTTGTG‑ACAACTA‑TTCATG‑ATTTTGTGTAAACGGTTTCTTAATTCCGTTTCAGCATCGGCCTTTTTTCGTCACGTCGACTGATAACAACTa                                                                                                                       >  1:190143/1‑250 (MQ=255)                                                                                                                            ccACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGGGACCGGTTTCCCAAATA‑TAAAAAATGAAAAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑CTTTTTTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCgcc                                                                                                                                   >  2:302646/1‑235 (MQ=255)                                                                                                                            ccACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGAc                                                                                                                                   <  1:302646/237‑1 (MQ=255)                                                                                                                                 tCCAGCATCGTCGTCATCATTTTCATCCCTTTTCACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTa                                                                                                                 <  1:452669/250‑1 (MQ=255)                                                                                                                                 tCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTa                                                                                                                 <  1:99518/250‑1 (MQ=255)                                                                                                                                                atcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑Atttt                                                                                                                                                                                          <  2:93057/162‑1 (MQ=255)                                                                                                                                                atcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAg‑g                                                                                                <  1:536212/251‑1 (MQ=255)                                                                                                                                                atcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAAATA‑TTCATG‑Atttt                                                                                                                                                                                          >  1:93057/1‑162 (MQ=255)                                                                                                                                                  catTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGAAAAAAAATTCAGGGGTCCGGGTTCACAAATA‑TAAAAAATGAACAATTCATCCTCTTGTCTCTTTATTG‑ACAACCA‑TTCATG‑ATTTTGTGAAAACGGGTTTTTAATTCCGTTTCAGCATCGGCATTTTTCTGGCACGTTGACTGATAAAAACTACCACTACCCTCCTgcataaaggt                                                                                                >  1:704757/1‑241 (MQ=255)                                                                                                                                                      ttCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATA‑TAAAAAATGAACAAATTACTCTCTTGCTTATTTAGTG‑ACAAACA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTTATTCCGTTTCAGCATTGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATaaa                                                                                           >  1:538399/1‑250 (MQ=255)                                                                                                                                                        cATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGGTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACCACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCCTCGGCATTTTTCCGTCACGTCGACTGATAACAAATACATCTTACCTACTGATAACAG‑GATAAAATc                                                                                        >  1:669837/1‑251 (MQ=255)                                                                                                                                                        cATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCCGTTTCAAAAAAA‑TAAAAAATGAAAAATTCACCCTCTTTCTTATTTAGTG‑ACAACTA‑TTCCTG‑ATTTTGTGGAAACGGTTTCTTTATTCCGTTTCAGACTCGGCATTTTTTCGGCACGTCGACTGATAACAACTAAATCTACCCTACTGATAAAAG‑GATaaaaac                                                                                        >  1:452718/1‑249 (MQ=255)                                                                                                                                                         aTCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATcc                                                                                       <  1:434048/251‑1 (MQ=255)                                                                                                                                                                       aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTAAGGGGACCGGGTTCACAAATA‑AAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑GCAACTA‑TTAATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGCTTCAGCATCGGCATTTTTTCGTCACGTCGACTGTTAAAAACTCCATCTTCCCTACTTATAACAG‑GATAAAATCCGATGGCCAAAAAtt                                                                         >  2:205449/1‑251 (MQ=255)                                                                                                                                                                          cAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATg                                                                      <  1:230180/251‑1 (MQ=255)                                                                                                                                                                                      cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGTACGGGTTCACAAATA‑TAAAAAGAGAACAAATTACTCTCTTGGTTATTTAGTGGACA‑CCC‑TTCATG‑AATTTTTGCAACCGGTTTCTTTATTTCGTTTCAACAACGGCATTTTTTCCTCACCTCGACTGATAACAACATCATTTACCCTACTGATAACAGAGA‑AAAATTCGATTGCCAAAAATTTTGAGGCGCTGGCAc                                                          >  1:439901/1‑251 (MQ=255)                                                                                                                                                                                       tGCAACGTTTGCCGATTGTACTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCCCCAATA‑AAAAAAATGAACAAATCACCCTCTTGCTTATTTATTG‑ACAAATA‑TTCAAG‑ATTTTTTGCAACCGGTTTCTTAATTTCGTTTCAGTATCAGCATTTTTACATCCCGTCGCCTTATAACAAATACATCTACCCCACCGATAACAG‑GATAAAATCCAATGGCCCAAAAtcttacggccctggcccc                                                         >  2:81401/1‑234 (MQ=255)                                                                                                                                                                                          aCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATaac                                                                                                                            <  1:522251/182‑1 (MQ=255)                                                                                                                                                                                          aCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATA‑TAAAAAATGAAAAATTCACTCTCTTGCTTATTTTGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGAATGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCAATGGCCAAAAATTATTCGGCGCTGGCACGCTc                                                      >  1:126088/1‑251 (MQ=255)                                                                                                                                                                                          aCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGTCCGGGTTCACAAAAA‑TAAAAAAAGAAAAATTCACTCTCTTGCTTATTTATAT‑ACAACAA‑TTTATT‑ATTTTTTTAAACCGGTTTCTTAATTTCGTTTTAGCATCCGCATTTTTCCGTCACGTCGACTGATaac                                                                                                                            >  2:522251/1‑182 (MQ=255)                                                                                                                                                                                                 cccgccgattgtggCACAATCGGCGGGAAAAATATTCAGGGGACAGGGTTCACAAATA‑TAAAAAATGAAAAATTCACACTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTTTGAAACCGGGTTCTTAAATCCGGTTCAGCATCGGCATTTTTCCGTCACGTCGCATTATAACAACTACATCAACCCTACTGATAACAG‑GATAAAATCCGCTGTCCAAAAAATATTCGGCGCGGGCACGCTGGGTGAt                                                >  1:373492/14‑250 (MQ=255)                                                                                                                                                                                                   cGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGCCCGGGTTCACAAATA‑TAAAAAAACAACAATTTACTCTCTTGCTTATTTAGGG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCCTTTTTCCGTCACGTCGACTGATAACAACTACATCTCCCCCACTGTTAACAG‑GATAAAATCCGATGGCCAAAAAATATGCGGCGCCGGCACGCTCGGGGATAgg                                             >  2:591587/1‑250 (MQ=255)                                                                                                                                                                                                       tGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGc                                          <  1:18397/250‑1 (MQ=255)                                                                                                                                                                                                        gTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACTCTCGGTGATAGCGGCAc                                        <  1:471382/251‑1 (MQ=255)                                                                                                                                                                                                          ccTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACt                                       <  1:362276/250‑1 (MQ=255)                                                                                                                                                                                                          ccTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTg                                      <  1:401907/251‑1 (MQ=255)                                                                                                                                                                                                                         gggAAAAATATTCAGGTGACCGGTATCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAAc                       >  1:38340/1‑251 (MQ=255)                                                                                                                                                                                                                           gAAAAATATTCAGTTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCATCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATACCAt                     <  2:157597/251‑1 (MQ=255)                                                                                                                                                                                                                                        gTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACg        >  1:187876/1‑251 (MQ=255)                                                                                                                                                                                                                                               gTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTAACGCACTGt  >  1:353094/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                              AGTTTGGACGGTAACCGCTCTCTTCTACCGCCTGAAACACGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATA‑TAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTG‑ACAACTA‑TTCATG‑ATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAG‑GATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGT  >  NC_000913/2839120‑2839607

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