Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,934,882	G→A	28.6%	A369T (GCC→ACC)	rbsA →	D‑ribose ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,934,882	0	G	A	28.6%	17.9 / 6.5	14	A369T (GCC→ACC)	rbsA	D‑ribose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base A (2/2);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCACCGCAGGATGGCCTGGCAAACGGCATTGTGTATATCTCCGAAGACCGTAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACA  >  NC_000913/3934733‑3935013                                                                                                                                                      |                                                                                                                                    gTTCACCGCAGGATGGCCTGGCAAACGGCATTGTGTATATCTCCGAAGACCGTAAACGTGACGGTGTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATg                                                                                                                                     >  2:3865/1‑150 (MQ=255)          aGGATGGCCTGGCAAACGGCATTGTGTATATCTCCGAAGACCGTAATCGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAAc                                                                                                                            <  2:123393/150‑1 (MQ=255)                   tGGCAAACGGCATTGTGTATATCTCCGAAGACCGTAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATACCGATGAACAGCAGGCtg                                                                                                                   >  2:76221/1‑150 (MQ=255)                   tGGCAAACGGCATTGTGTATATCTCCGAAGACCGTAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATACCGATGAACAGCAGGCtg                                                                                                                   >  2:76223/1‑150 (MQ=255)                      cAAACGGCATTGTGTATATCTCCGAAGACCGTAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTgtga                                                                                                                <  2:205130/150‑1 (MQ=255)                                      tatCTCCGAAGACCGTAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCt                                                                                                <  1:53458/150‑1 (MQ=255)                                                     tAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAAc                                                                                 >  2:235042/1‑150 (MQ=255)                                                        aCGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTcc                                                                              <  1:82105/150‑1 (MQ=255)                                                               ggTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATg                                                                       <  2:159992/150‑1 (MQ=255)                                                               ggTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATg                                                                       <  2:73986/150‑1 (MQ=255)                                                                        ttGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGAGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCa                                                              >  1:115764/1‑150 (MQ=255)                                                                                gTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCAATTGGTCt                                                      >  2:62432/1‑150 (MQ=255)                                                                                                         aCAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCaaaaa                             >  1:30548/1‑150 (MQ=255)                                                                                                                                    ggcggcAGTTTGAAGCATACCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGaca  <  1:76221/150‑1 (MQ=255)                                                                                                                                    ggcggcAGTTTGAAGCATACCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGaca  <  1:76223/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                    GTTCACCGCAGGATGGCCTGGCAAACGGCATTGTGTATATCTCCGAAGACCGTAAACGTGACGGTTTAGTGTTGGGCATGTCAGTAAAAGAGAACATGTCGCTGACAGCGCTGCGCTACTTCAGCCGCGCTGGCGGCAGTTTGAAGCATGCCGATGAACAGCAGGCTGTGAGTGATTTCATTCGTCTGTTTAATGTGAAAACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACA  >  NC_000913/3934733‑3935013

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