Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,330,803	A→G	36.2%	I261I (ATT→ATC)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,330,803	0	A	G	36.2%	40.1 / 41.4	69	I261I (ATT→ATC)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (10/34);  new base G (20/5);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 8.60e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCAT  >  minE/1330736‑1330868                                                                    |                                                                  tgatCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCa                                                                   >  1:452487/1‑68 (MQ=255)   atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                                <  1:3038899/69‑1 (MQ=255)   atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                                <  1:1870599/69‑1 (MQ=255)   atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                                <  1:3180679/69‑1 (MQ=255)   atCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                                <  1:318182/69‑1 (MQ=255)             cGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATa                                                                <  1:240957/59‑1 (MQ=255)              gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                                     <  1:545185/69‑1 (MQ=255)              gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                                     <  1:3095917/69‑1 (MQ=255)              gggATTTTTTTCGAGGAACGCTTCTGCCAGCAAACCGACAATGTAGTAACCCTCGATAAACTGCCCttt                                                     <  1:2381860/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:2559496/69‑1 (MQ=255)                            ggAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGa                                       <  1:2123316/69‑1 (MQ=255)                            ggAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGa                                       <  1:597229/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:1413144/1‑69 (MQ=255)                                            cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc                       >  1:1417242/1‑69 (MQ=255)                                            cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc                       >  1:3115419/1‑69 (MQ=255)                                            cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc                       >  1:2860400/1‑69 (MQ=255)                                            cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc                       >  1:784869/1‑69 (MQ=255)                                            cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc                       >  1:2619519/1‑69 (MQ=255)                                            cAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATc                       >  1:3245679/1‑69 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:1566612/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:2647944/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:2625904/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:147622/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:2238521/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:356813/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:723193/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:818698/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:890152/69‑1 (MQ=255)                                              ggCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGc                     <  1:977129/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:3057031/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:3269569/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1274904/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1168340/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1016839/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:3160071/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1878218/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:2901752/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1856433/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:2503757/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1979716/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:2432324/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1712944/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:2216921/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1857569/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:2167794/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:1777826/1‑68 (MQ=255)                                                   aCAATGTAGTAACCCTCGATAAACTGCCCTATTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATc                 >  1:6702/1‑68 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:1169831/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:807984/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:805761/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:782707/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:1718396/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:1173719/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:344326/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:1191739/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:3206164/69‑1 (MQ=255)                                                        gtagtaGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:2960904/69‑1 (MQ=255)                                                          agtagCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:880948/67‑1 (MQ=255)                                                                  tCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCat  <  1:2073120/68‑1 (MQ=255)                                                                   cAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGc           <  1:1351754/58‑1 (MQ=255)                                                                    |                                                                  TGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCAT  >  minE/1330736‑1330868

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