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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,160	0	T	C	98.0%	150.6 / ‑4.3	51	G186G (GGT→GGC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/20);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

CGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAC  >  NC_000913/4619920‑4620388                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      cGGTGTCCCGGTTCTGTTTGCGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAAc                                                                                                                                                                                                                             >  2:260691/1‑251 (MQ=255) cGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAAc                                                                                                                                                                                                                             >  1:136303/1‑251 (MQ=255)         cGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGCAGGAACCTTTGGCCTGGATAAACTCTACGa                                                                                                                                                                                                                      >  1:83605/1‑250 (MQ=255)          ggTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAAc                                                                                                                                                                                                                    >  2:80353/1‑251 (MQ=255)          ggTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCAACGAAGGAACCTTTGGCCTGGATAAACTCTACGAAc                                                                                                                                                                                                                    >  1:313040/1‑251 (MQ=255)                gTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCACGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGCTAAACTCTACGAACTGTGCg                                                                                                                                                                                                              >  1:366428/1‑251 (MQ=255)                                    ccGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACGGGTCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGc                                                                                                                                                                                            <  1:74481/249‑1 (MQ=255)                                     cGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCCGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTg                                                                                                                                                                                          >  1:408426/1‑250 (MQ=255)                                             aaagcAGCTTCCAGCAAGAGATGCTGGGTAAACTGGTCGAACGCGCTACTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGACGATTTTCTATACCTCCGCTGACACCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGCACTGAGCGAAATCGCCCGTGAAGAGCTGACCACcggc                                                                                                                                                                                 >  2:51825/5‑251 (MQ=255)                                                 cAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGAACAAcggcggat                                                                                                                                                                             >  1:454879/1‑249 (MQ=255)                                                      tCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGACGGCCACaac                                                                                                                                                                        >  1:185596/1‑250 (MQ=255)                                                              ggtctctgggaaaagagggtagAACGCGCTAACCTGCGGGGTTACCCGGGTAACTGCCACCCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCg                                                                                                                                                                <  2:14294/230‑1 (MQ=255)                                                                     tgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCACGTGAAGAGCTGACCAAAggcggcaaacacaccgggcgtagtat                                                                                                                                                          >  2:332951/1‑230 (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)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      CGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTA‑ACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAC  >  NC_000913/4619920‑4620388

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