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	92.9%	75.0 / NA	29	A261A (GCG→GCA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (13/13);  minor base C (2/0);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGA  >  NC_000913/4620147‑4620629                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     aaggaaccttaggccaggATAAACCCTACGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCACCGGCGGCGACAACATCGGTCGGGTTATCGCTCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGcc                                                                                                                                                                                                                                          <  2:100362/235‑1 (MQ=255)            gtgccggATAAACTCTACGAACTGCTCGAAATCGCTCTTTGAGTGCTGACCCACGGCGGCTACAATTTCGGTCGTGTTATCGCTCGTCCGTTTATGGGCGACAAAGCCGGTAACTCCCATCGTACCGGTAACCGTCATGACCTGGCTGTTGAGCCGCGAGCACCGACCGTGCGGCAGAAGCTGGTTGATGAAAAACACGGGCAGGTGGTTTCTGTCGGGAAAATTGCGGACATCTACGCCAACTGCGGTAt                                                                                                                                                                                                                               <  2:165880/246‑1 (MQ=255)             gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTCAAATTGCAGACATCTACGCCAACTGCGGTATc                                                                                                                                                                                                                              >  1:311224/3‑251 (MQ=255)                 ggATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCAcc                                                                                                                                                                                                                           >  1:339207/1‑250 (MQ=255)                    tAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGCCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTTCCGACATCTACGCCCACTGCGGTATCAccc                                                                                                                                                                                                                          >  1:159809/1‑247 (MQ=255)                      aaccctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaagaaa                                                                                                                                                                                                                     <  1:374797/247‑1 (MQ=255)                       aCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAg                                                                                                                                                                                                                    <  1:161927/251‑1 (MQ=255)                            aCGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGaaa                                                                                                                                                                                                               <  1:76168/251‑1 (MQ=255)                                                              cAACGGCGGCTACAACATCTGTCGTGTTATCGCTCGTCCGTTTTTCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGTCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGtcacca                                                                                                                                                                             <  1:69963/251‑1 (MQ=255)                                                                     ggcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCaa                                                                                                                                                                         >  1:53314/1‑248 (MQ=255)                                                                                        ttATCGCTCGTCCGTTTAGCGGTGACAAAGCCGTTAACTTCCAGCGTACCGGTACCCGCCACGCCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGa                                                                                                                                                    <  2:314005/250‑1 (MQ=255)                                                                                             gCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTAACGACCTGGCTGTTGAGCCGCCAGCACCGACAGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTATGTCGGTAAAATTGCAGAAATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAGGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGGGATAacac                                                                                                                                              >  1:170132/1‑251 (MQ=255)                                                                                                cGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACAcc                                                                                                                                             <  1:190136/249‑1 (MQ=255)                                                                                                    cgtttaccggtaaCAAAACCGGTACTTTCCCGCGTACCGGTAACCGTCACGACCTTGTTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCTCCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACAATCGTc                                                                                                                                       <  2:92808/239‑1 (MQ=255)                                                                                                                aaaaaGCCTGTAACTTCCTGGTTACCGGTAACCGCAACGGCCTGGATGGTGAGCCGCCAGTACCGACCGTGCTGCATAAACTGGTTTATGAAAAACACGGTCAGGTGGTGTCTGTCGGTAAAAGTGCAGACATCTACGCCATCTGCGGTATCACCAAGAAAGTGAAAGCAACTGGGCTGTACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATACCACCATAGTCTTCACCAACTTc                                                                                                                           <  2:56004/249‑1 (MQ=255)                                                                                                                       cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACt                                                                                                                    >  1:12333/1‑251 (MQ=255)                                                                                                                                            tAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCAcc                                                                                               >  1:291890/1‑251 (MQ=255)                                                                                                                                               ccGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTc                                                                                            <  1:67843/251‑1 (MQ=255)                                                                                                                                                             cTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGAGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGCCCACCGTCGCGACGTCGCCgtt                                                                              >  2:11820/1‑249 (MQ=255)                                                                                                                                                                 tGAGCCGCCAGCAACGACAGTGCTGCAGAAACTGGTTGATGAAAAACACACTCAGGTGGTTTCTGACGGTACAATTGCAGACAACTACGCCAACTGCGGCATCACCACGCACGAGAAAGCAACTGGCCTGGACGCGCTGTTTGCCGCCACCATCAAAGAGATGAAAGACGCGGGTGATAACACCATCGTCTTAACCAAATTCGTTGACTTCGACATTTCCTGGAGCCACCGTCGCGACGgcgcacggtagt                                                                          >  2:168401/1‑243 (MQ=255)                                                                                                                                                                      cgccAGCACCGACCGTGCTTTAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCtgcg                                                                     <  1:189831/251‑1 (MQ=255)                                                                                                                                                                              ccgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACAGCATCGGCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGa                                                             <  1:192047/251‑1 (MQ=255)                                                                                                                                                                                gaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGagc                                                           >  1:219151/1‑249 (MQ=255)                                                                                                                                                                                  ccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTAGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTg                                                         >  1:19631/1‑251 (MQ=255)                                                                                                                                                                                                                                    gatAAAATTGCCTACATATACGCCAACTGCGGTATCACAAAGAACGAGAATGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGCTAACACCATCGTCTTCACCACCTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGCCGCAGGTTATATTTCGGGTATGGAGCTGGTCGACCACAGTATGCCGGAGCTGACGTCTCTTCTTCGCGATgcagacg        >  2:340567/3‑244 (MQ=255)                                                                                                                                                                                                                                         aaTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCg                                                                                >  1:287743/1‑173 (MQ=255)                                                                                                                                                                                                                                         aaTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCg                                                                                <  2:287743/173‑1 (MQ=255)                                                                                                                                                                                                                                         aaTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTg   >  1:171088/1‑250 (MQ=255)                                                                                                                                                                                                                                         aaTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGGTTGACGCCACCATAAAAGAGATGAAAGACGCGGGGGATAACACCATCGTCTTCACCAACTTCGTTGACTTAGACTCTTCCTGGGGCCTCCGTAGCGACGTCGACGTTTATGCAGCGGGTCTGGAGCTGTTCGACCGACGTCTGCCGGAGCTGATGGCTCTGCTGCGCGagtacgccctacggt  >  2:4978/1‑237 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     AAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGA  >  NC_000913/4620147‑4620629

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