Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,028	G→A	S251S (AGC→AGT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,028	0	G	A	96.9%	198.4 / NA	65	S251S (AGC→AGT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/23);  minor base C (1/1);  total (40/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 = 1.00e+00

GCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCC‑TGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCA  >  NC_000913/64799‑65272                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gcgcATGTTTATCTTCGGCGTAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTTTCGATAGATTTTCCTTTGCCTCATTACTCC‑TGAGCGCCGGTTTCCCCCGAGAATTACAAGAAATTCCAGCACGCGGGTTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTATCCTGGGCCAAAAAGACGCCGTTTTTTACGCCGTTCTCTCGCCCCTCCCCCTCCCTGTTTTCGCGccccagcccc                                                                                                                                                                                                                                 >  2:52521/1‑244 (MQ=255)        tttATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGa                                                                                                                                                                                                                           >  1:228785/1‑250 (MQ=255)             cTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGactacgg                                                                                                                                                                                                                     >  1:498561/1‑245 (MQ=255)               tCGGAGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACCCC‑TGAGCGTCGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGGCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACg                                                                                                                                                                                                                      <  2:80321/248‑1 (MQ=255)                       aCGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCCAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCTCCCCAGCCGCAACGGAATACGGTAACGCTCgg                                                                                                                                                                                                           >  1:429507/1‑251 (MQ=255)                            ggTCAATTTCGTCCATTCGCTCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCCACGGAATACGGTAACGCTCGGCATGt                                                                                                                                                                                                      >  2:648453/1‑251 (MQ=255)                               cAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGtttt                                                                                                                                                                                                   <  1:190035/251‑1 (MQ=255)                               cAATTTCGCCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGtttt                                                                                                                                                                                                   <  2:114403/251‑1 (MQ=255)                                    cacccgccgctcttccgatctttATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAAGGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgc                                                                                                                                                                                                <  2:207087/228‑1 (MQ=255)                                          agtctaccccACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCTACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCTGCATGTTTTTGCAGCATTCg                                                                                                                                                                                        <  2:171947/244‑1 (MQ=255)                                            ccgaacccACGGGTTATCGATAGAGTTTCCTTAGCTAAATAACCCC‑TGAGCGACGGTTTCCAGCGAGAAGGAAGAGAAAGTCCAGAACGCGGATATCAGCGCCTCGATACCGCCGATAATTAGCCGTCTGTTAGCCTGGGCAAAAAAGACGCCGCTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCa                                                                                                                                                                                      <  2:363444/246‑1 (MQ=255)                                             cGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGc                                                                                                                                                                                       >  1:515426/1‑249 (MQ=255)                                              gATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGc                                                                                                                                                                                       >  2:91577/1‑248 (MQ=255)                                                     cGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCA‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGCTAATTAGCCGACCTTTAGCCTGGGCACAAAAGACGCCGTTTTTAAAGCCGGTCTCGCGACACTCAAGCTAACTGTTATCGCGCCCCAGCCCCAGAGGAATCCGGGAACGCTCGACATGTTTTTGCAGAATTCGACGATCGAACt                                                                                                                                                                             >  2:167115/1‑251 (MQ=255)                                                         ttATCGATAGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgc                                                                                                                                                                                                >  1:207087/1‑228 (MQ=255)                                                              gaacgATTTTCCTTCGCCTAATA‑CTCCCTGCGCGTCGGTTTCCAGCGGGAATGAAGAGAAATTCCAGAACCCGGATTTCAGCGCCCCGATACCGTTGATAATTAGCTGTCCTTTAGCCTGGGAAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACGCCAGCTCACTGTTATTGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCCCCACGt                                                                                                                                                                    <  2:563692/247‑1 (MQ=255)                                                                 aGATTTTCCTTCGCCTAATAACTCC‑TGAGCGACGGTTCCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTcc                                                                                                                                                                 <  2:515426/251‑1 (MQ=255)                                                                     tttCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAAcc                                                                                                                                                             >  1:309764/1‑251 (MQ=255)                                                                       tCCTTCGCCTAATAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCg                                                                                                                                                            >  1:45765/1‑250 (MQ=255)                                                                                 aaTAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATcaca                                                                                                                                                   <  1:630101/249‑1 (MQ=255)                                                                                 aaTAACTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATcaca                                                                                                                                                   <  1:126111/249‑1 (MQ=255)                                                                                  ataacTCC‑TGAGCGACGGTTTCCAGCGAGAATGAAGAGTAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTTACGCTCGGCATGTTTTTGCAGCATTCGCCGCTCGAACTGCACCACGTTCCAACCGATGATCCCATc                                                                                                                                                 >  2:258972/5‑250 (MQ=255)                                                                                             aGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGt                                                                                                                                      >  1:376303/1‑251 (MQ=255)                                                                                              gCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATCCGGTAACGCTCGGCATGTTTTTGCCGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGTTCGTAGtt                                                                                                                                     >  2:167830/1‑251 (MQ=255)                                                                                              gCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                                                                                                     <  2:146901/251‑1 (MQ=255)                                                                                               ggacggTTTCCAGCGAGAATAAAGAGAAATTCCAGAACGCGGATTTCATCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTAGCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTg                                                                                                                                    <  2:434528/246‑1 (MQ=255)                                                                                               cGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTg                                                                                                                                    >  1:96670/1‑251 (MQ=255)                                                                                                gacggTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                                                                                                   >  1:576702/5‑251 (MQ=255)                                                                                                  cacgacgctcgctcgatacagagAAAGTCCAGACCCCGGATTTCAGCGCTTCGAGCCCTTATATAATTAGCCGACCTTTTGCCTGGGCAAAAAAGACGCCGTTTTGAAAGCCGGGCACGCGCCACTCCAGCTCTCTGTTAGCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGAGCACATCAGGATCGTAGTTGGCa                                                                                                                                 <  2:368130/232‑1 (MQ=255)                                                                                                   ggTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgcag                                                                                                                                                                                              >  1:17692/2‑189 (MQ=255)                                                                                                   ggTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgcag                                                                                                                                                                                              <  2:17692/188‑1 (MQ=255)                                                                                                        ccAGCGAGAATGAAGAGAAATTCCAGACCGCGGATTTCAGCGCCTCTGTACCGTCTATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGAGCACAGCAGGATCGTAGTTGGCAAACCa                                                                                                                            <  2:364266/250‑1 (MQ=255)                                                                                                        ccAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAg                                                                                                                           <  1:466846/251‑1 (MQ=255)                                                                                                                 aTGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGt                                                                                                                  >  1:545036/1‑251 (MQ=255)                                                                                                                    aaGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAATCCAGGCGTTGAGtttt                                                                                                               <  1:574104/251‑1 (MQ=255)                                                                                                                     agagAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATCCGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCa                                                                                                                                 >  1:368130/1‑232 (MQ=255)                                                                                                                            ttCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAgt                                                                                                         >  1:13384/1‑248 (MQ=255)                                                                                                                               cAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAgtagctg                                                                                                    >  1:379666/1‑245 (MQ=255)                                                                                                                                  aaCGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCg                                                                                                                                                                                                            >  1:425848/1‑144 (MQ=255)                                                                                                                                   aCGCGGATTTCAGCGCCTCGATACCGGCGATTATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCACAGCCGCAGCGGTATACGGTCACGCTCGGCATGTTTTTGCAGCCTTCGCAGATCGAACTGCACCACGTTCCAACCGCTGATCACATCAGGATCGTAGTTGGAACACCAGGCGTTGAGTTTTTCCGGTAGCTgcggg                                                                                                >  2:60542/1‑251 (MQ=255)                                                                                                                                     ccgGATTTCAGCGCCTCGTTACCGTCATTAATTAGCCGCCCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCGCGCGCCACTCCAGCTCCCTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTgcgggcg                                                                                              <  2:403346/250‑1 (MQ=255)                                                                                                                                       ggATTTCAGCGCCTCGATACCGTCGCTAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGGGTTATTCCAGTAGCTGCGggcggc                                                                                            >  2:278024/1‑251 (MQ=255)                                                                                                                                                          ccGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTc                                                                         >  1:32553/1‑251 (MQ=255)                                                                                                                                                                 tAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCCCTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGAATTCGCAGATCGAACTGCACCACGTTCCACCCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGGGACGTATTCCGGTTCAAAATcg                                                                  >  2:2371/1‑250 (MQ=255)                                                                                                                                                                  aaTTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCCGCATTCGCAGATCGAACTGCACCACGTTCCACCCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTcaacatcgg                                                                 >  2:316419/1‑243 (MQ=255)                                                                                                                                                                   aTTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATcg                                                                  >  1:710125/1‑248 (MQ=255)                                                                                                                                                                       cccGACCTTTAGCCGGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCCCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCTTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCgcg                                                            <  1:665911/250‑3 (MQ=255)                                                                                                                                                                             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)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCC‑TGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCA  >  NC_000913/64799‑65272

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