Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,385	G→A	A261A (GCG→GCA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,385	0	G	A	98.6%	228.2 / ‑5.0	72	A261A (GCG→GCA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (37/34);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAT  >  NC_000913/4620151‑4620630                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      aacctttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGTCGGCTACAACATCGGTCAAGTTATCGCTCGACCGTTTATCGGAGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGATGCATAAACTGGTTGATGAAAAAAAAGGACAGGGGGTTTCTGTCGGTAAAATTGAAGACATATACGCCAACt                                                                                                                                                                                                                                       >  1:397843/11‑251 (MQ=255)     tttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGATACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAAGACCTGGCTGTTGAGCCGCCAGCACCGACCGGGCTGCCGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCACACTGCgg                                                                                                                                                                                                                                   >  1:208218/7‑251 (MQ=255)       tggccTGGATAAACCCTACGAACTGGGCGAAATCGCCCGTGTAGAGCTGACCAACGGCGGCTACAACAGCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGt                                                                                                                                                                                                                                  <  2:345062/246‑1 (MQ=255)         gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGACGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCACCACCGCCCGTGCTGCAGAAACTGGTTGATGGAAAACACGGTCAGGTGGTTTCTGTAGGTAAAATTTCAGACATCTAAGACAACTCCGGtgtc                                                                                                                                                                                                                               >  1:105480/3‑248 (MQ=255)          ccTGGATTAACTCTCCGAACTGTGCGAAATCGCGCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGATAAAGCCGGTAATGTCCAGCGTACCGGTAACCGTCACGACCTGGTTGTTGAGCCGCCATCACCGACCGTGCTTCAGAAACTGGTTGATGAAAAACACTGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                                                                                                                                                                                              <  1:617375/250‑1 (MQ=255)          ccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                                                                                                                                                                                              <  1:446750/250‑1 (MQ=255)              gtaaaaccctccGAACTGTTCGAAACTGCCAGTGAAGTGCTGCCCAACGGCGGCTACAACTTCGTCCGTTTGATCGCTCGTCCGTTTATCGGTTCCAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                                                                                                                                                                                          <  2:172842/240‑1 (MQ=255)              gATAAACTCTACGAACTGGGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                                                                                                                                                                                          <  1:501122/251‑1 (MQ=255)                 aaaCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaaga                                                                                                                                                                                                                        >  1:513447/1‑248 (MQ=255)                               gtgCGAAATCGCCCGTGAAGAGCTGTCCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTg                                                                                                                                                                                                         <  2:363562/251‑1 (MQ=255)                               gtgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACt                                                                                                                                                                                                          <  1:231169/250‑1 (MQ=255)                                        cGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAc                                                                                                                                                                                                 <  1:168044/250‑1 (MQ=255)                                         gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAcgc                                                                                                                                                                                               >  1:156971/1‑251 (MQ=255)                                                aaGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCTGTGTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGTGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGttt                                                                                                                                                                                        <  2:388344/251‑1 (MQ=255)                                                      tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                                                  >  1:386476/1‑251 (MQ=255)                                                           aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                                                             <  1:294656/251‑1 (MQ=255)                                                                cggcTACAACATCGGCCGTGTTATCGCTCGCCCGTTTATCGGCGACAAAGCCGGTAAATTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACAGCGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGCGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTACCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGAAGCCACCAGCaa                                                                                                                                                                          >  2:293571/1‑249 (MQ=255)                                                                  gcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                                                                                                                                                        >  1:315824/1‑249 (MQ=255)                                                                   cTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                                                                                                                                                        >  1:34968/1‑248 (MQ=255)                                                                       aacatCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTAACGGTAACCGTCACGAACTGGCTGTTGAGCCGCCAGCACCGCCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGATTCTGTCGGCAAACTTGCAGACATCTACGCCACCTGCGGCATCACCACGAAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACACCACCATCAAAGAGCTGa                                                                                                                                                                  >  1:286307/4‑250 (MQ=255)                                                                            cgggcgttttATCGCGCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGACACGTCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAaagaag                                                                                                                                                            <  2:46657/243‑1 (MQ=255)                                                                                    ttATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGAt                                                                                                                                                    >  2:167342/1‑251 (MQ=255)                                                                                      atctcTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCg                                                                                                                                                                                                                                    >  1:343004/5‑169 (MQ=255)                                                                                      atctcTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCg                                                                                                                                                                                                                                    <  2:343004/165‑1 (MQ=255)                                                                                      agcGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCAGCTAGCACCGACCGTGCTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTTGTTTCTTTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAATAAAGTGAAAGCAACTGGCCTGGACGCGTTTTTTGACGCCTCCATCAAAGAGATGAAAGAAGCGGGTGATaa                                                                                                                                                  <  1:605667/249‑1 (MQ=255)                                                                                        cgcctgtcctttaaacggtgACAAAGCCGGTAACTTCCAGCTTACTGGTAACCATCACGCCCTGGCTGTTGATCCGCCAGCACCGATCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCATGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACCTCGGTATCACCCAGAAAGTGAAAGTAACTGGCCTGGACGCGCTGTTTGACGCCACCTTCAAAGAGATGAAATAAGCCGGTGATatca                                                                                                                                                <  2:78078/232‑4 (MQ=255)                                                                                        cGCTCGTCCGTTTATCGGTGACAAAGCCGTTAACTTCCAGCGTACCTGTAACCTTCACGACCTGGCTTGTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCAAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGTGATGAAAGAAGCGGGTGATAac                                                                                                                                                 <  1:400470/250‑1 (MQ=255)                                                                                        cGCTCGTCCGTTTACCGGTATCAAAGCCGGTAGCTTCCAGCGTGCCGGAAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGGAACTGGTTGATGAAAAACACGGTCAGGTGTTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTTGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAaca                                                                                                                                                <  1:285304/251‑1 (MQ=255)                                                                                            cGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCAt                                                                                                                                            >  1:562635/1‑251 (MQ=255)                                                                                              tCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCg                                                                                                                                          >  1:591204/1‑251 (MQ=255)                                                                                                  tttATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCATACATCTACGCCAACTGCGGCATCACCAAGAAAGAGAAGGCAACTGGCCTGGACGCGCTGTTTGACGCCCCCATCAAAGAGATGAAAGAATCggtgaataaccccaacgcctt                                                                                                                                      >  2:256456/1‑232 (MQ=255)                                                                                                    tATCGGTGAAAAAGCCGGAAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTc                                                                                                                                     <  2:156971/250‑1 (MQ=255)                                                                                                       cggtgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCAcc                                                                                                                                  >  1:573123/5‑250 (MQ=255)                                                                                                       cggtgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCAcc                                                                                                                                  >  1:387987/5‑250 (MQ=255)                                                                                                                  ccGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACACCACCCTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGAc                                                                                                                      >  2:261468/1‑251 (MQ=255)                                                                                                                   cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTAATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACt                                                                                                                     <  1:169618/251‑1 (MQ=255)                                                                                                                    ggTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACtt                                                                                                                    >  1:267888/1‑251 (MQ=255)                                                                                                                     gTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTc                                                                                                                   <  1:48938/251‑1 (MQ=255)                                                                                                                           cccctcgtcccGGTAACCGTCACGACCTGGCTGTGGAGCCGCCAGCACCGTCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGGCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGActct                                                                                                             <  2:387987/242‑1 (MQ=255)                                                                                                                             catcGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCACCAGCACCGACCGTGCTGCAGAAACTGGTTGCTGAAAAACACGATCAGGAGGCTTCTGTCGGTAACAATGCAGACATCTACGCCAACTGCGGTATCACCACGAAAGGGAAAGCAACTGGCCTGTACGCGCTGTTTGACGCCACCATCAAAGAGATACAAGAAGCTCGAAATAACACCCACGACTTCCCCCAATTCGTTGACTTCGActctac                                                                                                           >  2:60062/4‑249 (MQ=255)                                                                                                                                cGTACCGGTAACCGTCACGACCGGGCTGTTGAGCCGCCAGCACCGACAGTGCAGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGGCGGAAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCTGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                                                        <  2:265870/251‑1 (MQ=255)                                                                                                                                  tACCGATAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGCAGCGGGTGATAACACCATCGTCTTCACCAACATCGTTGACTTCGCCTCTTCCTggg                                                                                                      >  2:311708/1‑251 (MQ=255)                                                                                                                                                  gACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGt                                                                                      >  1:546951/1‑251 (MQ=255)                                                                                                                                                    ccTGGCTGTTGAGCCGCCAGCACCGAACGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCCGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACGCTTCCTGGGGCCACCGTCGCGACATCg                                                                                    >  2:457009/1‑251 (MQ=255)                                                                                                                                                             tGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTg                                                                                                                                                                                                    <  1:640468/130‑1 (MQ=255)                                                                                                                                                             tGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTg                                                                                                                                                                                                    >  2:640468/1‑130 (MQ=255)                                                                                                                                                                       gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                                                 <  1:399864/251‑1 (MQ=255)                                                                                                                                                                         accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                                                    >  1:189881/1‑246 (MQ=255)                                                                                                                                                                         accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                                               <  1:251500/251‑1 (MQ=255)                                                                                                                                                                         accgaccgTGCTGCAGAAACTGGGTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                                                    <  2:189881/246‑1 (MQ=255)                                                                                                                                                                               cgTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGAGATACCACCATCGTCTTCACCAACTTCGTTGACTACGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGagcttt                                                         >  2:505988/1‑246 (MQ=255)                                                                                                                                                                                     gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                    <  1:256456/250‑1 (MQ=255)                                                                                                                                                                                     gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                    <  1:440735/250‑1 (MQ=255)                                                                                                                                                                                              ggTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTg                                           <  1:277375/250‑1 (MQ=255)                                                                                                                                                                                               gTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGTCCGCCGTCTGcc                                         <  2:90776/251‑1 (MQ=255)                                                                                                                                                                                                  gatgaAAAACACGGACAGGTGGATTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACACTTCCTGGGGCCACCGTCGCTCCGTCGCCGGTTATGCTGCGGGCCTGGAGTTGTACCACCGCCGCCTGCCggg                                      >  2:276264/1‑250 (MQ=255)                                                                                                                                                                                                      aaaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTg                                  >  1:408635/1‑251 (MQ=255)                                                                                                                                                                                                       aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCTGGTGATAACACCATCGTCTTCACCAACGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGa                                 <  1:98925/251‑1 (MQ=255)                                                                                                                                                                                                        aaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGACAGTGAAAGCAACTGGCCTGGCCGAGCTGTTTGACGCCACCAACAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCCCCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCAGCCGCCGTCTGCCGGAGCTGAt                                >  2:601797/1‑251 (MQ=255)                                                                                                                                                                                                                  aGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGGGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTctgctg                      >  1:373386/1‑251 (MQ=255)                                                                                                                                                                                                                   ggaggtTACTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCGGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTTCGGGTCTGGAGCTGGTCGACCGAAGTCTGCCCGTGGTGATGtctttacttc                     >  2:571722/4‑244 (MQ=255)                                                                                                                                                                                                                    gtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcg                    <  1:675156/251‑1 (MQ=255)                                                                                                                                                                                                                    gtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCGGCCGGAGCTGATGTCTCTGCTgcg                    >  1:329055/1‑251 (MQ=255)                                                                                                                                                                                                                     tggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcgc                   >  1:519052/1‑251 (MQ=255)                                                                                                                                                                                                                         ttCTGTAGGTAAAATTGCAGACATATACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGCCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGCTCTGGAGCTGTTCGACCGCAGTCTGACGGAGCTGATGGCTCTGCTGCGCGATg               >  2:292990/1‑251 (MQ=255)                                                                                                                                                                                                                          tCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATg               <  1:311708/250‑1 (MQ=255)                                                                                                                                                                                                                            tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCTGCTGCGCGATgacg            >  2:240837/1‑251 (MQ=255)                                                                                                                                                                                                                                ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt        <  1:9619/251‑1 (MQ=255)                                                                                                                                                                                                                                   aaaaTTGCAGACATCTACGCCAACTTCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCAGCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                                                        <  2:404623/152‑1 (MQ=255)                                                                                                                                                                                                                                   aaaaTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                                                        >  1:404623/1‑152 (MQ=255)                                                                                                                                                                                                                                    aatttGCGGACATCTACGTCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACATCATCGTCTTCACCAACTCCGTTGACTTCGTCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACGGTTCGACCGCCGTCTGCCGGAGCTGATGGCTCGGCCGCGCGATGACAGCATCCTg    >  2:290173/4‑251 (MQ=255)                                                                                                                                                                                                                                      aTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAt  >  1:342506/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      AACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAT  >  NC_000913/4620151‑4620630

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