Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,415	A→G	K271K (AAA→AAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,415	0	A	G	96.5%	170.7 / ‑3.8	57	K271K (AAA→AAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (36/19);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

ACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAC  >  NC_000913/4620169‑4620660                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        aCTCTACGAACTGTGCGAAATCGCCCGTGACGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAg                                                                                                                                                                                                                                                    <  1:59256/251‑1 (MQ=255)   tctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCCACTGCGGTATCACCAAG‑AAAGTg                                                                                                                                                                                                                                                  >  1:111492/1‑251 (MQ=255)     tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTCCAACCTCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGTCCGTGCTGCCGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGa                                                                                                                                                                                                                                                 >  1:189518/1‑250 (MQ=255)       cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGGTGATGAAAAAAACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTTAGCCAACTGCGGTATCACCAAG‑AAAGTGaaa                                                                                                                                                                                                                                               >  1:409257/1‑250 (MQ=255)       cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCCGCACCGACCGGGCTGCAGAAACTGGGTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCCGACATCTACCCCCACTGCGGGATCACCACG‑AAAGTGaaa                                                                                                                                                                                                                                               >  2:113761/1‑250 (MQ=255)       cGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAg                                                                                                                                                                                                                                              <  2:189518/251‑1 (MQ=255)         accTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCCCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGca                                                                                                                                                                                                                                            <  2:17005/249‑2 (MQ=255)                       gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCCCGTCCGTTTATCGGTGACACAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAG‑AAAGGTAAAGCAACTGGcacggaagc                                                                                                                                                                                                                              >  2:155146/1‑243 (MQ=255)                             gAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTACGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGtt                                                                                                                                                                                                                        >  1:447211/1‑251 (MQ=255)                                gagcggaacaagggaggcTAAAAAATCGGTCGGGTTATCGCGAGGCAGTGTAACGGTGAAAAAGACGGTAACGTCCAGCGTACCGGTAACCGGAACGACCGGGCGGTTGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGGTGAGGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAAATGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGa                                                                                                                                                                                                                     <  2:163428/236‑1 (MQ=255)                                    tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                                                                                 >  1:86622/1‑251 (MQ=255)                                         aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGTTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                                                                                            <  1:46538/251‑1 (MQ=255)                                                                      cGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGCGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAac                                                                                                                                                                                >  1:25325/1‑250 (MQ=255)                                                                             ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                                                                                        >  1:61734/1‑251 (MQ=255)                                                                             ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                                                                                        >  1:33991/1‑251 (MQ=255)                                                                                        tgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGACTTCACCA‑ACt                                                                                                                                                            >  2:98955/2‑251 (MQ=255)                                                                                                 cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACt                                                                                                                                                   >  1:241862/1‑251 (MQ=255)                                                                                                     aaCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCg                                                                                                                                                <  1:260362/250‑1 (MQ=255)                                                                                                     aaCTTCCAGCGTACCGGTAACCGGCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCg                                                                                                                                                <  2:233418/250‑1 (MQ=255)                                                                                                      aCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAACGAGATGAAAGAAGCGGGAGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGAc                                                                                                                                              >  2:192085/1‑251 (MQ=255)                                                                                                          ccAGCGAACCGGTAACCGTCACGACATGGCTGTTGTGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑CAAGTGAAAGCAGCTGGCCGGGACGCGCTGTTTGAGGCCACCATCAAAGAGATGAACGAAGAGGGAGATAACACCATCGGCTTCACCA‑ACTTCGTTGACTTCGActttt                                                                                                                                          >  2:183883/1‑248 (MQ=255)                                                                                                                    ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGTCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCC‑ACTTCGTTGACGTCGACTCTTCCTgggccc                                                                                                                                 >  1:268854/1‑247 (MQ=255)                                                                                                                    ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACAATCGTCTTCACCA‑ACTTCGTTGACTTCTACTCTTCCTGGGGcc                                                                                                                                 >  1:264961/1‑250 (MQ=255)                                                                                                                                        tGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGtt                                                                                                            >  2:134740/1‑251 (MQ=255)                                                                                                                                                    aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                                                                                      <  2:213141/134‑1 (MQ=255)                                                                                                                                                    aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                                                                                      >  1:213141/1‑134 (MQ=255)                                                                                                                                                     gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                                                                               >  1:130684/1‑251 (MQ=255)                                                                                                                                                       accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAATCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCAACGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                                                                             <  1:415435/251‑1 (MQ=255)                                                                                                                                                             cgTGCTGCAGAAACTGGTTGTTGACACACAAGGTCACGTGGTTTATGTCGGTAAAATTGCCGGCATCTACGCAAACTGCGGTATCACAAGG‑AAAGCGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAGGCGGGTGCTAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTCTGCTGCGGGTATGGAGCTGt                                                                                       >  2:301998/1‑251 (MQ=255)                                                                                                                                                                gctgaAGAAAATGGTTGAAGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATGGCAGACATCGACGCCAACTGCGGTAACACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCg                                                                                    <  2:373389/251‑1 (MQ=255)                                                                                                                                                                  tgcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                  >  1:320015/1‑251 (MQ=255)                                                                                                                                                                     aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                                                               >  1:370290/1‑251 (MQ=255)                                                                                                                                                                     aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                                                               >  1:492654/1‑251 (MQ=255)                                                                                                                                                                                   gaAAAACTCGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGGGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCt                                                                 <  1:459015/251‑1 (MQ=255)                                                                                                                                                                                   gaAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGc                                                                  <  2:307761/250‑1 (MQ=255)                                                                                                                                                                                     aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtt                                                                >  2:457614/1‑249 (MQ=255)                                                                                                                                                                                       aacacgttcAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATg                                                             <  1:432244/243‑1 (MQ=255)                                                                                                                                                                                        acacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtgttgt                                                            >  2:157239/1‑251 (MQ=255)                                                                                                                                                                                         cacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGActc                                                                                                                                            >  2:14200/1‑170 (MQ=255)                                                                                                                                                                                         cacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGActc                                                                                                                                            <  1:14200/170‑1 (MQ=255)                                                                                                                                                                                                 ggtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCtgctgc                                                   >  1:351695/1‑251 (MQ=255)                                                                                                                                                                                                 ggtggatcatgacGGTAAAATTGCAGACATCTCCGACAACTGCGGTATAACCAAG‑ACAGACAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTC‑CCACACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTTTGGGGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCtgctgc                                                   >  2:493516/13‑251 (MQ=255)                                                                                                                                                                                                          tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCTTCTGCCGTGGCTGATGTCTCTTCTGCGCGagtacg                                          >  2:348502/1‑246 (MQ=255)                                                                                                                                                                                                             cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                     <  2:361863/237‑1 (MQ=255)                                                                                                                                                                                                             cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                     >  1:361863/1‑237 (MQ=255)                                                                                                                                                                                                              ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                      >  2:406045/1‑251 (MQ=255)                                                                                                                                                                                                              ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                      >  1:319837/1‑251 (MQ=255)                                                                                                                                                                                                                  aaaTTGCAGACATCTACGCCAACTGCGGTATCACCACG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTAGGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTg                                  >  2:73568/1‑251 (MQ=255)                                                                                                                                                                                                                       gcagACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCt                             <  1:87975/248‑1 (MQ=255)                                                                                                                                                                                                                          gACAGCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACCCCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCAc                          <  2:319837/251‑1 (MQ=255)                                                                                                                                                                                                                               cTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCActgctg                     >  1:373871/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGAc                   >  1:145314/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCCCTGCTGAcccc                >  2:390609/1‑249 (MQ=255)                                                                                                                                                                                                                                     cAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACCTCCTGCTCCTCACTGCTGACCACg               >  2:401039/1‑251 (MQ=255)                                                                                                                                                                                                                                        acgtcgttataCCAAAGAAAGTTAAAGCAACGTGCCTGGACGCGCTTTTTGACGCCACCATCAAAGAGAGGAAAGAGGCGGGTGATATCACCATCGTCTTCTCCA‑ACTGCTTTGACTTCGACGCTTCCTGGGGCCACCGGCGCGCCGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGt             <  2:268854/241‑1 (MQ=255)                                                                                                                                                                                                                                                 cACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCg    <  2:373871/250‑1 (MQ=255)                                                                                                                                                                                                                                                  aCCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGAc  <  1:73568/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        ACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAC  >  NC_000913/4620169‑4620660

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