Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,223	T→C	N207N (AAT→AAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,223	0	T	C	96.3%	146.9 / ‑6.4	55	N207N (AAT→AAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (28/24);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCAT  >  NC_000913/4619979‑4620459                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                aaggccgggtggaaaaccTGGTCGAACGCCCTAACCGCCCGGGTTACCCGGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                                                                                                                                                                                           <  2:370539/234‑1 (MQ=255) aaGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGACGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATACACTCTACGAACTGTGCGCAATCGCCCGTGAAGCGCTGACCACCGGCGGCTAAAACATCggg                                                                                                                                                                                                                                           >  2:87975/1‑250 (MQ=255)  agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTc                                                                                                                                                                                                                                          >  1:388624/1‑251 (MQ=255)           ggATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCg                                                                                                                                                                                                                                 >  1:17005/1‑251 (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)                                                                                                                                                                      acgaaggatccttgggccTGGACAAACTCTACGAATTGTGCGAAATCGCCCGTGACGAGCGGACCACAGGCGGCTTTAATATCGGTGGTGTTATCGCTCGTCCCTTTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACTGCCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTAc                                                                   <  2:377518/234‑1 (MQ=255)                                                                                                                                                                                 tttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGCAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTAt                                                                                                                                                                                                                                   >  1:450729/7‑83 (MQ=255)                                                                                                                                                                                 tttggccTGGATAAACTCTACGAACTGTGCGAAAGCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACAGCGGTCGTGTTAt                                                                                                                                                                                                                                   <  2:450729/77‑1 (MQ=255)                                                                                                                                                                                           aTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGACGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTTACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGCGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAAG‑GGGTTTATGTCGGTAAAATTGCAGGCATCTAAGCCCACAGCGGTAttcaca                                                >  2:305268/1‑244 (MQ=255)                                                                                                                                                                                            tAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCa                                                >  1:420452/1‑248 (MQ=255)                                                                                                                                                                                              aaCTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCAGACCTCTACGCCAACTGCGGTATCAcaaagaa                                            >  1:163271/1‑244 (MQ=255)                                                                                                                                                                                               aCTCTACGAACTGTGCGAAATCGCCCGTGACGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAg                                          <  1:59256/251‑1 (MQ=255)                                                                                                                                                                                                 tctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCCACTGCGGTATCACCAAGAAAGTg                                        >  1:111492/1‑251 (MQ=255)                                                                                                                                                                                                   tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTCCAACCTCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGTCCGTGCTGCCGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGa                                       >  1:189518/1‑250 (MQ=255)                                                                                                                                                                                                     cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGGTGATGAAAAAAAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTTAGCCAACTGCGGTATCACCAAGAAAGTGaaa                                     >  1:409257/1‑250 (MQ=255)                                                                                                                                                                                                     cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCCGCACCGACCGGGCTGCAGAAACTGGGTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCCGACATCTACCCCCACTGCGGGATCACCACGAAAGTGaaa                                     >  2:113761/1‑250 (MQ=255)                                                                                                                                                                                                     cGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAg                                    <  2:189518/251‑1 (MQ=255)                                                                                                                                                                                                       accTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCCCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGca                                  <  2:17005/249‑2 (MQ=255)                                                                                                                                                                                                        aCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACCACATCGGTCGTTTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAGAAAAAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCAGACATCATCGCCAACTGCGGTATCACCacgaaaaggaaagca                                  >  2:258683/1‑236 (MQ=255)                                                                                                                                                                                                                     gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCCCGTCCGTTTATCGGTGACACAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAGAAAGGTAAAGCAACTGGcacggaagc                    >  2:155146/1‑243 (MQ=255)                                                                                                                                                                                                                           gAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTACGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGtt              >  1:447211/1‑251 (MQ=255)                                                                                                                                                                                                                              gagcggaacaagggaggcTAAAAAATCGGTCGGGTTATCGCGAGGCAGTGTAACGGTGAAAAAGACGGTAACGTCCAGCGTACCGGTAACCGGAACGACCGGGCGGT‑TGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGGTGAGGAAAAACAC‑GGTCAGG‑GGGTTTCTGTCGGTAAAAATGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGa           <  2:163428/236‑1 (MQ=255)                                                                                                                                                                                                                                  tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc       >  1:86622/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGTTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt  <  1:46538/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                AAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGT‑TGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCAT  >  NC_000913/4619979‑4620459

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