Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,251,586	G→A	Q249Q (CAG→CAA)	nfo →	endonuclease IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,251,586	0	G	A	98.4%	206.2 / NA	63	Q249Q (CAG→CAA)	nfo	endonuclease IV
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/39);  minor base C (1/0);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCGTGGAAGATAAATCCCGCGTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGAAAAGGAGATCCTTGCAATTTTACGGCGTAACCCGCTGATTCAGCAGAACGAAATTGCGGACATGCTGCAAATCAGCCGTTCGCGCGTTGCGGCGCATATTATGGATTTAA  >  NC_000913/2251339‑2251823                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               gCGTGGAAGTTAAATCCCGCTCCGGCGTCTGCATGCATCCCCGCCATGCTTTCGCTGCCGGGTATGATGTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAATTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTGGACCGCCATCATAGCCGCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCGGGATCATGCAAgac                                                                                                                                                                                                                                             <  2:244478/251‑1 (MQ=255) gCGTGGAAGATAAATCCCGCGTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgac                                                                                                                                                                                                                                             <  1:283068/251‑1 (MQ=255)    tGGAAGATAAATCCCGCGTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGAAATATCGGTCATGATGCGTTCCGCTGGATCATGCACgacgac                                                                                                                                                                                                                                          >  2:49836/1‑251 (MQ=255)        aGATAAATCCCGCGTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGtt                                                                                                                                                                                                                                      >  2:448772/1‑251 (MQ=255)          aTAAATCCCGCGTCGGCGGCTGTATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTGGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCGCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                                                                                                                                                                                                    <  2:391935/251‑1 (MQ=255)             aaTCCCGCGTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGACTGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCAGTATGCGTTCCGCTGGATCATGCAAGACGACCGTTTccac                                                                                                                                                                                                                                 >  2:526565/1‑248 (MQ=255)                cccGCGTCGACGTCTGCATTGATACATGACATGCTTTCGCTGCCGGGTATGATTTGCGCACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCCGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATTCAACACGACCGTTTCGACGGc                                                                                                                                                                                                                              >  2:555344/1‑251 (MQ=255)                     gTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTcc                                                                                                                                                                                                                         >  1:354825/1‑251 (MQ=255)                        ggCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGACGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCt                                                                                                                                                                                                                      >  2:174096/1‑251 (MQ=255)                           gTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCAACTATGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCACTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCTACGGCATTCCGCTGAt                                                                                                                                                                                                                   >  2:450846/1‑251 (MQ=255)                               gCATTGATCCCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                                                                                                                                                                                               <  1:163466/251‑1 (MQ=255)                                 attgactccTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTTCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGTATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                                                                                                                                                                                             <  1:164028/244‑1 (MQ=255)                                  ttGATCCCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                                                                                                                                                                                             <  1:174096/250‑1 (MQ=255)                                     aTACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATGCCGCTGATCCTCGAAAcc                                                                                                                                                                                                         <  2:157393/251‑1 (MQ=255)                                        ccTGCCATGCTTTCGCTGTCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCCTTCCGCTGATCCTCGAAACCATc                                                                                                                                                                                                      >  2:481745/1‑251 (MQ=255)                                                 gtttagctgccggggaggaggtGCGTACTCCNGCCGAATGCGAGAAAACAGTCGGGGATTTTGCCCGTACGGGCGGCTTTAAGTATCTTCGCGGGATGCCCCTTAACGATGCGAAATGCCCCTTTGGCAGCGGCGTTGACCGCCATCATAGCCGCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAATCCATCAACCCGGat                                                                                                                                                                                             <  2:272170/230‑1 (MQ=255)                                                  tttCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                                                                                                                                                                                               <  2:148373/232‑1 (MQ=255)                                                  tttCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                                                                                                                                                                                               >  1:148373/1‑232 (MQ=255)                                                   ttctgtgcggtgTATTATTTGCGTCCTCCGGCCGAATGCGTTAAAACATTCGCGGATTTTGCCCTTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAGCGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGat                                                                                                                                                                                             <  2:419695/238‑1 (MQ=255)                                                      gCTGCCGGGTATGATTTTCGTACTCCAACCGAATGCGGGAATACCTTCGCGGATTGTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCCCCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTg                                                                                                                                                                                        <  2:309099/251‑1 (MQ=255)                                                         gCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTgg                                                                                                                                                                                       <  1:140001/249‑1 (MQ=255)                                                         gCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGc                                                                                                                                                                                     <  1:448772/251‑1 (MQ=255)                                                           cGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCa‑g                                                                                                                                                                                  <  2:94914/251‑1 (MQ=255)                                                                     ttGCGTACTCCAGCCGAATGCGAGAAAACATTCGTGGCTTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCCATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGCCGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATATGGGCA‑GCAGAGAtagc                                                                                                                                                                        >  2:400356/1‑248 (MQ=255)                                                                         gTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTgg                                                                                                                                                                    <  1:263713/251‑1 (MQ=255)                                                                           aCTCCCGCCGAACGCGCGAAATCATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGTGCTGCAGGCACCTTAACGATTCGAAAAGCACCTTTGGCAGCCGCGTTGTCCGTCATAATAGCCTCGGTGAAGGCAATTGCGGTCATGATGCGTTCCGCTGGATCATTCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGTACCCATCAACCCGTATATCTGGGCA‑GAAGTGATCGCCTGGc                                                                                                                                                                   <  2:337921/250‑1 (MQ=255)                                                                           aCTCAAGCCGAATGCGAGAAGACATTCGCGTATTTTGCCCGGATTGTCGGTTGTAATTAGCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTTGCAGCCGCGTTGACCGCAATCATAGCCTCGGTGAAGTCAATATCGGTCATGATGCGTTCCTCTGGAGCATGCAAGACGACCGTTTCGTCGGCATTCCGCTTATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCt                                                                                                                                                                  <  2:515435/251‑1 (MQ=255)                                                                             tCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAACGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTg                                                                                                                                                                 >  1:504022/1‑250 (MQ=255)                                                                                     aaTGCGAGAAAGCATGCACGGATTTTGCCCGTACTGTCGGCTTTAAGTCTCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACATTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGGGACGGAAATAACGGTCATGATACGTTACGCTGGATCATGCAAAACGACCGTTTCGAAGGCATTCCGTTGATCCTCGAACCCATCAAACCGGAATTCAGGGCAATAACA‑ATCGCCTGGCTAAAAGCGca                                                                                                                                                         >  2:366757/1‑250 (MQ=255)                                                                                     aaTGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGAGGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGATGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGAC‑GCCGAGATCGCCTGACTGAAAgagcc                                                                                                                                                         >  2:412524/1‑246 (MQ=255)                                                                                             aaaaCATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGAGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATTATAGCCTAGGCGAAGGCAATATAGGTCATGATGCGTTCCGCTGGAGCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGCTATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTg                                                                                                                                                 >  2:533880/1‑250 (MQ=255)                                                                                                      gcgGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGt                                                                                                                                       >  1:6054/1‑251 (MQ=255)                                                                                                          aTTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGcc                                                                                                                                   <  1:341706/251‑1 (MQ=255)                                                                                                                        gTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCAAAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCAGGCAAGACGACCGTTTCGACGGCATTCCGCTGAGCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAAAGATCGCCTGGCTAAAAGCGCAACAACATGAAAAAGCGGTAGCCTGAAGAgcaagtcc                                                                                                                     >  2:311752‑M1/1‑243 (MQ=255)                                                                                                                        gTCGGCTTTAAGTATCGGCGCGGGATGCACCTTAGCGATGCGAAAAGCTCCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGATGTCAATATCGGTCATGATGCGTGCCGCTGGATCAGGCAAGACGTCCGTTTCGACGGCATGCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAATGCGGTCGCCTGAAGAgcacga                                                                                                                       <  2:249115‑M1/249‑7 (MQ=255)                                                                                                                                tAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCCCCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggta                                                                                                              <  1:206521‑M1/250‑16 (MQ=255)                                                                                                                                 aaGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaa                                                                                                            <  1:190642‑M1/251‑18 (MQ=255)                                                                                                                                      ccggcgcgggaggcgcagtaacgaggcgaaaagcaccCTTGGACGACGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAACAGCGGTCATGATGCGTGACGCTGGATCATGCAAGACGACCGTTGCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCGGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaatt                                                                                                       <  2:107444‑M1/233‑23 (MQ=255)                                                                                                                                            cgggatgcaaaaaaaagatgcgaaaagcaccttTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGAA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaac                                                                                                 >  1:141677‑M1/17‑223 (MQ=255)                                                                                                                                              ggatgccaccaaacgatgcgaaaaGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccc                                                                                               <  2:48672‑M1/240‑31 (MQ=255)                                                                                                                                              ggATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccc                                                                                               >  2:575713‑M1/1‑221 (MQ=255)                                                                                                                                              ggATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccc                                                                                               <  1:403224‑M1/251‑31 (MQ=255)                                                                                                                                                     ccTTAACGATGCGAATACTACTGTTGGCAGCNGCGTTGCCCGCCATAATAGACTCGGGGAAGGCAATATCGGTCATGCTGCGTTCCGCTGGATCATGCAAGACGACCGTGTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGACTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaagtaccacccctgtcaat                                                                                        <  2:294197‑M1/251‑38 (MQ=255)                                                                                                                                                        tAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatcc                                                                                      <  1:533880‑M1/250‑40 (MQ=255)                                                                                                                                                           cGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaa                                                                                  <  1:400356‑M1/251‑44 (MQ=255)                                                                                                                                                                   aaaGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatat                                                                            >  1:224189‑M1/1‑200 (MQ=255)                                                                                                                                                                                cggccGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagac                                                             <  2:219150‑M1/249‑65 (MQ=255)                                                                                                                                                                                      cgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatacggcaaagacggttta                                                       <  1:188088‑M1/251‑71 (MQ=255)                                                                                                                                                                                                   catAGCCTCGGTGAAAGCAATATCGGTCATGATTCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGGTCCGCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCTGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccagccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttc                                            <  2:504022‑M1/249‑82 (MQ=255)                                                                                                                                                                                                    attgcctCGGGGAAGGCAATATCGGTCATGATGCGTGCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcata                                         <  2:141677‑M1/248‑85 (MQ=255)                                                                                                                                                                                                      aGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatt                                       >  1:320957‑M1/1‑165 (MQ=255)                                                                                                                                                                                                           cGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGCTATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaa                                   >  2:323386‑M1/1‑160 (MQ=255)                                                                                                                                                                                                           cGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGACAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcacaaacggtttaaacatgtgttcatatttaaac                                  >  2:112046‑M1/1‑160 (MQ=255)                                                                                                                                                                                                            ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaaca                                 <  1:237716‑M1/251‑93 (MQ=255)                                                                                                                                                                                                                 aGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaacatgtgc                            <  1:302420‑M1/251‑98 (MQ=255)                                                                                                                                                                                                                  ggCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaacatgtgct                           <  1:323386‑M1/251‑99 (MQ=255)                                                                                                                                                                                                                   gCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatcaacaagaatatagatactgcaaagacggtttaaacatgtgttaatatttaaacatgtgctc                          >  2:312812‑M1/1‑152 (MQ=255)                                                                                                                                                                                                                          cGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatctttaaacatgtgctcatagtat                   >  2:302244‑M1/1‑145 (MQ=255)                                                                                                                                                                                                                          cGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaacatgtgctcatagtat                   <  2:350438‑M1/251‑107 (MQ=255)                                                                                                                                                                                                                                 gatgcttgccgctggatcATGCAAGACGACCGTTTCGCCTGCATTCCGCTGAGCCTCGAATCCATCAACCCGGATATCTGTGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaatcatgtgttcacatttaaacatgtgctcatagtattcatgtg            <  1:98030‑M1/243‑114 (MQ=255)                                                                                                                                                                                                                                  atgCGGTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaacatgtgctcatagtatgcatgttt           <  1:369477‑M1/251‑115 (MQ=255)                                                                                                                                                                                                                                     cGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGAgcacgtcctctggtaaaaaattaccaaccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaacatgtgctcatagtatgcatgttttac        <  1:555344‑M1/251‑118 (MQ=255)                                                                                                                                                                                                                                           gCTGGATAATGCAAGCCGCCCGTTGCGACGGCATTCCGCTGATCCGCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAACCTGAAGAgcacgtcctctggtaaaaaattaccatccctgtcaatccacaagaatatagatactgcaaagacggtttaaacatgtgttcatatttaaacatgtgctcatagtatgcatgttttacgtaaag  <  2:356303‑M1/251‑124 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GCGTGGAAGATAAATCCCGCGTCGGCGTCTGCATTGATACCTGCCATGCTTTCGCTGCCGGGTATGATTTGCGTACTCCAGCCGAATGCGAGAAAACATTCGCGGATTTTGCCCGTACTGTCGGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCA‑GAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGAAAAGGAGATCCTTGCAATTTTACGGCGTAACCCGCTGATTCAGCAGAACGAAATTGCGGACATGCTGCAAATCAGCCGTTCGCGCGTTGCGGCGCATATTATGGATTTAA  >  NC_000913/2251339‑2251823

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