Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,213	A→C	S190A (TCG→GCG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,213	0	A	C	95.3%	192.8 / ‑6.2	65	S190A (TCG→GCG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (28/33);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCAC  >  NC_000913/64970‑65459                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       cTTTAGCCTGGGCAAAACAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGAGCGGAGGCg                                                                                                                                                                                                                                                 >  1:154294/1‑251 (MQ=255)  tttAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGt                                                                                                                                                                                                                                                >  1:123459/1‑251 (MQ=255)  tttAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGt                                                                                                                                                                                                                                                >  1:634338/1‑251 (MQ=255)               aaaaaGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATtctctgg                                                                                                                                                                                                                                   <  1:170983/251‑4 (MQ=255)                  aaGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGCTCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGCAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGAGGCTAGCGACGTAATCCAGTTCAAAATCGATCGCGAAGGCGTCGCCATTATCTGGccc                                                                                                                                                                                                                                >  2:485439/1‑251 (MQ=255)                  aaGACGCCGTTTTTAAAGCCGTGCTCGAGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCGCACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGGTGGCGTCGTATTCCAGATCAAAATCGCGCGGTGAGGCGACGCCATtcctgggccc                                                                                                                                                                                                                                >  2:80995/1‑243 (MQ=255)                   aGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGGATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGcccc                                                                                                                                                                                                                               >  1:465526/1‑251 (MQ=255)                     ccGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAg                                                                                                                                                                                                                             <  1:507344/250‑1 (MQ=255)                                aaagcctgccccgcgCCACCCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCCCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGa                                                                                                                                                                                                                  <  2:576702/240‑1 (MQ=255)                                 aaGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGAt                                                                                                                                                                                                                 >  1:255463/1‑251 (MQ=255)                                   gCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACAAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGATGGCGACGTATTCCAGTTCAAAATCGCGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATgc                                                                                                                                                                                                               >  1:150555/1‑251 (MQ=255)                                        gCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGGTCAAAATCGAGCGCGGCGGCGTCGCAATACTCTGGCCCCCGCATAAAAACGATGCGCTGc                                                                                                                                                                                                          >  2:531674/1‑251 (MQ=255)                                              gcCACTCCAGCTCACTGTTATCGCGCCCCACCCGCAGCGGACTACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTTGCGCAGTATTCCATTTCAAAATCGGGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAACCGGTGCGCTccccccac                                                                                                                                                                                                    >  2:438738/1‑247 (MQ=255)                                               ccACTCCAGCTCACTGGTAACGCGCCCCAGCCGCAGCGGAATACGGTAACGCGCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGc                                                                                                                                                                                                   <  2:45765/251‑1 (MQ=255)                                                 aCTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCt                                                                                                                                                                                                 >  1:163984/1‑251 (MQ=255)                                                 aCTCCAGCTCACTGTTATCGCGCCCCACCCGCAGCGGCATACGGTAACGATCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGCTGATCCCATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCATAATCGAGCCACGAGGCGTCGCCATTATCATGCCCCAGCATATaacacatgccctcaccgcatact                                                                                                                                                                                                 >  2:303213/1‑230 (MQ=255)                                                    ccAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTc                                                                                                                                                                                               <  1:91577/250‑1 (MQ=255)                                                     cAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTGCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCa                                                                                                                                                                                             <  1:352892/251‑1 (MQ=255)                                                                 atCGCGCCCCAGCCGCAGCGGCATAAGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGCTCCAACCAATGATCACCTCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCAAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGCGCGAGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGt                                                                                                                                                                                 >  1:350361/1‑251 (MQ=255)                                                                 atCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGt                                                                                                                                                                                 >  1:705545/1‑251 (MQ=255)                                                                   cgcgCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTAc                                                                                                                                                                               <  2:13384/251‑1 (MQ=255)                                                                      gccccagccGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGc                                                                                                                                                                            >  1:36639/9‑251 (MQ=255)                                                                          cagccGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCAc                                                                                                                                                                        >  1:157707/5‑251 (MQ=255)                                                                          cagccGCAGCAGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCCCGTTCCAACCGATGCTCACATAGGGATCGTCGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAACATCGAGCGCGGAGGCGACGCCATTATCTAGCCCCAGCATATAACCGCTGAGCAGCCCGCAGCCTTCCAGGCCGCTGCAACACAGCTCAc                                                                                                                                                                        >  2:640415/5‑251 (MQ=255)                                                                               gCAGCGGAATTCGGTAACGCTAGGCATGTTTTTGCAGCATTCGCAGATCGAACTGAACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAGACCCGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGAGGTGGGGTCGACATTCTCTGGCCCCAGCCTATAAAAGATGCGCTGCCCGCAGCCTTCCAGGCCGAGGCAGTACAGCTCACCGTgg                                                                                                                                                                   >  2:81236/1‑251 (MQ=255)                                                                                 aGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCg                                                                                                                                                                 >  1:99458/1‑251 (MQ=255)                                                                                     gaataCGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGggtg                                                                                                                                                             >  1:522529/5‑251 (MQ=255)                                                                                          cGGTAACGCTCGGCATGTTTTTGCAGCAGTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAGACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTc                                                                                                                                                        <  1:81236/251‑1 (MQ=255)                                                                                           ggTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTc                                                                                                                                                        >  1:470575/1‑250 (MQ=255)                                                                                                 gCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGCAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGCGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATc                                                                                                                                                  >  1:637742/1‑250 (MQ=255)                                                                                                     ggCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAg                                                                                                                                             <  1:78215/251‑1 (MQ=255)                                                                                                          gtgtttGCAGCAATCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAcc                                                                                                                                        <  1:52521/248‑1 (MQ=255)                                                                                                          gTTTGTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAc                                                                                                                                         <  1:710118/250‑1 (MQ=255)                                                                                                             tttGCATCATTCGCAGATCGAACTGCCACGCGGTCCAACCGATGATCACAGCAGGAGCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCATTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCa                                                                                                                                      <  1:664030/250‑1 (MQ=255)                                                                                                                gcagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACtt                                                                                                                                   <  1:2371/250‑1 (MQ=255)                                                                                                                 cagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCCCACCGTGGCGGGGGGGTTCAATATCCACAGAAAACCACTTg                                                                                                                                  >  2:533494/1‑250 (MQ=255)                                                                                                                  agcttatgCAGATCTAACTGCCCCACGGTACAACGGTGGACCACATCAAGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGGGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGAGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAg                                                                                                                                <  2:228785/244‑1 (MQ=255)                                                                                                                    caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCCTATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTAAATATCAATACAAACACACTTGTgcg                                                                                                                              >  1:695179/1‑251 (MQ=255)                                                                                                                       tCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAgcggc                                                                                                                            <  1:486210/250‑1 (MQ=255)                                                                                                                         gCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAAAATCAAAAAAAACCCACTTGAGcggcggg                                                                                                                         >  1:579188/1‑250 (MQ=255)                                                                                                                            gATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACa                                                                                                                                                                              <  2:474230/195‑1 (MQ=255)                                                                                                                            gATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACa                                                                                                                                                                              >  1:474230/1‑195 (MQ=255)                                                                                                                                   cTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGATAGTCgg                                                                                                               >  1:488484/1‑251 (MQ=255)                                                                                                                                     gccccacGTTCCAACCGAGGAGCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTGCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCCCTTGAGCGGCGGACGATAGTCgg                                                                                                               <  2:137699/246‑1 (MQ=255)                                                                                                                                        ccacGGTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTGCCAGGAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCGCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGc                                                                                                          <  2:309764/251‑1 (MQ=255)                                                                                                                                               ccAACAGATGATCACATCAGGATCGTAGTTGGCCAACCGGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGTCGACGTATTCCAGATCAAACGCGAGCGCGGAGGCGTCGCAATTCTCTGGCCCAAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGAGAGTCGGGATGCGGTTTCa                                                                                                   <  2:113001/251‑1 (MQ=255)                                                                                                                                                         atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCAtt                                                                                         <  1:263346/251‑1 (MQ=255)                                                                                                                                                         atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGGTCAAAATCGAGCGCGGAGGCGGCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCAtt                                                                                         <  2:710125/251‑1 (MQ=255)                                                                                                                                                             caTCAGGATCTGAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACg                                                                                     <  2:545036/251‑1 (MQ=255)                                                                                                                                                                cAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATa                                                                                  <  1:648453/251‑1 (MQ=255)                                                                                                                                                                cAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATa                                                                                  <  1:167830/251‑1 (MQ=255)                                                                                                                                                                   gATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAgcg                                                                               <  1:47021/251‑3 (MQ=255)                                                                                                                                                                     cacccgacgctcctccgatctttGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAAGCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGcc                                                                             <  2:225542/230‑1 (MQ=255)                                                                                                                                                                        tAGTTGGCAAACCAGGCGTGTATTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCa                                                                            <  2:96670/249‑1 (MQ=255)                                                                                                                                                                               cacaCCAGCCGTTTAGTTTTTCAAGTAGCTGCGGGCGGCTGGACCCGTAGTCCAGTATAAAAGCGAGCGCGTAGCCGTCGCCATTCGCTGGCCCCAGCATATAAACGATGCGCTGCTCGCAGCCTTCCAGGCCGATGCAGTACTGCTCACCGTGGCGGGTTGTTTCAAGATCAATAGAAACCGGCTTTAGTGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCata                                                                   <  2:269291/248‑1 (MQ=255)                                                                                                                                                                                          ttGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCCTTAAAGATAGCGcccgatcgggagagcacacgtcc                                                        >  1:225542/1‑230 (MQ=255)                                                                                                                                                                                                 tttCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCgg                                                                                                                                                                >  1:105682/1‑140 (MQ=255)                                                                                                                                                                                                 tgtccagtagcTGCGGGCGGCTGGCGCCGTATTCCAGGTCAAAACCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCgg                                                                                                                                                                <  2:105682/130‑1 (MQ=255)                                                                                                                                                                                                     cagtagcttcgggcggCTGGCGACGTAGACCAGTTCAAAAGCGAGCGCGGAGGCGTCTCCATTCTCTGGCCCCAGCATATAAACGATGCTCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTTGACCCACACCggtg                                             <  2:470575/242‑1 (MQ=255)                                                                                                                                                                                                              cgggcggCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGGTAGTCGGGATGCGGTTTCAGACGGGCCTTAACGATAGCGCCATTGTGCATATCACCCTCGACCCCCCACGGTGAGGGGATa                                     >  1:575412/1‑250 (MQ=255)                                                                                                                                                                                                                    gcTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCt                              <  1:368271/251‑1 (MQ=255)                                                                                                                                                                                                                       ggCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTcc                            <  1:431363/250‑1 (MQ=255)                                                                                                                                                                                                                               attccatgtcaaaaTCCAGCGCGGAGGCGCCGCCATTCTCTGGCCCCCGCATATAAACGAGGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGGGGCGGGTGGTTTCAATATCTATAGAAACCCACTGGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACTGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATa                    <  2:376303/239‑1 (MQ=255)                                                                                                                                                                                                                                              cgaGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGGTTTAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGGGAGGTGGTAAACCGCTCCATCAGATAGCGGTCTGGCGGgcgc    >  1:119869/3‑251 (MQ=255)                                                                                                                                                                                                                                                aGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTTAATATCAAAAAAAACCCACTTGAGCGGCGGACGGGAGTCGGGGAGCGGGTTCTGAAGGGCGTTAACGATAACGCCCTTTTGGATATCCCCCCCGCCCCCCACCCGGGGGGGGAAAAAACGCTCCCTCAGATAACGTTCTGGCGGagcgcc  >  2:675583/1‑246 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       CTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCAC  >  NC_000913/64970‑65459

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: ‑