Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,512,108	C→T	39.7%	C141Y (TGC→TAC)	yicE ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,512,108	0	C	T	39.7%	45.4 / 77.3	68	C141Y (TGC→TAC)	yicE	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (15/26);  new base T (16/11);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 8.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AACGCCAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGTTTTC  >  minE/2512038‑2512176                                                                       |                                                                     aaCGCCAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGc                                                                      <  1:3246047/71‑1 (MQ=255)  aCGCCAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGc                                                                      >  1:1108289/1‑70 (MQ=255)  aCGCCAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGc                                                                      >  1:74167/1‑70 (MQ=255)    gCCAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAAc                                                                   >  1:2974357/1‑71 (MQ=255)      cAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTaa                                                                    >  1:1229338/1‑68 (MQ=255)        gAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGc                                                               >  1:363675/1‑71 (MQ=255)        gAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGc                                                               >  1:1771338/1‑71 (MQ=255)        gAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGc                                                               >  1:2003533/1‑71 (MQ=255)         aaaCCAGTGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGc                                                               <  1:2936444/70‑1 (MQ=255)         aaaCCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGcc                                                              <  1:1061977/71‑1 (MQ=255)         aaaCCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGc                                                               >  1:253129/1‑70 (MQ=255)         aaaCCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGc                                                               >  1:3695876/1‑70 (MQ=255)         aaaCCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGc                                                               >  1:1474897/1‑70 (MQ=255)           aCCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCa                                                             >  1:2386317/1‑70 (MQ=255)            ccAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCATCTTGCCAGc                                                           <  1:489650/71‑1 (MQ=255)            ccAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGc                                                           <  1:2026056/71‑1 (MQ=255)            ccAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGc                                                           <  1:657878/71‑1 (MQ=255)            ccAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGc                                                           <  1:3319696/71‑1 (MQ=255)             cAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCa                                                          <  1:2151899/71‑1 (MQ=255)               gcggcgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCAt                                                         >  1:1960688/1‑70 (MQ=255)                 ggcgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCaa                                                                    <  1:393122/57‑1 (MQ=255)                 ggcgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCaa                                                                    <  1:3404920/57‑1 (MQ=255)                 ggcgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCaa                                                                    <  1:3591727/57‑1 (MQ=255)                 ggcgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCaa                                                                    <  1:2438747/57‑1 (MQ=255)                 ggcgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCaa                                                      >  1:1491129/1‑71 (MQ=255)                   cgTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACt                                                                  >  1:416772/1‑57 (MQ=255)                    gTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGt                                                   <  1:3397496/71‑1 (MQ=255)                         aatGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCg                                               <  1:266794/70‑1 (MQ=255)                         aatGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCg                                               >  1:346551/1‑70 (MQ=255)                          atGCGGCGCGCCAGATTCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGa                                              <  1:2840941/70‑1 (MQ=255)                          atGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGa                                              <  1:1356626/70‑1 (MQ=255)                          atGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGa                                              <  1:2112655/70‑1 (MQ=255)                          atGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGa                                              <  1:2595439/70‑1 (MQ=255)                          atGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGa                                              <  1:3522615/70‑1 (MQ=255)                             cggcgCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACaa                                          <  1:2346587/71‑1 (MQ=255)                             cggcgCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCa                                                       <  1:2292695/58‑1 (MQ=255)                                   gcCAGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAgccgc                                    >  1:2699452/1‑71 (MQ=255)                                      aGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcat                                 >  1:258933/1‑71 (MQ=255)                                      aGATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcat                                 >  1:2496987/1‑71 (MQ=255)                                       gATGCAGAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcatc                                <  1:933810/71‑1 (MQ=255)                                       gATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAgccgc                                    <  1:1197027/67‑1 (MQ=255)                                           cAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATgg                            <  1:1771392/71‑1 (MQ=255)                                             gAACGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACaaa                                         <  1:3504434/56‑1 (MQ=255)                                             gAACGCGGGAGATCACCATCTCGGTGTAACTTGACAGCATCAACGTGCCGAACaaa                                         <  1:2289580/56‑1 (MQ=255)                                             gAACGCGGGAGATCACCATCTCGGGGTAACTTGCCAGCATCAACGTGCCGAACaaa                                         <  1:612277/56‑1 (MQ=255)                                              aaCGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcat                                 >  1:2996541/1‑63 (MQ=255)                                              aaCGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcat                                 >  1:3117430/1‑63 (MQ=255)                                              aaCGCGGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcat                                 >  1:3262882/1‑63 (MQ=255)                                               aCGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAgg                        <  1:2701290/71‑1 (MQ=255)                                                cgcgGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcatca                               >  1:2137474/1‑63 (MQ=255)                                                cgcgGGAGATCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAgg                        >  1:1767894/1‑70 (MQ=255)                                                    ggAGATCACCATCTCGGTGTAACTTGCCAGCATAAACGTg                                                  >  1:3531696/1‑40 (MQ=255)                                                         tCACCATCTCGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCcatc                                <  1:1899887/53‑1 (MQ=255)                                                             cATCTCGGTGCAACTTGCCAGCATCAACGTGCCGACCAAAgc                                       >  1:825331/1‑42 (MQ=255)                                                             cATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAg                         >  1:1205060/1‑56 (MQ=255)                                                             cATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGcaccac          >  1:2780329/1‑71 (MQ=255)                                                             cATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGcaccac          >  1:2567001/1‑71 (MQ=255)                                                             cATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGcaccac          <  1:1814471/71‑1 (MQ=255)                                                               tctcGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCg        >  1:3167434/1‑71 (MQ=255)                                                                ctcGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGcacc            <  1:3288235/66‑1 (MQ=255)                                                                ctcGGTGTAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGcacc            <  1:2692212/66‑1 (MQ=255)                                                                  cGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGtt     >  1:2529168/1‑71 (MQ=255)                                                                   ggTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCgg       <  1:212002/68‑1 (MQ=255)                                                                    gtgtAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGtttt   >  1:643980/1‑71 (MQ=255)                                                                    gTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTa                          <  1:2966047/48‑1 (MQ=255)                                                                    gTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTa                          <  1:3439206/48‑1 (MQ=255)                                                                    gTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGtttt   <  1:3393568/71‑1 (MQ=255)                                                                    gTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGtttt   <  1:3208060/71‑1 (MQ=255)                                                                     tGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTa                          >  1:141141/1‑47 (MQ=255)                                                                     tGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGtttt   >  1:531630/1‑70 (MQ=255)                                                                     tGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGTTTTc  <  1:3722379/71‑1 (MQ=255)                                                                      gCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGtttt   >  1:784385/1‑69 (MQ=255)                                                                       cAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAgg                        >  1:997877/1‑47 (MQ=255)                                                                       |                                                                     AACGCCAGAAACCAGCGGCGTAATAATGCGGCGCGCCAGATGCAGAACGCGGGAGATCACCATCTCGGTGCAACTTGCCAGCATCAACGTGCCGAACAAAGCCGCCATCATGGTAGGAACATCAGCACCACCGGTTTTC  >  minE/2512038‑2512176

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