Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	195,490	G→A	35.5%	I113I (ATC→ATT)	yaeF ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	195,490	0	G	A	35.5%	29.2 / 23.3	45	I113I (ATC→ATT)	yaeF	predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (20/9);  new base A (9/7);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

AGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTG  >  minE/195428‑195552                                                               |                                                               aGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTgg                                                          >  1:946600/1‑69 (MQ=255) aGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTgg                                                          >  1:3091197/1‑69 (MQ=255) aGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTg                                                           >  1:662097/1‑68 (MQ=255)            tCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGcc                                               <  1:1128871/69‑1 (MQ=255)            tCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGcc                                               <  1:2886216/69‑1 (MQ=255)            tCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGcc                                               <  1:1646384/69‑1 (MQ=255)             tgTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCt                                              >  1:2243725/2‑69 (MQ=255)             cGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCt                                              <  1:3072365/69‑1 (MQ=255)             cGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCt                                              <  1:2191653/69‑1 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCt                                                            >  1:962097/1‑52 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:2090345/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:2020087/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:1644026/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:2258966/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:1016572/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:254391/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:3122053/1‑69 (MQ=255)                aaGACGAAAAGCTTATCACGATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                                           >  1:2383825/1‑69 (MQ=255)                                 aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat                           <  1:1498399/68‑1 (MQ=255)                                 aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat                           <  1:1204502/68‑1 (MQ=255)                                 aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat                           <  1:2802817/68‑1 (MQ=255)                                 aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat                           <  1:1393675/68‑1 (MQ=255)                                    aTGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATc                       >  1:1345793/1‑69 (MQ=255)                                    aTGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATc                       >  1:3056864/1‑69 (MQ=255)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:1822739/56‑1 (MQ=38)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:3334691/56‑1 (MQ=38)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:3261085/56‑1 (MQ=38)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:3206085/56‑1 (MQ=38)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:301486/56‑1 (MQ=38)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:2692012/56‑1 (MQ=38)                                     tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc                                   <  1:1661434/56‑1 (MQ=38)                                       cTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCa                      >  1:523234/1‑67 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGt   >  1:1637102/1‑68 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGt   >  1:1149945/1‑68 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGt   >  1:755468/1‑68 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:2962330/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:2846788/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:2807482/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:1215075/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:141348/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:2139458/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:742225/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:903431/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:1096281/1‑69 (MQ=255)                                                         ggAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTg  >  1:1037878/1‑69 (MQ=255)                                                               |                                                               AGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCACCCAGAAAGATTGCAACGTG  >  minE/195428‑195552

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