Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,534,728	T→A	32.1%	V107V (GTT→GTA)	rfaG →	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,534,728	0	T	A	32.1%	54.3 / 36.1	53	V107V (GTT→GTA)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base T (16/20);  new base A (11/6);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTGAGCGA  >  minE/2534668‑2534794                                                             |                                                                   taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                                             >  1:3297807/1‑68 (MQ=255) taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                                             >  1:2617548/1‑68 (MQ=255) taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                                             >  1:2506974/1‑68 (MQ=255) taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                                             >  1:2395635/1‑68 (MQ=255)     aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                                         <  1:2979534/65‑1 (MQ=255)     aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                                         <  1:2885633/65‑1 (MQ=255)     aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                                         <  1:2703397/65‑1 (MQ=255)     aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                                         <  1:1460004/65‑1 (MQ=255)     aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                                         <  1:1605723/65‑1 (MQ=255)     aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                                         <  1:1910814/65‑1 (MQ=255)      aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                                       >  1:917974/1‑69 (MQ=255)      aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                                       >  1:881626/1‑69 (MQ=255)      aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                                       >  1:3335647/1‑69 (MQ=255)      aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                                       >  1:479151/1‑69 (MQ=255)       gATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAgg                                                      >  1:3159161/1‑69 (MQ=255)          gcctggcctgGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGtt                                                    >  1:3049928/1‑68 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:424304/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:23520/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:1971976/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:2435760/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:2899366/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:3380674/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:1081220/69‑1 (MQ=255)              ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGGTTTTta                                               <  1:91605/69‑1 (MQ=255)               gcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:1407054/68‑1 (MQ=255)                    gATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCg                                          >  1:1186593/1‑68 (MQ=255)                    gATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCg                                          >  1:255853/1‑68 (MQ=255)                      tGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGttt                                       <  1:2011566/69‑1 (MQ=255)                        tttATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTaa                                     <  1:493144/69‑1 (MQ=255)                        tttATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTaa                                     <  1:1366827/69‑1 (MQ=255)                        tttATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTaa                                     <  1:3414289/69‑1 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCgcc                                        >  1:2505501/1‑54 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:84776/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:816921/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:3286031/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:370144/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:2559078/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:2729634/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:2691857/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:1792444/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:1741560/1‑69 (MQ=255)                                  cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta                           >  1:1157282/1‑69 (MQ=255)                                  cGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGat                            <  1:2003054/68‑1 (MQ=255)                                      gATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGc                       <  1:1419719/69‑1 (MQ=255)                                            tgttACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                               <  1:1956919/39‑1 (MQ=255)                                            tgttACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGc                       >  1:696207/1‑63 (MQ=255)                                            tgttACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTAt                 >  1:847416/1‑69 (MQ=255)                                                       gaAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTga      <  1:2542783/69‑1 (MQ=255)                                                       gaAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTga      <  1:2076322/69‑1 (MQ=255)                                                       gaAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTga      <  1:1899643/69‑1 (MQ=255)                                                       gaAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTga      <  1:1830803/69‑1 (MQ=255)                                                           gTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTgagcga  >  1:2798889/1‑69 (MQ=255)                                                           gTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTgagcga  >  1:2386947/1‑69 (MQ=255)                                                             |                                                                   TAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTATGCCGCATTTGAGCGA  >  minE/2534668‑2534794

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