Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,789,192	T→C	69.5%	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	69.5%	92.3 / 52.1	72	V110A (GTG→GCG)	ygaF	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/9);  new base C (30/20);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAAT  >  minE/1789126‑1789253                                                                   |                                                              gCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAgcggcg                                                           <  1:4197871/71‑1 (MQ=255) gCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAgcggcg                                                           <  1:4191463/71‑1 (MQ=255)      aTGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACgg                                                      >  1:1938126/1‑71 (MQ=255)      aTGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACgg                                                      >  1:3506582/1‑71 (MQ=255)          tGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATc                                                  >  1:2562983/1‑71 (MQ=255)            gTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATc                                                  >  1:5296072/1‑69 (MQ=255)              tGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAgc                                              >  1:5380786/1‑71 (MQ=255)              tGCAACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCg                                                           >  1:5491487/1‑58 (MQ=255)                ccACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAgcg                                             <  1:1586615/70‑1 (MQ=255)                 cACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:2972544/71‑1 (MQ=255)                 cACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:3663057/71‑1 (MQ=255)                 cACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:372861/71‑1 (MQ=255)                  aCGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGa                                          <  1:3594118/71‑1 (MQ=255)                  aCGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGa                                          <  1:2555292/71‑1 (MQ=255)                   cGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTAt                                                   >  1:2200941/1‑61 (MQ=255)                   cGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTAt                                                   >  1:4294875/1‑61 (MQ=255)                   cGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTAt                                                   >  1:5270501/1‑61 (MQ=255)                     tCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAgc                                              >  1:2913928/1‑64 (MQ=255)                     tCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:2114533/67‑1 (MQ=255)                     tCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:5324647/67‑1 (MQ=255)                     tCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:3127728/67‑1 (MQ=255)                     tCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGcgcg                                           <  1:3109438/67‑1 (MQ=255)                         aTCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAAt                                        <  1:2130554/66‑1 (MQ=255)                         aTCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAAt                                        <  1:2057484/66‑1 (MQ=255)                         aTCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAAt                                        <  1:3683753/66‑1 (MQ=255)                           ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTg                                                              >  1:1173125/1‑42 (MQ=255)                           ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACgg                                                      <  1:5683538/50‑1 (MQ=255)                           ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATgg                                      >  1:3529615/1‑66 (MQ=255)                           ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATgg                                      >  1:4137851/1‑66 (MQ=255)                           ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATgg                                      >  1:3923649/1‑66 (MQ=255)                           ctcGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATgg                                      >  1:4921675/1‑66 (MQ=255)                             cGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTaa                                 >  1:1765374/1‑69 (MQ=255)                             cGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTaa                                 >  1:2694128/1‑69 (MQ=255)                             cGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTaa                                 >  1:4631469/1‑69 (MQ=255)                              gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGaa                                                         >  1:5636977/1‑44 (MQ=255)                              gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCgg                                                             >  1:5625687/1‑40 (MQ=255)                              gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTaa                                 >  1:572425/1‑68 (MQ=255)                              gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTaa                                 >  1:1295349/1‑68 (MQ=255)                              gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTa                                  <  1:4648420/67‑1 (MQ=255)                              gAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGc                              >  1:2689268/1‑71 (MQ=255)                               aaaTGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGt                                                     <  1:2568672/47‑1 (MQ=255)                                aaTGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGcc                             <  1:4730679/70‑1 (MQ=255)                                   ggAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAgcgc                                            >  1:5325823/1‑52 (MQ=255)                                     aaCGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcgacg                         >  1:1254326/1‑69 (MQ=255)                                     aaCGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcgacg                         >  1:326815/1‑69 (MQ=255)                                      aCGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAAc                      >  1:885584/1‑71 (MQ=255)                                       cGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACg                                                       <  1:1369525/37‑1 (MQ=255)                                             cgcgcCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcg                <  1:2169292/70‑1 (MQ=255)                                               cgcCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagc             <  1:5436293/71‑1 (MQ=255)                                               cgcCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagc             <  1:5265993/71‑1 (MQ=255)                                                   ttATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCtt                                   <  1:3720088/45‑1 (MQ=255)                                                    tATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGaa        >  1:1143708/1‑71 (MQ=255)                                                     aTGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc       >  1:3204023/1‑71 (MQ=255)                                                     aTGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc       >  1:3963685/1‑71 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGaa                       >  1:4355376/1‑55 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGaa                       >  1:2356651/1‑55 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGaa                       >  1:4468889/1‑55 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGaa                       >  1:4606276/1‑55 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGaa                       >  1:4505272/1‑55 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagcgcga         <  1:4734565/69‑1 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagcgcga         <  1:5246421/69‑1 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagcg            >  1:4639595/1‑66 (MQ=255)                                                     aTGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc       >  1:1190001/1‑71 (MQ=255)                                                      tGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCAGACGAACTGCGCGAGCGCGAAcc      >  1:3209595/1‑71 (MQ=255)                                                        ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCg     <  1:3798290/70‑1 (MQ=255)                                                        ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCg     <  1:712803/70‑1 (MQ=255)                                                         gAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGaa   <  1:4654878/71‑1 (MQ=255)                                                          aaCGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAat  <  1:483757/71‑1 (MQ=255)                                                          aaCGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAat  <  1:1308795/71‑1 (MQ=255)                                                             gCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAcc      >  1:2687268/1‑64 (MQ=255)                                                             gCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAcc      >  1:1809422/1‑64 (MQ=255)                                                             gCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAcc      >  1:3167071/1‑64 (MQ=255)                                                                caGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGaa   >  1:2344669/1‑64 (MQ=255)                                                                  gTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcg                >  1:3196378/1‑49 (MQ=255)                                                                   |                                                              GCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAAT  >  minE/1789126‑1789253

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