Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,067	A→G	Q155Q (CAA→CAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,067	0	A	G	100.0%	244.4 / NA	75	Q155Q (CAA→CAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/31);  total (44/31)

CGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAG‑AAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCC  >  NC_000913/4619823‑4620312                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ggggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGc                                                                                                                                                                                                                                                  <  2:239729/250‑1 (MQ=255)    gAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGaa                                                                                                                                                                                                                                               <  2:391434/251‑1 (MQ=255)      agggACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCTTCCGGTAAAGATACCCCGTCGGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAaga                                                                                                                                                                                                                                             <  2:535752/249‑1 (MQ=255)        ggACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGGTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGc                                                                                                                                                                                                                                           >  2:137081/1‑251 (MQ=255)                aCGCTGAAGTTATCGGCGAGTACGCATGGGCGCACGAAATGACATACGGAAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCAGGTACGGTCATTCTGGATCAGCTGGGAGAAGAGCACATGaa                                                                                                                                                                                                                                    >  1:175857/1‑250 (MQ=255)                  gCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAc                                                                                                                                                                                                                                 >  1:128811/1‑251 (MQ=255)                       aGTTATCGGCGCGTCCGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCAGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCa                                                                                                                                                                                                                            <  1:358236/251‑1 (MQ=255)                                  cgttcgcttGGGCGCACGAAATTTCATCCGGTAAAGATACCCCGTCTGGTGAATGGTAAATCGCCGGTGTCCCGGTTCTGTTTTAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCGGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                                                                                                                                                                                                                 <  2:613023/243‑1 (MQ=255)                                   gTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCt                                                                                                                                                                                                                >  1:115780/1‑251 (MQ=255)                                       gcaggtgcgcACTAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTTTCCCGGTTCTGTGTGAGTGGGGATATTTCTCCTATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGTGCTAATCTGCCGTGTTACCCCGGTACCTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCTAAGAGCACATTAAAACCGGCAAGCCGATTTTCTATAc                                                                                                                                                                                                            <  2:415634/245‑1 (MQ=255)                                         aTGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCt                                                                                                                                                                                                          >  2:159167/1‑251 (MQ=255)                                               gcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGGCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGATGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACCGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTg                                                                                                                                                                                                    >  2:148390/1‑251 (MQ=255)                                                cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTg                                                                                                                                                                                                    >  1:539168/1‑250 (MQ=255)                                                 aCGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAaga                                                                                                                                                                                                                                             <  2:154386/208‑1 (MQ=255)                                                 aCGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAaga                                                                                                                                                                                                                                             >  1:154386/1‑208 (MQ=255)                                                 aCGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGAc                                                                                                                                                                                                  >  1:454676/1‑251 (MQ=255)                                                  cGAAATGTCATCCGGTAAAGATACCCCGTCTGGCCACTGGGAAATAGCCGGCGTCCCGGTTCTGTTTGAGTGGCGATATTTCCCCGATCACGAACACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTACCTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAATACCGGCAAGCCGACTTTCTATACATCCGCTGACt                                                                                                                                                                                                 >  1:275360/1‑251 (MQ=255)                                                      atatCATCCAGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCACGATTCTGATTGAGTGGGGCTATTTCTCCGATCACGACAACAGCTTCCCGCAAGAGCTGCTGGATAACCTGGGCGAACGCGCTAATCTCCCGGGTTCCCTCGTTAACTACCACTCATCCGGCACGGCCGTCCTGGAACAGCTGGACGAAGCACACATGAAAACCGGCAAGCCGCTTTTCTCTACCTCCGATGCCTCCgg                                                                                                                                                                                             >  2:145338/4‑250 (MQ=255)                                                      aTGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGAGTCCCGGTTCTGTTTGAGTGGGGCTATTCCTCCGATCACGACACCAGCTTCCCGCCAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGACGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgt                                                                                                                                                                                             >  1:314044/1‑251 (MQ=255)                                                        gTCATCCGGTAAAGATACCCCGCCTGGTCACTGGGAAATCGCCGTTGTCCCGGTTCTGTTTGTGAGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCAGTTACGGTCATTCTGGATAAGCTCGGCGAAGAGCACACTAAAACCGGCACACCGATTTTCTATACCTCCGCTGCCTCCgtat                                                                                                                                                                                           >  2:512251/1‑249 (MQ=255)                                                         tCATCCGGTACAGATACCCCGTCTGGTCACTGGGCAATCGCCGGAGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGAGCACGAACACAGCTTCCCGCAAGAGCTGCTGGATAACCTGGTCGAACGAGCTAATCTGCCGGGTTACCGCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGATGGGCGAAGAGCACCTGAAAACCGGCAAGCCGATTTTCTATACCTCCCCTGCCTCCGAGtt                                                                                                                                                                                          >  2:136387/1‑251 (MQ=255)                                                         tCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGtt                                                                                                                                                                                          >  2:402771/1‑251 (MQ=255)                                                                  aaaGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgc                                                                                                                                                                                 >  2:493700/1‑251 (MQ=255)                                                                   aaGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcc                                                                                                                                                                                >  1:217192/1‑251 (MQ=255)                                                                      aTACCCTGTCTGGTCCCTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAATCAGCTTCCCGCAAGAGCGGCTGGATAAACTGGTCGACCGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctgc                                                                                                                                                                             <  2:199560/251‑1 (MQ=255)                                                                        aCCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAAGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctgcc                                                                                                                                                                            >  1:400090/1‑250 (MQ=255)                                                                             gTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGCTTTTCTATACCTCCg                                                                                                                                                                                                       >  2:266106/1‑218 (MQ=255)                                                                             gTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCg                                                                                                                                                                                                       <  1:266106/218‑1 (MQ=255)                                                                                        gggAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAAcctttg                                                                                                                                                            <  1:63729/249‑6 (MQ=255)                                                                                          gaaatcgacGGTGTCCAGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCAGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGACGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGCTTGCCTGCCACGAAG‑GAAcctttggcct                                                                                                                                                        >  2:542928/9‑242 (MQ=255)                                                                                            aatcgCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAAcctttggcctgg                                                                                                                                                      >  1:544973/5‑240 (MQ=255)                                                                                                gCCGTTGTCCCGGTTCTGTTGGAGTGGGGTGATTCCCCCGACCATGAAAACAGCTCCCCGCAAGTGCTGCTGGACAAACTGTCCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAAccgttggcctggataa                                                                                                                                                  <  2:134289/251‑16 (MQ=255)                                                                                                gCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATa                                                                                                                                                   >  2:158806/1‑250 (MQ=255)                                                                                                  cGGTGTACCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACTGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAAc                                                                                                                                                >  1:480110/1‑251 (MQ=255)                                                                                                     tgtCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAActct                                                                                                                                             >  2:139523/1‑251 (MQ=255)                                                                                                     tgtCCCGGATCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGCTTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAACctct                                                                                                                                             >  2:287228/1‑251 (MQ=255)                                                                                                              tCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACtgt                                                                                                                                    >  1:546333/1‑251 (MQ=255)                                                                                                                 gTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACtgt                                                                                                                                    <  1:343113/248‑1 (MQ=255)                                                                                                                   ttGAGTGGGGATATTTCCCCGATCACGAAACCAGCTTCCCGCAAGAGCTGCTGGATAAACGGGTCGAACGCGCAAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGa                                                                                                                                <  1:118039/250‑1 (MQ=255)                                                                                                                     gAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACGCTACGAACTGTGCGaaa                                                                                                                              <  1:169018/250‑1 (MQ=255)                                                                                                                      aGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATc                                                                                                                            <  1:158806/251‑1 (MQ=255)                                                                                                                          gggATATTTCTCCTATCACGAAAACATCTTCCCGGAAAGGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGTTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCTAAATCGccc                                                                                                                        <  2:310014/251‑1 (MQ=255)                                                                                                                           ggATATTTCTCCGATCAAGAAAAAAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTATGGCCTGGATAAACTCTACGAACTGTGCGAACTCGCCCg                                                                                                                       >  2:129480/1‑251 (MQ=255)                                                                                                                               atTTCTCCGATCACGAACACAGCTTCCCGCAAGAACTTCTGGATAAACTGGTCGAACGCTCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCAGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAGAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACAAAG‑AAACCTTTGGCCTGGATACACTCGACAAACTGTGCGAAATCACCCGTGac                                                                                                                   >  2:554644/1‑250 (MQ=255)                                                                                                                                  tctcCGATCACGAAACCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGCAATCGCCCGTGAag                                                                                                                  >  2:420927/1‑249 (MQ=255)                                                                                                                                      cgatcaagaccacaacTTCCCGCAAGAGCTGCTGGATAAACTGGTAGAAAGCGCTAATCTGCCGGGTAACCTCGGTAACTGCCACTATTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACAGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGACG‑GAACCTTTGGCCTGGCTAAAGTCTACGAACTGTGCGAAATCGCCCGTGACGAGcggc                                                                                                            >  2:91607/16‑248 (MQ=255)                                                                                                                                            cGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGTTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGTTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAATAGCTGa                                                                                                            <  1:513637/245‑1 (MQ=255)                                                                                                                                             gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGCAACTGCCACTCTTCCGGCACGGTCATTCCGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGCACTGTGCGAAATCGCCAGTGAAGAGCTGACCAAc                                                                                                       >  1:300382/1‑249 (MQ=255)                                                                                                                                               aaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcggcg                                                                                                   <  1:159167/251‑1 (MQ=255)                                                                                                                                                   aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTa                                                                                               <  1:19774/251‑1 (MQ=255)                                                                                                                                                         ccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACaacac                                                                                         >  2:508074/1‑248 (MQ=255)                                                                                                                                                          cGCAAGAGCTGCTGGATAAACTGGGCAACCGCGCTAATCTGCCGGGTTGCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATc                                                                                        <  1:508074/251‑1 (MQ=255)                                                                                                                                                                agCTGCGGGATAAACTGGTCGAACGCGCTAACCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGATCACATGAAAACCGGCAAGCCGATTTTCGATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCgt                                                                                  <  1:97798/251‑1 (MQ=255)                                                                                                                                                                 gctgcggtaTAAACTGGGCGAACGCGCTAATCGGACGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACGCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCgtg                                                                                 <  2:180710/243‑1 (MQ=255)                                                                                                                                                                  ctgctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGAACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTCTACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGACCGTGAAGAGCTGACCAACCGCGGCTATAACATCAGTCgtgt                                                                                >  1:56988/1‑251 (MQ=255)                                                                                                                                                                   ggctgaaaaaaCTGGTCGACCGCGCTAATCTGCCGGGTTCCCTCGGTACCTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGtt                                                                               <  2:366976/243‑1 (MQ=255)                                                                                                                                                                     ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCtcccacagcggtcgggtttt                                                                             >  2:497009/1‑232 (MQ=255)                                                                                                                                                                       ggATAACCCGGTCCAACGCGCAAACTTGCCGGGTTCCCGCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCCCCGCGGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTGGGCCTGGATAAACTCTACGAACTATGTGAAATCGCCCGTGAAGAGCTGACCAACTGCGGCTACAACATCGGTCGTGTTAt                                                                             <  2:630243/249‑1 (MQ=255)                                                                                                                                                                        gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATAc                                                                                                                                                                                                            <  1:422072/122‑1 (MQ=255)                                                                                                                                                                        gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGACACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATAc                                                                                                                                                                                                            >  2:422072/1‑122 (MQ=255)                                                                                                                                                                         aTAAACTGGTCGAACGCTCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCt                                                                         <  1:285786/251‑1 (MQ=255)                                                                                                                                                                          tAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTc                                                                        >  1:156083/1‑251 (MQ=255)                                                                                                                                                                           aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGGTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTc                                                                        <  1:252287/250‑1 (MQ=255)                                                                                                                                                                           aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCg                                                                       >  1:130787/1‑251 (MQ=255)                                                                                                                                                                           aaaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTc                                                                        <  1:452505/250‑1 (MQ=255)                                                                                                                                                                                    gAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGACGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GACCCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGATCGGGTAATCGCTCGTCCGTTTAt                                                              >  1:56211/1‑251 (MQ=255)                                                                                                                                                                                     aaCGCGCTACTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATc                                                             >  1:173186/1‑251 (MQ=255)                                                                                                                                                                                      aCGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGCTCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGCCGTGTTATCGCCAGTCCGTTTATCg                                                            >  2:614451/1‑251 (MQ=255)                                                                                                                                                                                          gcTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCggtga                                                        <  1:493700/251‑4 (MQ=255)                                                                                                                                                                                            tAATCTGCCGGGTTAACTCGGAACCTGCCACTCTTCCGGTACGGTCATTCTGGATAAGCTGGGAGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTCTACATCAGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAACtttggcc                                                                                                                                                          >  2:82943/1‑148 (MQ=255)                                                                                                                                                                                                          aCCTCGGTAACTGCCACTCTTCCGGTACGGACATTCTGGATCAGCTGGGCGAAGAGCACATGCAAACCGGCAAGCCGATTTTCTATACCTCCGCAGACCCCGTGTTCCAGCTTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACCTCGGTCGCGTTATCGCTCGACCCGTTATCGGCGCCAAATCCGGTTACTTc                                        >  2:278604/1‑251 (MQ=255)                                                                                                                                                                                                                            ctTCCGGTACGGTCATTTGGGCTCAGCTGGGCGAAGAGCACATGAAATCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTTGTCCAGAGTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGTCCTACGGCGGCTACAACATCGGTCTTGTTATCGCTCGTCCGTTTATCGGTGACAAGGCCGGTAACTTCCATCGTACCGGTatcc                        <  2:360515/249‑4 (MQ=255)                                                                                                                                                                                                                              tCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGt                      <  1:137081/249‑1 (MQ=255)                                                                                                                                                                                                                               acGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAAAATCGGTCGTGTTATAGCTCGTCCGTTTTTCGGTGACCAATCCGGGAACTTCCAGCGTACCGGTAAACGTCa                    >  2:247827/2‑250 (MQ=255)                                                                                                                                                                                                                                                gATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAG‑GAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCCGCGTCCCGGTAACCGTCACGACCTGGCTGTTGAgcc  >  1:270021/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAG‑AAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCC  >  NC_000913/4619823‑4620312

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