Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,150	A→G	E183G (GAA→GGA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,150	0	A	G	96.2%	146.8 / ‑3.0	53	E183G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (27/24);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

ACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGC  >  NC_000913/4619907‑4620384                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              cctggaaaaccgCCGGTGTCCCGGTTCTGTTTGAGTGGTGACATTTCTCCAATCACGAAAACAGCTTCCTGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAACCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAccttt                                                                                                                                                                                                                                        <  1:92634/240‑5 (MQ=255)     ggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATTAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTc                                                                                                                                                                                                                                    >  1:100512/1‑251 (MQ=255)      aaaaTCGCCGGGGTCCCGGTTTTGTTTGAGTGGGGTTATTTCCCCGATCCCGAAAACAGCTTCCCGCAGGACCTGCTGGATAACCTGGT‑CGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcct                                                                                                                                                                                                                                   <  2:320594/250‑10 (MQ=255)       aaaTCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcctg                                                                                                                                                                                                                                  >  1:174773/1‑241 (MQ=255)                   ccccGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTCCCCGCAAGAGCTGCTGGATAAACTGTC‑CGAACGCGCAAATCTGCCGGTTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTAc                                                                                                                                                                                                                      <  1:267797/250‑1 (MQ=255)                     ccGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTACCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACg                                                                                                                                                                                                                     >  1:450035/1‑250 (MQ=255)                        gTTCTGTTTGAGAGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACt                                                                                                                                                                                                                 >  1:340907/1‑251 (MQ=255)                            tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATCAACTATACGAACTGTGc                                                                                                                                                                                                             >  1:476691/1‑251 (MQ=255)                                       ggtttttcctcagtacACGAAAACGGCTTCCCCCAAGGTCGGCTGGTGAAACGGGT‑CTACCGGCCTAACCTGCCGGTTTCCCTGGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGTGAAGAGCACATGAAAACCGGCAATCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCg                                                                                                                                                                                                  <  2:381693/236‑1 (MQ=255)                                         atatTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCAGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATTAAAACCGGCAAGCCGATTTTCTATACCTCAGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGAAACACTCAACGAACTGTGCGAAATCGCCCggg                                                                                                                                                                                                >  2:314170/1‑249 (MQ=255)                                             ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGa                                                                                                                                                                                               >  1:231622/1‑248 (MQ=255)                                                       aCGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGAGATACTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                                                                                                  >  1:108897/1‑251 (MQ=255)                                                         gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                                                                                                >  1:126720/1‑251 (MQ=255)                                                         gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                                                                                                >  1:474348/1‑251 (MQ=255)                                                         aaaaaCATCTTCCCGCAAGTGCTGCTGGAAAAAC‑GTTACGAACCCTCTATTCGGCCGGGTTACCCCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcg                                                                                                                                                                                 <  2:474348/249‑1 (MQ=255)                                                                         aaggccgggtggaaaaccTGGT‑CGAACGCCCTAACCGCCCGGGTTACCCGGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                                                                                                                <  2:370539/234‑1 (MQ=255)                                                                         aaGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGACGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATACACTCTACGAACTGTGCGCAATCGCCCGTGAAGCGCTGACCACCGGCGGCTAAAACATCggg                                                                                                                                                                >  2:87975/1‑250 (MQ=255)                                                                          agagCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTc                                                                                                                                                               >  1:388624/1‑251 (MQ=255)                                                                           gagCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAGGGAACCTTTGGCCTGGATAAACTCTACGACCTGTGCGAAATCGCCCGTGAAGAGCTGCCCAACggaggcacacacacccgtcg                                                                                                                                                              >  2:260362/1‑237 (MQ=255)                                                                                   ggATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCg                                                                                                                                                      >  1:17005/1‑251 (MQ=255)                                                                                      tAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGCCTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGATAccacctccgtctggatatcaccc                                                                                                                                                   >  2:459015/1‑229 (MQ=255)                                                                                                       ccTAACCTGCGGGGTTACCTCGGTAACTGCCACCCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTACTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCTAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTCATCGCTCGTCCGCTTATCggtga                                                                                                                                   <  2:454817/250‑4 (MQ=255)                                                                                                             cggCCGGTTTACCCCGGCACCTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGc                                                                                                                             <  2:388624/249‑1 (MQ=255)                                                                                                                       acctcgtaacctgccaacctTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGTCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCTTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACtt                                                                                                                    <  2:225887/232‑1 (MQ=255)                                                                                                                            ggTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCg                                                                                                              <  1:99616/251‑1 (MQ=255)                                                                                                                             gTAACTGCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGt                                                                                                             <  1:233809/251‑1 (MQ=255)                                                                                                                                    ccACACTTCCGGTACTGTAATTATGGATCAGCTGGGCAAAGAACACATGAAAACCGCCAAGCCGAGTTTCTATACCTCCGCTGACTCCGTGTCCCAGATTGCCTGCCACGAAGGCACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCATCGGCGGCTACAACATCTGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCgtta                                                                                                      <  2:116511/251‑4 (MQ=255)                                                                                                                                          ttCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGCTTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGACGGTAACTTCCTGCGTACCGGTaaaccac                                                                                                >  2:78625/1‑246 (MQ=255)                                                                                                                                            ccGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAc                                                                                              <  1:258111/251‑1 (MQ=255)                                                                                                                                            ccGGTACGGGCATTCTGGAGTAGTGGGGCTAAGAGCCCTTGAAAACCGTCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCCCGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGTCGGTA‑CTTCCAGCGTACCGGTAACCGTCACg                                                                                             <  2:373706/251‑1 (MQ=255)                                                                                                                                             cGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCATCGTACCGGTAACCGTCAc                                                                                              >  1:373389/1‑250 (MQ=255)                                                                                                                                                   ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCggtga                                                                                                                                   >  1:124045/1‑204 (MQ=255)                                                                                                                                                            gtaccaggggggCGAAGAGCCCATGACAACCGTCAAGCCGTTTTTCTAAACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAg                                                                              <  2:340907/241‑1 (MQ=255)                                                                                                                                                                  gcggggCGAAGAGCACAGGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGAGTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAg                                                                       <  2:7576/248‑1 (MQ=255)                                                                                                                                                                      ggCGAAGAGCACATGAAACCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCacc                                                                   <  1:448435/251‑1 (MQ=255)                                                                                                                                                                        cGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCAACAACATCGGGCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAAGACCTTGCTGT‑TGAgacgcccgaccccg                                                                 >  2:441806/1‑243 (MQ=255)                                                                                                                                                                             ggccaaTGAAAACCGGCAAGCCGATTTTCTGTACCCCCGCTGACTCCGTGTTCCAGAGTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAGCCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGtg                                                            <  2:231622/246‑1 (MQ=255)                                                                                                                                                                             agCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGATGACTCCGTGTTCCAGATTGCCTGCCACGACGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGCCGGACACAACCTCGGACGTGTTATCGCTCGTCCCATTATCGGTGACAAAGCCGGTAAAATCAAGCGTACCGGTAACCGTCACCACCTGGccatttaggaccccacaaccaccgtg                                                             >  2:67650/1‑226 (MQ=255)                                                                                                                                                                                acaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACATCATCGGTCGTGTTATCGCTCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGtgctg                                                         <  1:246634/251‑1 (MQ=255)                                                                                                                                                                                       aaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTACGATTATCGGTGACAAAGCCGGTAACTTCCACCGTACCGGGAACCGTCGCGACCTTGCAGT‑TGAGCCGCCAGCCCAGACCGCGCTGCCGaaa                                                   >  2:484740/1‑250 (MQ=255)                                                                                                                                                                                          cGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGGTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                               <  2:126720/251‑1 (MQ=255)                                                                                                                                                                                           ggCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCAGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGAGACAAAGCCGGGAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCCCCGCCCGTGCTGCAGAAACTggg                                              >  1:67907/1‑250 (MQ=255)                                                                                                                                                                                               aGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGCCCAACGGCCGCCACAACATCGGTCGTGATATCGCTCGTCCGTTTATCCGTGACAAAGCCGGTAACTTCCAGCGTACCGGAAACCGTCACGACCTGGATGT‑TGAGCCGCCAGCACCGACCGAGCTGCAGAAACTGGTtgat                                          >  1:73554/1‑251 (MQ=255)                                                                                                                                                                                                      tttcattaCCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGTCCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGaaaaa                                    <  2:397852/244‑1 (MQ=255)                                                                                                                                                                                                      tttCTATCCCTCCGCTGGCTCCGCGTTCCAGCTTGCCTGCAAAGAAGGAACCTTTGGCCTGGCTACACTCTCCGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCACCGGCGGCTACAACCTCGGACGTGGTTTCGCTAGTCCGTTTCTCGGTTACCAAGCCGGGCACTTCCCGCGGAGCTGGCACCGTCCCGACCAGGC‑GTGTGTGCCGCCAGCCCCGCCCGGGCTGCAgctacttgtgggagaacgac                                   >  2:220221/1‑232 (MQ=255)                                                                                                                                                                                                         ccttccctcccatggcTCCGTGATCCAGATGGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACAACCTGTGCGATATCGTCCGTGAAGAGCGGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGCAACTTCCAGCGTACCGGGAACCGACACGACCTGGATGT‑TGAGCCGCCAGCCCCGACCGTGCTTCAGAAACTTGTTGATGAAAAACAC‑g                               <  2:31290/236‑1 (MQ=255)                                                                                                                                                                                                                 ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAgg‑t                       <  1:78625/250‑1 (MQ=255)                                                                                                                                                                                                                   gCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAgg‑tggt                    >  1:307761/1‑251 (MQ=255)                                                                                                                                                                                                                      tatcccGTGTTCCAGTTTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGACCTGTGCGAAACCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGttt                  <  2:503852/246‑1 (MQ=255)                                                                                                                                                                                                                          ccGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAAAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTGT‑TGAGCCGCCAGCACCGGCCGTGCTGCAGAAACTGGTTGATGGAAA‑CACCGGGCAGG‑GGGTTTCTGTc             >  1:52598/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑GGGTTTCTGTCGGTaaa       >  1:233418/1‑250 (MQ=255)                                                                                                                                                                                                                                  cAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCCCGACCTGGCTGT‑TGGACCCCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGGAAAACAC‑GGTCAGGGTGGTT‑CTTTCGGGAAAAt     >  2:393229/1‑251 (MQ=255)                                                                                                                                                                                                                                      ttgccgtccATGAAGAAACTTTCGGTCTGGATAACCTCTCCGAACTGTGCGAAATCGCCCTTGAAGAGCTGACCAACGGCGGCGACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGc  <  2:100512/243‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGT‑CGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGC  >  NC_000913/4619907‑4620384

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