Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,507,099 T→C 28.5% L55L (TTG→CTG)  yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,507,0990TC28.5% 38.8 / 15.2 42L55L (TTG→CTG) yicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (14/16);  new base C (2/10);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 9.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCATTTTCAGGGGGCGC  >  minE/2507039‑2507162
                                                            |                                                               
tATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCttt                                                         >  1:2196990/1‑69 (MQ=255)
tATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCttt                                                         >  1:2158811/1‑69 (MQ=255)
tATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCttt                                                         >  1:1753723/1‑69 (MQ=255)
tATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCttt                                                         >  1:3411328/1‑69 (MQ=255)
tATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCtt                                                          >  1:731582/1‑68 (MQ=255)
tATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCtt                                                          >  1:1569550/1‑68 (MQ=255)
  tgctgcCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCttttt                                                       >  1:1129952/1‑69 (MQ=255)
      gcCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTtctc                                                    >  1:1449921/1‑68 (MQ=255)
                  gtgCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGt                                       <  1:662044/69‑1 (MQ=255)
                  gtgCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGt                                       <  1:1784217/69‑1 (MQ=255)
                  gtgCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGt                                       <  1:278883/69‑1 (MQ=255)
                  gtgCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGt                                       <  1:3339651/69‑1 (MQ=255)
                      gTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTg                                   >  1:1201877/1‑69 (MQ=255)
                      gTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTg                                   >  1:2218537/1‑69 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:2051798/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:2831462/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:693523/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:2933476/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:2120292/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:1391637/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:953795/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:1350787/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:962872/69‑1 (MQ=255)
                         aaCGTACCTGGCAGCTTGATACGCCTTTATTTACGCTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCg                                <  1:989485/69‑1 (MQ=255)
                               ccTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGtgtg                           >  1:718356/1‑68 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:2716636/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:837891/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:2931593/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:2727390/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:1051940/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:2330995/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:2208125/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:2088213/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:1496082/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:1352321/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCg                         <  1:1108293/68‑1 (MQ=255)
                                 tGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTAATGTCGGTGTGCgg                        <  1:72488/69‑1 (MQ=255)
                                         tGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCa                >  1:31803/1‑69 (MQ=255)
                                         tGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCa                >  1:3269137/1‑69 (MQ=255)
                                         tGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCa                >  1:1118200/1‑69 (MQ=255)
                                                        taCGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCATTTTCAGGGGgcgc  >  1:317382/1‑68 (MQ=255)
                                                        taCGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCATTTTCAGGGGgcgc  >  1:63320/1‑68 (MQ=255)
                                                            |                                                               
TATGCTGCCCCCCGTGATGTGCGTGAACGTACCTGGCAGCTTGATACGCCTTTATTTACGTTGCGCTTTTTCTCCCCACAGGAAGGTATTGTCGGTGTGCGGATTGAGCATTTTCAGGGGGCGC  >  minE/2507039‑2507162

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: