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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,223	0	T	C	97.0%	185.1 / NA	69	N207N (AAT→AAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (35/30);  minor base A (2/0);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGA  >  NC_000913/4619981‑4620465                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    ggtctctgggaaaagagggtagAACGCGCTAACCTGCGGGGTTACCCGGGTAACTGCCACCCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCg                                                                                                                                                                                                                                             <  2:14294/230‑1 (MQ=255)             aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCg                                                                                                                                                                                                                                 >  1:384306/1‑251 (MQ=255)                      gAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTAt                                                                                                                                                                                                                        <  1:98987/251‑1 (MQ=255)                         cgcgCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGGTATCGCTCGTCCGTTTATCgg                                                                                                                                                                                                                     >  1:46873/1‑251 (MQ=255)                         ccgctaaacttGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAAAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATGGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCgg                                                                                                                                                                                                                     <  2:249613/241‑1 (MQ=255)                           cgcTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCCCGTCCGTTTATCggtg                                                                                                                                                                                                                   >  2:335448/1‑249 (MQ=255)                                    gCCGGGTTACCTCGGTAACTGCCACTCTTCCGGAACGGACATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTCCAAGATCGGTCGTGTTATCGCTCGTCCGTTTagaggtggcaaagccg                                                                                                                                                                                                          >  1:156041/1‑236 (MQ=255)                                               tCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATAGCCCGTGAAGAGCTGACCAACGGCGGCCACAACATCAGTCGTGGTATCGCTCGTCCGTTTTTGGGGGGCAAAGCCTGTA‑ACATCCa                                                                                                                                                                                               >  2:83539/1‑250 (MQ=255)                                                                 ccGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCAc                                                                                                                                                                            >  1:323539/1‑251 (MQ=255)                                                                    gTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACggc                                                                                                                                                                         >  1:86987/1‑249 (MQ=255)                                                                      aCGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAGTTTCTATCCCGCCGCTGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACAGCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCt                                                                                                                                                                       <  1:48344/251‑1 (MQ=255)                                                                       cggtcattcggtaccagcggggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTg                                                                                                                                                                      <  2:12085/232‑1 (MQ=255)                                                                         gTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTCCCGGAAACAGTCACGACCTGGc                                                                                                                                                                    >  1:296977/1‑251 (MQ=255)                                                                         gTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAAAACATCGGGCGGGTTATGCCTCGTCCGGTTATTGGTGTCAAAACCGGTACATTTCC‑TCGTACCGGTGACCGTCGCGAcccgg                                                                                                                                                                     >  2:298971/1‑247 (MQ=255)                                                                           cATTCTGGAGCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGGGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTg                                                                                                                                                                  <  2:194744/251‑1 (MQ=255)                                                                            aTTCTGGATCAACGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATGCCGCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGt                                                                                                                                                                 <  1:12073/251‑1 (MQ=255)                                                                              tCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAAACTCGGTCGTGTTATCACTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACAGGGAACCGTCACGACCTGGCTGttt                                                                                                                                                               >  2:327556/1‑250 (MQ=255)                                                                                tGGATCAGCTGGGCGAAGAGCACATGTAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGa                                                                                                                                                              <  2:323539/250‑1 (MQ=255)                                                                                 ggATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATTGGTCGTGTTATCGCTCGTCCGTTTTTCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCGGTTTAgc                                                                                                                                                            >  2:451566/1‑251 (MQ=255)                                                                                      agcTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAgccgccc                                                                                                                                                       >  1:15262/3‑250 (MQ=255)                                                                                         tGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCGACAAAATCGGGCGTGTTATCGACCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGGACCGGTAAACGTCACGACCTGGCTGGTGAgcagaacga                                                                                                                                                     >  2:321383/1‑243 (MQ=255)                                                                                          gggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCTAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCa                                                                                                                                                    <  1:68559/250‑1 (MQ=255)                                                                                          gggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACCACATAGGACGTGTTATCGCTCGTACGTTTTTCGGTGACCAAGCCGGTA‑ACTTCCCACGGACCGGGAACCCTCACGCCCTGGCTGTTGAgccccccgccc                                                                                                                                                   >  2:383607/1‑243 (MQ=255)                                                                                                agaacacaTGACAAACGGCAAGCCGATTTTGTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGAAGGGACCTTTGGCCTGGCTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACACCATCGGGCGTGTTATCGCGCGTCCGGTTATAGGGGTCAAAGCCGGTA‑AATTCCACAGTACCGGTGACCGCCACGAACTGGCTGTTGAgacgccgccaccgcccc                                                                                                                                             >  2:55017/5‑240 (MQ=255)                                                                                                  agCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGTTTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtg                                                                                                                                           <  2:384306/251‑1 (MQ=255)                                                                                                  agCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtg                                                                                                                                           >  1:353786/1‑251 (MQ=255)                                                                                                      caTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctgc                                                                                                                                       <  1:371709/251‑1 (MQ=255)                                                                                                       aTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCa                                                                                                                                      <  1:15157/251‑1 (MQ=255)                                                                                                        agAAAACCGCCAAGCCGATGTTCGTTCCCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACACCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAg                                                                                                                                     <  2:246306/250‑1 (MQ=255)                                                                                                           aaaCCGGCAAGCCGGTTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGCACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGAGGTATTGCCCGTCCGTTTCTCGGTGCCACAACCGGTA‑ACTTCCAGCGTACCGGTAAACGTGAAGACCTGGCTGTTGAGTCGCCCGCACCGCACGTGCTGCCGaa                                                                                                                                   >  2:222263/1‑250 (MQ=255)                                                                                                             aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGCTTGCCTGCCTCGACGGAACCTTTGGCCTGGATACAATCTACGAACTGTGCGAACTAGCCCGTGAAGCGCTGACCAACGCCGGCCACAACATCGGTCGAGCTATAGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAACGTACCCGTAACCGTCACGACCTGGCTGTTTAGCCGCCAGCACCGACCGTGCTGCAGAAACt                                                                                                                                >  1:23302/1‑251 (MQ=255)                                                                                                                           tttcgatccccccGCTGATTCCGTGGTCCAGAGTTCCTGCCCTGAAGGAACCTTGGGCCTTGATAAACTCTACGAACTGTTCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTGATCGCTCGTGCTTTTATCGGTGACAAAGCCGGTA‑ATGCCCAGCTTACCGGTAACCGTCACGGCTTGGCTGTTGAGCCGCCAGCACCGATCGTGCTGCAGAAACTGGGTTATGaaaaa                                                                                                                   <  2:477720/239‑1 (MQ=255)                                                                                                                                   ccTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAg                                                                                                          <  1:456986/251‑1 (MQ=255)                                                                                                                                    cgccgcggcctctgtgtTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGACCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGGACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTGGATGAAAAACACGGTCAgg                                                                                                         <  2:434737/238‑1 (MQ=255)                                                                                                                                           gACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAATTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTc                                                                                                  <  1:421057/251‑1 (MQ=255)                                                                                                                                           gACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCCGCaccgaccg                                                                                                                                             >  1:400452/1‑208 (MQ=255)                                                                                                                                           gACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccgaccg                                                                                                                                             <  2:400452/208‑1 (MQ=255)                                                                                                                                                     cccAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGAAAAACTCTACGAACTTGGCGAAATCGCCCGGGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCGCGTCCGTTTAGCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                         <  2:146323/249‑1 (MQ=255)                                                                                                                                                              cctgccACGAAGGAACCTTTGGCCTGGATAAACTCTTCGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGAc                                                                               <  1:383607/251‑1 (MQ=255)                                                                                                                                                                   cacgaaggaacctttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTa                                                                          <  1:451566/233‑1 (MQ=255)                                                                                                                                                                          gaacctttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGGTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCCGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCCCCCGCACCGACCGTGCTGCGGAAACTGGTTGATGAAAAACACGGTCCGGGGGTTTCTGTCGGTGACATTGCAGACATCTCCGCCAAc                                                                   >  2:150484/12‑251 (MQ=255)                                                                                                                                                                                    ccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                         <  1:471687/250‑1 (MQ=255)                                                                                                                                                                                      tGGATAAACTCTACGAACTGTGCGCAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACCCCATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCATCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCCCCAGCACCCACCCTTCTGCCGAAACCGGTTTATGTAAAACCCGGTCCGGTGGTTTCTGTCGGTAAAATTGCAGACATTTCCCCCAACTGCGGTATcccc                                                       >  2:242039/1‑248 (MQ=255)                                                                                                                                                                                      tGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGGTGATGGAAAAAACGGTCAGGGGGTTTCTGTCGGTAAAATTGCAGACATCTACCCCCACTTCGGTATCAcc                                                       >  2:312130/1‑251 (MQ=255)                                                                                                                                                                                           aaaCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGAACGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGGCGGTAAAATTGCAGACATCTACGCCAACTGCGGAATCAccca                                                     >  2:61316/1‑246 (MQ=255)                                                                                                                                                                                            aaCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTTGTGGTATCACTCGTCCGTTTATCGGTTCCAAAAACGTTA‑ACTTCCATCGTACCGGAAACCGTCACGACCTGGGTGTTCTGCCGCCCACACCGGCCGGGCTGCCGAAACCGGTTGATGCAAAAAACGGTCAGGGGGATTCTGTCGGTAGAATTGCAGACCTCTACGCCTACTGccggctgatcaagga                                                  >  1:240846/1‑236 (MQ=255)                                                                                                                                                                                               tctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGCCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGGATCACCAAGAAAGTg                                              >  1:111667/1‑251 (MQ=255)                                                                                                                                                                                                ctACGAACTGGGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTACCGGTGACAAAGCCGGTA‑ACTGCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGt                                               <  1:298971/249‑1 (MQ=255)                                                                                                                                                                                                  ccGAACTGTGCGAAATCGCCCGTGAAGAGCTGTCCACCGGCGGCTACAACATCGGTCGTGTTATCGCTCTTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGaaa                                           <  2:146232/250‑1 (MQ=255)                                                                                                                                                                                                   cGAACGGTGCGAAATCGCCCGTTAAGAGCTGCCCAACGGCGGCTACAACATCGTTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCTGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAg                                          <  1:280328/251‑1 (MQ=255)                                                                                                                                                                                                     aaCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                                                              >  1:407556/1‑165 (MQ=255)                                                                                                                                                                                                     aaCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                                                              <  2:407556/165‑1 (MQ=255)                                                                                                                                                                                                       cTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGGTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGGCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGAGGTTTCTGTCGGGAAAATTGCCAACATCTACGCCAACTGAGGTGTCCCCAAGAAAGTGaaatcaac                                      >  1:129709/1‑246 (MQ=255)                                                                                                                                                                                                          tgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTgg                                   >  1:35219/1‑251 (MQ=255)                                                                                                                                                                                                           gctaaaTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTTGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGc                                  <  1:203260/248‑1 (MQ=255)                                                                                                                                                                                                             gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGCCCGGGCTGCAGCAACTGGTTGATGAAAAACCCGGTCCGGTGGTTTCTGTCGGTAAAATTGCAGACATCTTCCCCAACTGCGGGATCACCAAGAAAGTGAAAGCAACTGGCCt                                >  1:385645/1‑251 (MQ=255)                                                                                                                                                                                                                   gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAcgc                          >  1:158733/1‑251 (MQ=255)                                                                                                                                                                                                                   gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCCCCAGCACCGACCGTGCTGCAGAAAATGGTTGATGAAAAACACGGTCCGGGGGTTTCTGTCGGTAAAAATGCAGACATCTACCCCAACTGCGGTTTCACCCAGAAAGTGAAAGCAACTGGCCTGGaagc                          >  2:273171/1‑248 (MQ=255)                                                                                                                                                                                                                     ccGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACgcgc                        >  1:41073/1‑251 (MQ=255)                                                                                                                                                                                                                        tGAAGATCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACCAAGCCGGTA‑ACTTCCAGCGTACCGGGAACCTTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGGAAAACACGGGCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCCAGACAGAGAAAGCCACTTGCCTGGGCGCCCTgg                     >  2:135919/1‑250 (MQ=255)                                                                                                                                                                                                                        tGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCCCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGCGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGt                     >  1:109635/1‑251 (MQ=255)                                                                                                                                                                                                                        ggAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGt                     <  2:158733/250‑1 (MQ=255)                                                                                                                                                                                                                         gAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGtt                    <  1:260691/251‑1 (MQ=255)                                                                                                                                                                                                                          aaGACCTGACCAACGGCGGCTACAACACCGGTCGTGTTATCGCTCGCCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTAGCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGttt                   <  2:385645/251‑1 (MQ=255)                                                                                                                                                                                                                             agCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGAc                <  1:61316/251‑1 (MQ=255)                                                                                                                                                                                                                               cTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGc              >  1:253278/1‑251 (MQ=255)                                                                                                                                                                                                                                 gACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTTATGAAAAACAAGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCAACAAGAAAGTGAACGCAACTGGCCTGGACGCGCTGTTTGGCGccc            >  1:267383/1‑250 (MQ=255)                                                                                                                                                                                                                                           ggcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGCCAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCaa     >  1:186092/1‑248 (MQ=255)                                                                                                                                                                                                                                            gcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCCGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACCGCGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGCCCTCTACGCCAACTGCGGTATCACCACGAAAGTGAAAGCAACTGGCCTGGAGGCGCCGTTTGACGCCCCCATCAAaga  >  2:325075/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGA  >  NC_000913/4619981‑4620465

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