Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,826,478	T→C	18.3%	I30T (ATC→ACC)	srlE →	glucitol/sorbitol‑specific enzyme IIB component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,826,478	0	T	C	18.3%	111.0 / 2.0	59	I30T (ATC→ACC)	srlE	glucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base T (33/11);  new base C (2/10);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 3.98e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGCTGGTTGGTCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATT  >  NC_000913/2826266‑2826710                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             aTCTGCTGGTTGGTCTGGTCACCACTTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGCAACAGGAACGGGTGGCTGGGGC‑GGCCCGCGT‑GAGCTGAAAGACAC‑GCCGGGAAAAAAAA‑ACTACG‑TA‑AT‑CC‑CCCCCGGCACCCAGCCTGcccagt                                                                                                                                                                                                                  >  1:886104/1‑234 (MQ=255) aTCTGCTGGTTGGTCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGCTGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGTTGGGGC‑GGCCCGCTT‑GCGCTGCAAGCCAC‑GCCGGGAAAAAAAA‑TCTTCT‑AT‑TT‑TA‑CCGCCGGTACCCGGCCTCCGCttg                                                                                                                                                                                                                  >  2:772998/1‑239 (MQ=255) aTCTGCTGGTTGGTCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGAttg                                                                                                                                                                                                                  >  2:753159/1‑239 (MQ=255) aTCTGCTGGTTGGTCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGCCT‑AT‑AT‑CA‑CCCCCGGTACCCGGCCTGCGAttg                                                                                                                                                                                                                  >  2:1293415/1‑239 (MQ=255)       tggttggtCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GACGGGCAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACACGGCCTGCGATGGTTGACa                                                                                                                                                                                                            >  2:259429/1‑239 (MQ=255)       tggttggtCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGACTATTCGGATCGAAAAAGGCACGGCGGGCTGGGGC‑GGCCCGATT‑TGGCTGAAAGCCAA‑GCCGGGAAAAAAAA‑TCGTCT‑CT‑AT‑CC‑CAGCCGGTACCCGGCCTGCGATTGTTGaaa                                                                                                                                                                                                            >  1:1018594/1‑237 (MQ=255)         gttggtCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCACCTTGAACAAAAAGATCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAACAGGAACGGGTTGCTGGGGT‑GGCCCGCTT‑GGGCTAGAAGCCAC‑GCCGGGAAAAACAA‑TCGTGT‑AT‑TC‑AA‑CCGACGGGAACAGCCCTGCGAttggagaaaa                                                                                                                                                                                                           >  1:1094541/1‑231 (MQ=255)         gttggtCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACaaa                                                                                                                                                                                                          >  2:215560/1‑239 (MQ=255)         gttggtCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑CCGTCT‑AT‑AT‑CA‑CCGCCGGTCCCCGGCCTGCGATTGTTGACaa                                                                                                                                                                                                           >  2:792217/1‑238 (MQ=255)                     aCCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCACCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCtt                                                                                                                                                                                              >  2:388864/1‑239 (MQ=255)                                 cgcgGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACCGCTTActggctctg                                                                                                                                                                                    >  2:122944/1‑234 (MQ=255)                                   cgGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCCTATTCGGATCGAAAAAGGAACGGGTGGCTAGGGC‑GGCCCGCTT‑GAGCTGAAAGCCCC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCCCCGGTACCCGGCCTGCGATTGTTTACAAACGGGCACAGCTTACTCGGTGGCA‑GGc                                                                                                                                                                               >  1:930465/1‑239 (MQ=255)                                   cgGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCAAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AC‑AA‑CCGCCGGTACCCGCCCTCCGCTTGTTGACAAACCGGCACACCTAACTGGCGGGAA‑gac                                                                                                                                                                               >  1:336678/1‑237 (MQ=255)                                        gggTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCgg                                                                                                                                                                                                                            >  2:653671/1‑190 (MQ=255)                                        gggTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTAttt                                                                                                                                                                          >  2:1036269/1‑238 (MQ=255)                                                     gACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTTACGGATTTaaaaa                                                                                                                                                             >  2:243546/1‑237 (MQ=255)                                                     gACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTaaa                                                                                                                                                               >  2:124588/1‑237 (MQ=255)                                                           cGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGACCGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGCTT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTGTGCAAACTGGCACAGCTTCATGGCTGGCA‑GGCGATTGACGGATTTAAAAA‑Agggag                                                                                                                                                      >  1:904957/1‑236 (MQ=255)                                                           cGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGGTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCCA‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCCGCGATGTTGGACAAAATGGCACAGCTTACTGGC‑GGCAGGGCTAGTGAGGGGTTTAAAGA‑AGGGGa                                                                                                                                                      >  1:862541/1‑239 (MQ=255)                                                                tttttGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAAccc                                                                                                                                                  >  2:164666/1‑238 (MQ=255)                                                                  ttgaaaaaaagGAGGGGATTTAAACTTAAAAAAAAATTTCACCGGGCAGGAACAACATCATG‑CCGCATATCGGGAGCAAAAAAGGAGCGGGTGGCTGGGGG‑GGCCCGCTG‑GAGCTAAAAGCCAC‑GCGGGAAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCcgcg                                                                                                                                               <  2:475857/233‑1 (MQ=255)                                                                        aaaaaGATGGGCATTCAACTTGAACAAAAAGTTCACCTTGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGACCC‑GCCGGGAAAAAAAA‑TCGTCT‑CT‑AT‑TA‑CCGCGGGGACCCGGCCTGCGATTGTTGCGAAATGTGAAAAGATTACTGGCTGGCC‑GGCGATTGACGGCGTTTAAGACAGGT‑AACCCACGGAGGCg                                                                                                                                         >  1:781315/1‑239 (MQ=255)                                                                                       cAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCCCGGAGGCGGAAATTGGTGTcgcg                                                                                                                          >  2:1198450/1‑239 (MQ=255)                                                                                              aaCAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGCCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCTCGGAGGCGGAAATTGGTGCCCCGGAAATc                                                                                                                    >  2:741194/1‑238 (MQ=255)                                                                                                        tCACCTGGGAGGGGCAAAAT‑ATGTGCGCATAGTCGGAAACAAAAAAGGACGGGGGGGGGGGGCGGGCCCGTTT‑GGGCTAAAACCCCC‑GCGGGAAAAAAAAA‑CCGTCT‑AT‑AA‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTggcg                                                                                                          <  2:336678/239‑1 (MQ=255)                                                                                                               gCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGCAATTGGTGTCGCGGTAATCGACTGTGGCGGCCCAtt                                                                                                   >  2:291140/1‑238 (MQ=255)                                                                                                                   gAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCCCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACCTTAc                                                                                                 >  2:523169/1‑236 (MQ=255)                                                                                                                       caaaatcaacgcgaattgtCGGACAGAAAAGGGAACGGGTCGGTGGGGC‑GCCCCGTTG‑GGGCAGAAAGCCAC‑GCCGGAAAAAAAAC‑TCGTCT‑AT‑AC‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGc                                                                                            <  2:1018594/220‑1 (MQ=255)                                                                                                                         catcatG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGGTTTAAAAA‑AGGGGAACCCGCGGAGGCGGGAAATGGTGGCGGGGTAAACCACTGGGGGGGGCCAATACGcgtgtg                                                                                          >  1:1298691/1‑232 (MQ=255)                                                                                                                            atcggagcaattttgcggaaaaaaaaGGACCGGGGGGGGTGGGG‑GGCCCGCTTGGAGGTAAAACCCCC‑CCCGGAAAAAAAAC‑CCGTCA‑AA‑AC‑CC‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCAt                                                                                       <  2:1136274/218‑1 (MQ=255)                                                                                                                              tG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGTCGGCCA‑GGCTATTGACGGATTTAAAGA‑AGGTGACACCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGGGGCGGGCAATTACggcggggacactc                                                                                    >  1:905400/1‑226 (MQ=255)                                                                                                                              tG‑ACGCATATTCGCATCGAAAACGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGCAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGCTTGTTGACAAACTGGCACAGCTTACTGGCTGCAA‑GGCTGTTGACGGTTTTAAGAA‑AGGTCAACCCCCGCAGCGGGAAATGCGGCTCGCGGGACTCGACGTTggcttccccggtcaccatggtgtcgc                                                                                    >  1:264796/1‑215 (MQ=255)                                                                                                                              gtgccGCATATGCGGACAAAAAAAGGAAGGGGGCGCGGGGGC‑GGCCGGTTT‑GGCCGAAAACCCAC‑GCCGGAAAAAAAAA‑TCTTCT‑AT‑AC‑CC‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCa                                                                                        <  2:997885/231‑1 (MQ=255)                                                                                                                                  cGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCACCTATCCg                                                                                >  2:1061930/1‑239 (MQ=255)                                                                                                                                     atatTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGCGAACCCGCGGGGGCGGAAATTGGTGTCGCGGAAATCGACTGTGGAGGCACATTaactcggggctttcatccgaa                                                                              >  2:49805/1‑218 (MQ=255)                                                                                                                                             aTCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACCTTACCCTGGGGCATCTATCCGGAacgaaggtt                                                                     >  2:1263285/1‑234 (MQ=255)                                                                                                                                             aTCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAAA‑AGGTGAACCCCCGGAGGCGGAAGTTGGGGTCGCGGGATACAACTGTGGCGGCACCTTACACTGCGGGACATGTCCGAAAAGAGGTa                                                                      >  1:935591/1‑238 (MQ=255)                                                                                                                                                aaaagaggggcGGGTGGGGGGGGG‑GGCCCGTTG‑GAGCTAAACGCCGC‑CGGGGAAAAAAAAA‑TCGTAT‑AT‑AT‑GC‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTcc                                                                  <  1:523169/229‑1 (MQ=255)                                                                                                                                                aaaaaagacaggggggtgtGGGGC‑GCCCCGCTG‑GGTTGGAAAGCCCA‑GACGGAAAAAAAAG‑TCTTGT‑AT‑AT‑CG‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACGGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTcc                                                                  <  1:1293415/221‑1 (MQ=255)                                                                                                                                                  aagagaaggggggggGTGGGGG‑GGCCCGTTT‑G‑GGTGAAAGCCCCGGCCGGAAAAAAAAC‑TCGTCA‑AA‑AC‑CG‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCcca                                                                <  2:876571/230‑1 (MQ=255)                                                                                                                                                        aagggggggggggggC‑AGCCGGTTG‑GAGCTGAAATCCAA‑GCGGGAAAAAAAAA‑TCGTTT‑AT‑CT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCAtt                                                           <  2:535930/227‑1 (MQ=255)                                                                                                                                                                            ccGCTT‑GAGATGAACCGCCC‑GCCGGAAAAAAAAA‑ACTTCT‑ATCCC‑CC‑CCGCCGGTCCCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGCCA‑GGCTATGGACGGATTTAAAGA‑AGGTGATCCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCAGCTAGCCGAAACGACGTATTCCCACCATTAATATCCTCGcggcggct                                         <  1:861876/238‑3 (MQ=255)                                                                                                                                                                                     gaagtaaaacccaaggcggAAAAAAAA‑GCGTTTGCT‑AC‑CG‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAgtccgg                                 <  2:656473/222‑1 (MQ=255)                                                                                                                                                                                     agaacaaacaccagtcgaaaaaaaaga‑TCGTCT‑AT‑AC‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAgtccgg                                 <  2:160211/219‑1 (MQ=255)                                                                                                                                                                                         aacgCCAC‑GCGGGAAAAAAAAC‑GCGTAT‑AA‑TA‑CC‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGc                            <  2:170876/235‑1 (MQ=255)                                                                                                                                                                                         aaccgcccccgataaaaaaaaaa‑TCCT‑T‑AT‑AG‑CCCCCGCCGGTTCCCGGGCTGGGAGTGGTGACAAACTGGCCAAGCTTACGTGCTGGCA‑GGCTATTGACGGATTAAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGc                            <  1:164666/224‑1 (MQ=255)                                                                                                                                                                                         aaaaCCAC‑GGGGGGAAAAAAAG‑CCGT‑T‑AA‑ATCCA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGc                            <  2:1187984/238‑1 (MQ=255)                                                                                                                                                                                          aacgccggcgcGGAAAAAAACA‑TCGTAC‑AA‑CC‑CA‑CCGCCGGTACCCGGCCTGCGATTGTCGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGTTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGc                            <  2:346967/228‑1 (MQ=255)                                                                                                                                                                                          aaGCCCC‑GCCGGGAAAAAAAACTCTTTT‑AT‑AT‑CT‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACTGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGc                            <  2:260089/238‑1 (MQ=255)                                                                                                                                                                                           aGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCCTTAATATCCACTCGACGGGCCAgtccggtccccgt                          >  2:861876/1‑234 (MQ=255)                                                                                                                                                                                              ccggcggagaaaaaaaac‑TCGTAT‑AT‑AT‑CG‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGgcg                       <  2:714842/229‑1 (MQ=255)                                                                                                                                                                                              cAC‑GCGGGAAAAAAAAA‑TTGTCT‑AC‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATGGGTGACAAACTGGCACAGCTTACTGGCGGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAAGCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGgcg                       <  2:691661/238‑1 (MQ=255)                                                                                                                                                                                               cccggtggaaaaaaata‑TCGTTT‑TA‑AT‑CC‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGGCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGgcgc                      <  2:644007/230‑1 (MQ=255)                                                                                                                                                                                                cgcggggaaaaaaaaa‑GCGTAT‑AT‑AC‑CG‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCa                     <  2:802411/231‑1 (MQ=255)                                                                                                                                                                                                 ccgggggAAAAAAAA‑CCGGCTAAT‑AT‑CC‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCa                     <  2:1204075/235‑1 (MQ=255)                                                                                                                                                                                                           aaaaa‑TCGTCT‑AT‑AC‑CA‑CCGCCGGTACCCGGTCTGGGATTGGTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTgg          <  2:405075/239‑1 (MQ=255)                                                                                                                                                                                                             aaatctctCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGa         <  2:862541/232‑1 (MQ=255)                                                                                                                                                                                                                   gTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCCACGGGCAAGACCGGTCCGCGGGCGCAGTCACTTTTGGAAGatat   >  2:284836/1‑239 (MQ=255)                                                                                                                                                                                                                    tCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCCCAGTAAATTTTGGAAGAAAtt  >  2:351698/1‑239 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             ATCTGCTGGTTGGTCTGGTCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATG‑ACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGC‑GGCCCGCTT‑GAGCTGAAAGCCAC‑GCCGGGAAAAAAAA‑TCGTCT‑AT‑AT‑CA‑CCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCA‑GGCTATTGACGGATTTAAAGA‑AGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGGCACATTACGCTGCGGCATCTATCCGAAACGACGTATTCCCACCATTAATATCCACTCGACGGGCAAGTCCGGTCCGCTGGCGCAGTACATTGTGGAAGATATT  >  NC_000913/2826266‑2826710

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