Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	328,149	G→A	12.5%	G98G (GGC→GGT)	aes ←	acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	328,149	0	G	A	12.5%	184.2 / 15.6	64	G98G (GGC→GGT)	aes	acetyl esterase
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (4/4);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACAAAA  >  minE/328080‑328214                                                                      |                                                                  aaTACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTg                                                                   <  1:2527572/70‑1 (MQ=255) aaTACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTg                                                                   >  1:471185/1‑70 (MQ=255) aaTACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGc                                                                  >  1:1066599/1‑71 (MQ=255)    aCCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGcc                                                                 >  1:261462/1‑69 (MQ=255)    aCCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGcc                                                                 >  1:2752480/1‑69 (MQ=255)     ccAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTgac                                                                 >  1:561341/1‑66 (MQ=255)     ccAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGcc                                                                 >  1:2453525/1‑68 (MQ=255)     ccAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGcc                                                                 >  1:1020629/1‑68 (MQ=255)        aTCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           <  1:1028259/71‑1 (MQ=255)         tCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           >  1:3069826/1‑70 (MQ=255)         tCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           >  1:2757825/1‑70 (MQ=255)         tCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           >  1:2686388/1‑70 (MQ=255)         tCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           >  1:1130110/1‑70 (MQ=255)         tCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                                                          <  1:840752/71‑1 (MQ=255)          cACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                                                          >  1:2418826/1‑70 (MQ=255)           aCCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaa                                                        <  1:24801/71‑1 (MQ=255)             cGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                                                          <  1:12527/67‑1 (MQ=255)              gTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAAcc                                                     >  1:275173/1‑71 (MQ=255)               tACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAAcc                                                     >  1:1714688/1‑70 (MQ=255)                   ttGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           <  1:921648/60‑1 (MQ=255)                   ttGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                                                           <  1:2438630/60‑1 (MQ=255)                   ttGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                                                          <  1:2341231/61‑1 (MQ=255)                    tGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                                                          <  1:2597145/60‑1 (MQ=255)                    tGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTc                                                >  1:2403591/1‑70 (MQ=255)                     ggCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCa                                              >  1:1466514/1‑71 (MQ=255)                      gCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCa                                              <  1:1990508/70‑1 (MQ=255)                          tAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaa                                                         <  1:2452030/55‑1 (MQ=255)                          tAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaa                                                         <  1:1708705/55‑1 (MQ=255)                            gCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaa                                        >  1:601769/1‑70 (MQ=255)                                gCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCAt                                             >  1:477854/1‑61 (MQ=255)                                gCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaatga                                     <  1:1875424/69‑3 (MQ=255)                                gCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataaa                                    <  1:1419352/70‑1 (MQ=255)                                gCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataa                                     <  1:497124/69‑1 (MQ=255)                                 ccAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataaaat                                  <  1:2704043/71‑1 (MQ=255)                                 ccAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataaaa                                   <  1:2547678/70‑1 (MQ=255)                                 ccAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataa                                     <  1:87501/68‑1 (MQ=255)                                    gcagGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAAc                                                      >  1:2533646/1‑48 (MQ=255)                                    gcagGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGc                               >  1:915114/1‑71 (MQ=255)                                                tgCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTa                   <  1:190798/71‑1 (MQ=255)                                                tgCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTa                   <  1:2868066/71‑1 (MQ=255)                                                tgCGATCGTGGGTATCGAGATTACCGAGAATAAAACCGCCTCCATGCaaataaa                                    <  1:2997924/54‑1 (MQ=255)                                                tgCGATCGTGGGTATCGAGATTACCGAGAATAAAACCGCCTCCATGCaaataaa                                    <  1:1416874/54‑1 (MQ=255)                                                   gATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGTCTATCt                >  1:1911530/1‑71 (MQ=255)                                                   gATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCt                >  1:2367118/1‑71 (MQ=255)                                                          ggTATCGAGATTACCGAGAATAAAACCGCCTCCATGCAAATAAAATAGc                               <  1:2390281/49‑1 (MQ=255)                                                           gTATCGAGATTACCGAGAATAAAACCGCCTCCATGCaaata                                      <  1:1605376/41‑1 (MQ=255)                                                              tCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataa                                     <  1:1187221/39‑1 (MQ=255)                                                              tCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataa                                     <  1:3332249/39‑1 (MQ=255)                                                               cGAGATTACCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTgg                       >  1:3016965/1‑52 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAg                                >  1:2801189/1‑42 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAg                                >  1:2658459/1‑42 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATc                 >  1:1651585/1‑57 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    <  1:3033152/70‑1 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    <  1:2543994/70‑1 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    <  1:2660394/70‑1 (MQ=255)                                                                gagaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    >  1:2611310/1‑70 (MQ=255)                                                                gagaTTACCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTggg                      >  1:3122286/1‑52 (MQ=255)                                                                gagaTTACCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTggg                      >  1:1359885/1‑52 (MQ=255)                                                                gagaTTACCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTggg                      >  1:1877050/1‑52 (MQ=255)                                                                 agaTTGCCTAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    <  1:2956024/69‑1 (MQ=255)                                                                  gaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATGGCGTCGCTGGGCTATCTGGCTGCGGACaa    >  1:3306452/1‑68 (MQ=255)                                                                  gaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    >  1:3052106/1‑68 (MQ=255)                                                                  gaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    >  1:284975/1‑68 (MQ=255)                                                                  gaTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaa    >  1:459017/1‑68 (MQ=255)                                                                   aTTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaaaa  >  1:99711/1‑69 (MQ=255)                                                                    ttGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTa                   <  1:1647326/51‑1 (MQ=255)                                                                     tGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACaaa   >  1:2735395/1‑66 (MQ=255)                                                                      gCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGc                               >  1:2510659/1‑37 (MQ=255)                                                                      |                                                                  AATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACAAAA  >  minE/328080‑328214

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