Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,330,797	G→A	35.1%	G263G (GGC→GGT)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,330,797	0	G	A	35.1%	49.9 / 38.7	71	G263G (GGC→GGT)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (15/31);  new base A (20/5);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTGGATGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGC  >  minE/1330729‑1330859                                                                     |                                                               tcgtggatGATCTTCGCGCCGGGGTTTTTTTCGAGGAATGCTTCTGACAACAGGCCGACAATGtagtag                                                                <  1:167318/69‑1 (MQ=255) tcgtggatGATCTTCGCGCCGGGATTTTTTTCGAGGAATGGTTCTGCCAACAGGCCGACAATGtagtag                                                                <  1:3269845/69‑1 (MQ=255) tcgtggatGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGtagtag                                                                <  1:2370216/69‑1 (MQ=255) tcgtggatGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGtagtag                                                                <  1:1151498/69‑1 (MQ=255) tcgtggatGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGtagtag                                                                <  1:2658664/69‑1 (MQ=255)      gatgatCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTc                                                           >  1:1981383/1‑69 (MQ=255)      gatgatCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTc                                                           >  1:458243/1‑69 (MQ=255)      gatgatCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTc                                                           >  1:3309216/1‑69 (MQ=255)      gatgatCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTc                                                           >  1:3263189/1‑69 (MQ=255)        tgatCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCa                                                          >  1:452487/1‑68 (MQ=255)          atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                       <  1:318182/69‑1 (MQ=255)          atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                       <  1:3038899/69‑1 (MQ=255)          atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                       <  1:1870599/69‑1 (MQ=255)          atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                       <  1:3180679/69‑1 (MQ=255)                    cGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                       <  1:240957/59‑1 (MQ=255)                     gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                            <  1:2381860/69‑1 (MQ=255)                     gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                            <  1:3095917/69‑1 (MQ=255)                     gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                            <  1:545185/69‑1 (MQ=255)                     gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                            <  1:1956642/69‑1 (MQ=255)                        aTTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                            <  1:3264574/66‑1 (MQ=255)                                   ggAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGa                              <  1:142984/69‑1 (MQ=255)                                   ggAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGa                              <  1:597229/69‑1 (MQ=255)                                   ggAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGa                              <  1:2123316/69‑1 (MQ=255)                                   ggAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGa                              <  1:2559496/69‑1 (MQ=255)                                     aaTGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCGCTTTTTCGTCAAACAGGa                              <  1:2218604/67‑1 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCGCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:462958/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTc                                         >  1:2124637/1‑42 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCaaaa                >  1:1067448/1‑67 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:2860400/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:784869/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:3115419/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:2619519/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:3245679/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:1413144/1‑69 (MQ=255)                                                   cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc              >  1:1417242/1‑69 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:977129/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:890152/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:818698/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:723193/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:2625904/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:2647944/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:356813/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:147622/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:1566612/69‑1 (MQ=255)                                                     ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:2238521/69‑1 (MQ=255)                                                         gACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc            <  1:301194/65‑1 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCaa                                    >  1:201916/1‑40 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCaa                  >  1:1658574/1‑58 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1016839/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:2216921/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:2432324/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:2503757/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1168340/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:2167794/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1274904/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1857569/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1979716/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:2901752/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:3269569/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1856433/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:3160071/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1777826/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1712944/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:3057031/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:1878218/1‑68 (MQ=255)                                                          aCAATGTAGTAACCCTCGATAAACTGCCCTATTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc        >  1:6702/1‑68 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:1191739/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:1173719/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:1169831/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:3206164/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:782707/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:1718396/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:805761/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:807984/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:344326/69‑1 (MQ=255)                                                               gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:2960904/69‑1 (MQ=255)                                                                 agtagCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc  <  1:880948/67‑1 (MQ=255)                                                                     |                                                               TCGTGGATGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGC  >  minE/1330729‑1330859

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