Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,789,192	T→C	15.7%	V110A (GTG→GCG)	ygaF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,789,192	0	T	C	15.7%	120.0 / 15.9	51	V110A (GTG→GCG)	ygaF	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/24);  new base C (5/3);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATC  >  minE/1789122‑1789256                                                                       |                                                                 tGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGt                                                                  >  1:415693/1‑71 (MQ=255)   cGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTgg                                                                >  1:250291/1‑71 (MQ=255)   cGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTg                                                                 <  1:4057298/70‑1 (MQ=255)       aaGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGaa                                                            <  1:2881183/71‑1 (MQ=255)       aaGATGATGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGaa                                                            <  1:1058834/71‑1 (MQ=255)         gATGCTGGTTGCCACGTCCGATCTCTAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTg                                                                 <  1:3035189/64‑1 (MQ=255)          aTGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGaa                                                            <  1:3410932/68‑1 (MQ=255)           tGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTg                                                                 >  1:1625373/1‑62 (MQ=255)           tGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAAc                                                           >  1:3024596/1‑68 (MQ=255)                gTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCg                                                    <  1:2521660/70‑1 (MQ=255)                gTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCg                                                    <  1:2233710/70‑1 (MQ=255)                  tGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAgc                                                 >  1:827972/1‑71 (MQ=255)                     cACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACgg                                                         <  1:2185397/60‑1 (MQ=255)                     cACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACgg                                                         <  1:659947/60‑1 (MQ=255)                     cACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGagcg                                              <  1:1553738/71‑1 (MQ=255)                       cGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGaa                                            >  1:428108/1‑71 (MQ=255)                         tCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACg                                                          <  1:945102/55‑1 (MQ=255)                           cGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGcgcg                                              <  1:1322477/65‑1 (MQ=255)                             aTCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGcgcg                                              <  1:1660313/63‑1 (MQ=255)                             aTCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGaa                                            >  1:1295543/1‑65 (MQ=255)                               ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTaa                                    >  1:1791904/1‑71 (MQ=255)                                tcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGa                                                             <  1:2441980/45‑1 (MQ=255)                                tcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGcgcg                                              <  1:1957397/60‑1 (MQ=255)                                tcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAAc                                   >  1:1341832/1‑71 (MQ=255)                                tcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAAc                                   >  1:3429616/1‑71 (MQ=255)                                  gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGTCGAACGGTATCGAGCGCGa                                             >  1:2407809/1‑59 (MQ=255)                                  gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCtt                                      <  1:1677469/66‑1 (MQ=255)                                  gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTaa                                    >  1:3349582/1‑68 (MQ=255)                                  gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTaa                                    >  1:3451303/1‑68 (MQ=255)                                  gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGaa                                            >  1:893749/1‑60 (MQ=255)                                  gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAAc                                   <  1:3472976/69‑1 (MQ=255)                                   aaaTGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGcc                                <  1:3947357/71‑1 (MQ=255)                                      tGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTa                                                       <  1:3105354/45‑1 (MQ=255)                                       ggAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGaa                                            >  1:3820358/1‑55 (MQ=255)                                        gAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcgacg                            <  1:2349712/70‑1 (MQ=255)                                        gAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcgacg                            <  1:368606/70‑1 (MQ=255)                                         aaCGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATg                                          <  1:628262/55‑1 (MQ=255)                                            ggATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAg                                                  >  1:2018047/1‑44 (MQ=255)                                            ggATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAg                                                  >  1:3614573/1‑44 (MQ=255)                                                 cgcgcCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgc                    >  1:1702483/1‑69 (MQ=255)                                                 cgcgcCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcg                   <  1:3247920/70‑1 (MQ=255)                                                   cgcCTTATGGGAACGCACAGCGGCGAACGGTAACGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagc                >  1:2708986/1‑71 (MQ=255)                                                        tATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGaa           >  1:1588212/1‑71 (MQ=255)                                                         aTGGGAACGCACAGTGGCGAACGGTATCGAGCGCGaa                                            >  1:1032134/1‑37 (MQ=255)                                                         aTGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc          >  1:202237/1‑71 (MQ=255)                                                         aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc          >  1:2178513/1‑71 (MQ=255)                                                         aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc          >  1:3260558/1‑71 (MQ=255)                                                            ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCg        <  1:2318060/70‑1 (MQ=255)                                                            ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCg        <  1:1140257/70‑1 (MQ=255)                                                             gAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGaa      <  1:1769285/71‑1 (MQ=255)                                                                 gCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATc  <  1:117711/71‑1 (MQ=255)                                                                      gTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc          >  1:4031857/1‑58 (MQ=255)                                                                       |                                                                 TGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATC  >  minE/1789122‑1789256

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: ‑