Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,256	C→T	G218G (GGC→GGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,256	0	C	T	94.1%	192.2 / ‑5.4	69	G218G (GGC→GGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (32/32);  minor base G (3/1);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCG  >  NC_000913/4620027‑4620494                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                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                                                                                                                                                           <  2:407556/165‑1 (MQ=255)                                                                                                                                                       aaCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                                                                                           >  1:407556/1‑165 (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)                                                                                                                                                                                                       tCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCAGTCGGTAAAATTGCAGACATCTACGCCAACTGCTGTATCACCAAGCAAGTGAACGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATAAAAGAGATGAaaga                      >  1:323095/1‑250 (MQ=255)                                                                                                                                                                                                            cGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCggg                <  1:46168/251‑1 (MQ=255)                                                                                                                                                                                                                ttatcgcactcccGTTTAGCGGTGACAAAGCCGGGA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTGGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATGGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTgt             <  2:267383/239‑2 (MQ=255)                                                                                                                                                                                                                      cccGTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaagaaa                                                                              <  2:222799/177‑1 (MQ=255)                                                                                                                                                                                                                         gTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTGGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTTGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATc   <  1:113557/251‑1 (MQ=255)                                                                                                                                                                                                                         gTCCGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGGCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATc   <  2:253278/251‑1 (MQ=255)                                                                                                                                                                                                                          cccGTTTATCGGTGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCg  <  2:41073/250‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                TCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCG  >  NC_000913/4620027‑4620494

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