Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,388,213	A→G	35.1%	G67G (GGT→GGC)	yigL ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,388,213	0	A	G	35.1%	48.0 / 41.9	74	G67G (GGT→GGC)	yigL	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (22/4);  total (47/27)
Fisher's exact test for biased strand distribution p-value = 5.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02

CGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGCCCCACATC  >  minE/2388146‑2388272                                                                    |                                                            cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                                            <  1:971682/69‑1 (MQ=255) cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                                            <  1:2372789/69‑1 (MQ=255) cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                                            <  1:1735913/69‑1 (MQ=255) cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                                            <  1:1571577/69‑1 (MQ=255) cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                                            <  1:1455187/69‑1 (MQ=255)  gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAcc                                                           <  1:3008474/69‑1 (MQ=255)  gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAcc                                                           <  1:1436329/69‑1 (MQ=255)  gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                                            <  1:1896685/68‑1 (MQ=255)             cgcgATCCAGGTTATGATCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                                 >  1:2578235/1‑68 (MQ=255)             cgcgATCCAGGTTATGAGCAAATATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                                 >  1:1190348/1‑68 (MQ=255)             cgcgATCCAGGTTATGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                                 >  1:1134455/1‑68 (MQ=255)             cgcgATCCAGGTTATGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                                 >  1:1880586/1‑68 (MQ=255)             cgcgATCAAGGTTATGAGCAAAAATCAGGTTGCCATCAAGATCGTGAACGCGCGCGCCATTGGAGGTa                                                 >  1:276227/1‑68 (MQ=255)              gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                               >  1:41172/1‑69 (MQ=255)              gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                               <  1:1945759/69‑1 (MQ=255)              gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                               <  1:1874949/69‑1 (MQ=255)              gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                               >  1:2995231/1‑69 (MQ=255)              gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                               >  1:1017582/1‑69 (MQ=255)                   ccAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCAt                                            >  1:3142476/1‑67 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:2484541/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:2985614/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:1722473/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:865891/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:2906638/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:1355763/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:2640553/68‑1 (MQ=255)                      ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                                        <  1:3064321/68‑1 (MQ=255)                        ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                                      >  1:2604691/1‑68 (MQ=255)                        ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                                      >  1:1732665/1‑68 (MQ=255)                        ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                                      >  1:1484006/1‑68 (MQ=255)                         tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCa                                             >  1:234902/1‑60 (MQ=255)                         tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGAc                                     >  1:1152082/1‑68 (MQ=255)                         tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGAc                                     >  1:1853521/1‑68 (MQ=255)                          aTGAGCAAGAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:1341656/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAAt                                               >  1:1419630/1‑57 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATc                                              >  1:1795305/1‑58 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:2750752/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:2968155/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:3143561/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:521874/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:429727/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:2418287/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:1901317/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:1115152/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:137106/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:1738687/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:1761502/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:255217/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:1948842/1‑59 (MQ=255)                          aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                                             >  1:2253128/1‑59 (MQ=255)                          aTGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACtt                                   <  1:1985291/69‑1 (MQ=255)                             aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                                <  1:1119772/69‑1 (MQ=255)                             aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                                <  1:3063921/69‑1 (MQ=255)                             aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                                <  1:3309369/69‑1 (MQ=255)                             aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                                <  1:3297449/69‑1 (MQ=255)                                      cAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGAtt                       >  1:2107349/1‑69 (MQ=255)                                            aCCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATcgcg                 >  1:2492610/1‑69 (MQ=255)                                            aCCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATcgcg                 >  1:2745946/1‑69 (MQ=255)                                            aCCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATcgcg                 >  1:2056260/1‑69 (MQ=255)                                               aTCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGaa               <  1:447772/68‑1 (MQ=255)                                               aTCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGaa               <  1:492879/68‑1 (MQ=255)                                                tCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGaa               <  1:503228/67‑1 (MQ=255)                                                tCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAAtt             <  1:3107821/69‑1 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGATGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:3028982/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAATACTTAATCTCCAGATTATCGCGAATTTGcc         >  1:2331711/1‑67 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:3128404/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:2116520/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:39080/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:400053/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:1410660/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:664192/1‑68 (MQ=255)                                                      tCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGccc        >  1:2025770/1‑68 (MQ=255)                                                       cGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGCCCca      >  1:2653605/1‑69 (MQ=255)                                                           aaCGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGCCCCACATc  >  1:549764/1‑69 (MQ=255)                                                                    |                                                            CGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCGAATTTGCCCCACATC  >  minE/2388146‑2388272

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