Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,502,757	C→T	18.5%	E226K (GAA→AAA)	insG ←	IS4 transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,502,757	0	C	T	18.5%	124.5 / 5.4	62	E226K (GAA→AAA)	insG	IS4 transposase
Reads supporting (aligned to +/- strand):  ref base C (4/46);  new base T (10/2);  total (14/48)
Fisher's exact test for biased strand distribution p-value = 5.31e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGA  >  NC_000913/4502514‑4502987                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            ccacaTCCCCCTTTTCCCCTCCGGGGAAGCGCAGGGGGCCCGTCACCGACGCCAGCAGTTGGAAGCTTTTTCC‑TTTGCCGGTCCCGGTCACCAGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACC‑‑‑‑‑‑‑CTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑Ttctc                                                                                                                                                                                                                                           <  2:610719/235‑1 (MQ=255)        ccccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCTTGCCTGCGGGCTGGTTTTCAGCTTACCCCGATGATCGCCTTTACCCAGTTTTCTGGTCTCTTCATATT‑GCCCTCCCTTCCTCGGGGGATTCATCCAGTGGCGGGT‑Ttct                                                                                                                                                                                                                                            >  2:81351/1‑238 (MQ=255)        ccccATTTCCCCACCGGGGAAGCGCAGGGCGTCCGTCATGGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑Ttctc                                                                                                                                                                                                                                           <  1:292430/239‑1 (MQ=255)          ccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTcc                                                                                                                                                                                                                                          >  1:117236/1‑238 (MQ=255)          ccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACC‑‑‑‑‑‑‑CTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGG‑GTTTTCCGCCCgg                                                                                                                                                                                                                                   >  1:610719/1‑238 (MQ=255)              ttCCCCCCCGGGGAAGCCCATGGCGGCCGTCAGCGGCGTCAGTAGATGGCAGTCTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCa                                                                                                                                                                                                                                     <  1:724208/239‑1 (MQ=255)               tctcCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGCTGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑TCGCTCCCTTTCTGAGAGGTATCATCCAGTGGGGGCG‑GTTTCTCGcccc                                                                                                                                                                                                                                    >  2:750506/1‑237 (MQ=255)                    tccGGGGAAGCGAGTGGCGTACGTCGACGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTAACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCCCTTTTTGCGGGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTcc                                                                                                                                                                                                                               <  1:628689/239‑1 (MQ=255)                    cccGGGGAAGCGCATGGCCGCCGTCATCGCCGGCATCAGATGGCAGACTTTTCC‑TTTGCGCGTCCCGGTCAGCGGGCGGGCAGGCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTcc                                                                                                                                                                                                                               <  2:824813/238‑1 (MQ=255)                             gcgcATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATAATCCAGTGGCGGGGTTTTTCCCCCCAgggcccaaggatta                                                                                                                                                                                                                       >  2:415310/1‑227 (MQ=255)                              cccctGGCGTCCGTTCTCGGCGCCAGCAGATGGCCGCTTTTTCC‑TTTGCGCGTCCCGGTCCGCCGGCGGGCAGTCATTTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                     <  2:395462/235‑1 (MQ=255)                               gcgtGGCGCCGGTCTTCGTCGTCAGCAGATGGCAGACTGTGCC‑TTTGTGCGTCACGGTCAGCTTGCGGGCAGTCACTTCATTCCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGTGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                     <  1:351081/235‑1 (MQ=255)                               ccATGGGGTCCGGCATCGACGTTCGCAGTTTGCAAACTTTTCC‑GTGGCGCGTCACGGTCAGGAGGCGGGGCGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTCCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTTTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                     <  2:30030/237‑1 (MQ=255)                                       gccgtcAACGGCGTCAGGAGATGGCAGATTTTTCC‑TTCGCGCCGCACGGGCAGCAGGCGGGCAGTCACTTCATTTCCTAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAg                                                                                                                                                                                                            <  1:61481/238‑1 (MQ=255)                                        tcgccaccggccgcagcagATGGCAGACTTTTCC‑TTGGCGGGTCACGGTCAGCGGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                           <  2:737891/226‑1 (MQ=255)                                        ccgtcATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                           <  1:340211/239‑1 (MQ=255)                                         cgtcATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGGATAATCCAGTGGCGGTG‑TTTTCCCCCCCGGTCCCCGGCCATTTACAGTGCCCGTg                                                                                                                                                                                                          >  1:855322/1‑239 (MQ=255)                                         cgccATCGACGTCAGCAGTGTGCAGCCTTTTCC‑TTTGGGCGTTCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                           <  1:582100/235‑1 (MQ=255)                                               tccctccggcaggtgggcgtcTTTTCC‑TTTGCCCGTCCCGTTCAGCCGCCGGGCTGTCCCTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCTGTTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGtaat                                                                                                                                                                                                    <  2:271304/219‑1 (MQ=255)                                               cggcgtcgtcagATGGCATAGTTGTCT‑TTTCGGTCTCACTGTCCTCAGGCGGGCAGTCACGTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGGCTGTGGGCTGGTTTTCAGCTTCACCAGGTGATCGCTTTGGCCAAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCATGCATTTAACAGTCCCAGTGAGtaac                                                                                                                                                                                                    <  2:151120/230‑2 (MQ=255)                                               cggcGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCCGGCGGGCAGTCCCTTAATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGtaat                                                                                                                                                                                                    <  1:889627/236‑1 (MQ=255)                                                        cagtttgcctgtttttCC‑TTTGCCCGTCCCGGGCAGCAGGGGGGCAGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTAt                                                                                                                                                                                          <  2:117236/227‑1 (MQ=255)                                                         agATGGCAGACTTTTCC‑TTTGCGCGTCCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTAt                                                                                                                                                                                          <  1:16262/238‑1 (MQ=255)                                                              ctcgttttttCCTTTTGCG‑GTCCCGGTCCACGGGCGGGCATTCCCTTCTTTTCCCAGCCCCGGCCACTTTTTTCGGGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTa                                                                                                                                                                                   <  2:261719/233‑1 (MQ=255)                                                                    ttttCC‑TGTGCGCGTCGCGGTCAGCAGGGGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTGCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCa                                                                                                                                                                             <  1:243073/239‑1 (MQ=255)                                                                     aagagccgggggtgTCCCGGTCGGCAGGGGGGCAGTCTCTTCTTTTCCCAGTCCCGGCCACTTTTTTCGGGCTTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCa                                                                                                                                                                             <  2:18080/226‑1 (MQ=255)                                                                       tCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCCCCCGGCTCCAGGCACTTTACAGTCCCCGTGAGTAATAAA‑CTTTATCCA‑TTTACGTCataa                                                                                                                                                                          >  1:676470/1‑236 (MQ=255)                                                                                    cACGGTCAGCGGGCGGGCAGTCACTCCA‑TCCCCAGTCCCGGCCACTTTTTCCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGt                                                                                                                                                            <  1:506043/239‑1 (MQ=255)                                                                                     ctgttccggcggggggCAGGCTCTTCTTTTCCCAGCCCCGGCCACTTTTTTCGTGCCTGGGGGCTGGTTTTCAGCTTCACCAGATCATCGTCTTTACCCAGTTTTCTGATCTCTGCATATT‑GCGCTCCCTTTCTGAGAGTTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGt                                                                                                                                                            <  2:432667/226‑1 (MQ=255)                                                                                      ctgtcaccgggggggcgtgcccTTCGTTTCCCTGTGCCGGCCCATTTTTACGTGCCTGCGGGCTTGTATTCCGCTTACCCAGATGACCTCCGTTACCAAGTTTTCTGATCTATTCATATC‑GCGCTCCCTTTCTGAGAGGTACCATCCATGGGCGGTT‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TGAACGTCAGAGT‑GTTATCGCCGGtt                                                                                                                                                           <  2:161933/218‑1 (MQ=255)                                                                                       ggTCCGCAGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCGGGTTTTCAGTTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGttt                                                                                                                                                          <  2:219754/239‑1 (MQ=255)                                                                                         tCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCGCCCAGCCCCACGCCTTTTACAATACCAGTGGGTAATAACACTTTTTCCC‑TT‑ACGTCAAAGG‑GTTTTCGCTGGtcctt                                                                                                                                                        >  1:18080/1‑235 (MQ=255)                                                                                           agcagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                     <  1:914524/179‑1 (MQ=255)                                                                                            gaagggggtaaGTCTCTTTATTTCCCAGTCCCGGCCACTTTTTGCGTGCCTGCGGGCTGGTTTGCAGCTTCACCAGATGATCGCCTTTACCCAGTTTGCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCt                                                                                                                                                     <  2:445415/229‑1 (MQ=255)                                                                                            ccgggggggCATTCTCTTCCTTTCCCAGCCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCt                                                                                                                                                     <  2:676470/233‑1 (MQ=255)                                                                                             cagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCGCCACGCTCCCGGCATTTTAAAGTACCCGTGAGTGAAAAA‑CTTTTTCCC‑ATAACCTCAGAAG‑GTTTTCTCCggggttttttt                                                                                                                                                    >  1:722299/1‑230 (MQ=255)                                                                                                          aCTTCATTTCCCAGGCCCGGCAACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCGGAGCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGc                                                                                                                                       <  1:81351/239‑1 (MQ=255)                                                                                                            ttCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCAttt                                                                                                                                                                                                                       <  1:919557/160‑1 (MQ=255)                                                                                                                      aGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTTCCGCCCAGCTCCAGGCCTTTAACAGTCCCCGTGAGGAATAAA‑CCTTTTTCA‑TTTACGTCAGAAA‑GGTTTTGCCGGGTTGTTTTTTAAAtttgttcgcaaagccatttt                                                                                                                            >  1:310494/1‑220 (MQ=255)                                                                                                                       ggcccGGCCCTTTTTGGCGTGCCTTGGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                          <  1:348699/237‑1 (MQ=255)                                                                                                                             gCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCCCCCCGCGCCCGGCATCTTACCGACCCCCGGAGTAAAAAC‑CTTTTTCCC‑TTAACACCACGGG‑GTTATAGCCCGGTTGTTTTTTAAATTtctccaacagaccattttttctgt                                                                                                                      >  2:169062/1‑217 (MQ=255)                                                                                                                                                cGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTTCCCGCCACGCTCCCGGCCTTTTAAAATCCCCGGGAGTAATAAC‑CCTTATCCC‑TTAACGTCAGAGT‑GTTTTCGCCGGGTTGTTCTATAAATTTCTCAGCAAACTCATTTTCTCTGTTCTTCATCGTGCCGAAGGCt                                                                                                  >  1:904830/1‑238 (MQ=255)                                                                                                                                                   gCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCa                                                                                               <  1:873907/238‑1 (MQ=255)                                                                                                                                                            ccgtttccccaagttatCGCCTTTACCCGGTTTGTTGATCTTTTCAGTTG‑TCGCTCCCGTTTTGGGGGTTATCACCCAGGGGCGGTG‑TTCTCCCGCCAGGCTCCAGGAATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTcag‑ca                                                                                    <  2:650777/224‑1 (MQ=255)                                                                                                                                                               cTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                          <  1:300333/199‑1 (MQ=255)                                                                                                                                                                 gcaccggtgttcccccTTTACCCTGTTTTCGGTTC‑CTTCACATTGGGGCTCCCTTTCCGAGAGGTTTCATCCGGTGGCGGTG‑TTCTCCCGCCAGGCTGCAGGCATTTAACAGTCCCAGTGAGGAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GGTAACGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATgg                                                                               <  2:651870/225‑1 (MQ=255)                                                                                                                                                                    ccAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTg                                                                            <  1:169062/239‑1 (MQ=255)                                                                                                                                                                        atgatCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCaga                                                                                                                                                                           <  1:384994/142‑1 (MQ=255)                                                                                                                                                                                      ccAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑g                                                         <  1:750506/239‑1 (MQ=255)                                                                                                                                                                                       ctgTTTTCTTTTCTTTTCATTTT‑GCCCTCCCTTTCTGGGAGGTACCATCCAGTGGCGGTG‑TTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑Ga                                                        <  2:343748/237‑1 (MQ=255)                                                                                                                                                                                          ttttCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTTACAGTCCCAGTGAGGAATAAA‑CTTTATCCC‑TTAACGTCAGAGT‑GGTATCTCCGGTTTGTTCTTTAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCTTTCCCAAGGCTCCAGCCGTCCC‑CAGATGGCTGGTCCATTCCCTCTTGcgggc                                                        >  1:670782/1‑233 (MQ=255)                                                                                                                                                                                          ttttCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                          <  1:571614/172‑1 (MQ=255)                                                                                                                                                                                          ttttCTGATCCCTTCATTTT‑GCGCTCCTTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGaa                                                                                                      <  2:776876/192‑1 (MQ=255)                                                                                                                                                                                                   ctctTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                          <  1:774867/163‑1 (MQ=255)                                                                                                                                                                                                     ctTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GGCCATTTTGACCTGc                                          <  1:75053/239‑1 (MQ=255)                                                                                                                                                                                                          tatT‑GCGCTCCCTTTCTGGGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTa                                     <  1:328941/239‑1 (MQ=255)                                                                                                                                                                                                                   tCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTTCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCg                                                                                                        <  2:152541/166‑1 (MQ=255)                                                                                                                                                                                                                   tCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGGTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GGTATTGCCGGGTTGTTCTATAAATTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTTCCTCTTGCG‑GACCCCTTTTACCTTCGGGGGGAGCGccc                             >  1:443496/1‑238 (MQ=255)                                                                                                                                                                                                                   tCCCTTGCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCg                             <  1:460279/239‑1 (MQ=255)                                                                                                                                                                                                                     ccTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCCCGCTCCCCGCATTTTACAGTCCCCGTGAGGAAAAAA‑CCTTTTCCACTTAACG‑CAGAGTGGTTTTCCCCGG‑TTGTTCTATAAATTTCTCCGCAAACCCCTTTTTTCTGTTTTTCTTCGTTGCGAAAGCAGCCACCCCCAGACGAATGGC‑GGTTCGTTCCCTCTGGCAGGCCCATTTTgtctccggggggaggcgacgg                            >  2:438700/1‑219 (MQ=255)                                                                                                                                                                                                                             gagggTATCATCCGGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑TTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGc                   <  1:415310/236‑1 (MQ=255)                                                                                                                                                                                                                                       ccccggggCGGTT‑TTCTCCCGCCCGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTgg         <  2:904830/233‑1 (MQ=255)                                                                                                                                                                                                                                             ggCGTTG‑TTCTCCCGCGAGGCTCCAGGCATTTAACAGTCCAAGTGAGTAATAAC‑CTTTATCCA‑GTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTAAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCgggg    <  1:139433/238‑1 (MQ=255)                                                                                                                                                                                                                                             ggCGGTG‑TTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCgggg    <  1:843129/238‑1 (MQ=255)                                                                                                                                                                                                                                               ctgtg‑tTTCCCCGTCAGGCTCCAGGCTTGTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGa  <  1:918345/236‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATT‑GCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTG‑TTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCA‑TTAACGTCAGAGT‑GTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAG‑CAGATGGCTGGTCAGTTCCATCTGGCA‑GACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGA  >  NC_000913/4502514‑4502987

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