Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,459,646 C→T 15.8% L86L (CTG→TTG)  asnC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,459,6460CT15.8% 82.9 / 5.3 50L86L (CTG→TTG) asnCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (36/6);  new base T (6/2);  total (42/8)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTTATGACGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCACTACAGCATCTT  >  minE/2459581‑2459713
                                                                 |                                                                   
tCGGTTATGACGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTg                                                                   >  1:3304488/1‑68 (MQ=255)
tCGGTTATGACGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTg                                                                   >  1:1905122/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:1181569/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:60380/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:629883/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:706839/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:2497020/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:2229608/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:2076100/1‑68 (MQ=255)
         aCGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTg                                                          >  1:2027073/1‑68 (MQ=255)
                      tATCGGTATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTtgatga                                            >  1:198725/1‑69 (MQ=255)
                      tATCGGTATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTtgatga                                            >  1:1997397/1‑69 (MQ=255)
                      tATCGGTATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTtgatga                                            >  1:1854762/1‑69 (MQ=255)
                      tATCGGTATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTtgatga                                            >  1:1981380/1‑69 (MQ=255)
                      tATCGGTATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTtgatga                                            >  1:2401870/1‑69 (MQ=255)
                      tATCGGTATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTtgatga                                            >  1:2337637/1‑69 (MQ=255)
                                  atTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCaaa                                                              >  1:2871254/1‑39 (MQ=255)
                                  atTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCaaa                                                              >  1:1934342/1‑39 (MQ=255)
                                  atTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCaaa                                                              >  1:372098/1‑39 (MQ=255)
                                  atTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCaaa                                                              >  1:3224996/1‑39 (MQ=255)
                                   ttAAAGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGcc                               <  1:8280/69‑1 (MQ=255)
                                   ttAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGcc                               >  1:759237/1‑69 (MQ=255)
                                    tAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGcc                               <  1:448810/68‑1 (MQ=255)
                                     aaaGAGCGCCAAAGACTACCCTTCCGCGTTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGcc                               <  1:148239/67‑1 (MQ=255)
                                     aaaGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCt                              <  1:3191984/68‑1 (MQ=255)
                                      aaGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGcc                               <  1:1701627/66‑1 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTcgatga                                            >  1:664994/1‑42 (MQ=38)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:2636683/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:330796/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:3352426/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:3103187/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:1584758/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:23140/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:76421/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:783889/1‑59 (MQ=255)
                                                 aGACTATCCTTCCGCGCTGGCAAAGCTGGAAAGCCTCGATGAAGTCACTGAAGCCTAtt                           >  1:2046835/1‑59 (MQ=255)
                                                  gACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGc                <  1:3184726/69‑1 (MQ=255)
                                                  gACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGc                <  1:2966850/69‑1 (MQ=255)
                                                  gACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGc                <  1:1211219/69‑1 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTa                          >  1:1996582/1‑56 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCa              >  1:12087/1‑68 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCAc             >  1:3385814/1‑69 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCAc             >  1:2907633/1‑69 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCAc             >  1:2835759/1‑69 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCAc             >  1:2389292/1‑69 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCAc             >  1:2078541/1‑69 (MQ=255)
                                                     tACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCAc             >  1:920158/1‑69 (MQ=255)
                                                          ttCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCACTa           >  1:673106/1‑66 (MQ=255)
                                                          ttCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCACTa           >  1:2265599/1‑66 (MQ=255)
                                                          ttCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCACTa           >  1:880327/1‑66 (MQ=255)
                                                                gcTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCACTACAGCATCtt  <  1:1307628/69‑1 (MQ=255)
                                                                 |                                                                   
TCGGTTATGACGTAGGCTGCTTTATCGGCATTATATTAAAGAGCGCCAAAGACTACCCTTCCGCGCTGGCAAAGCTGGAAAGCCTTGATGAAGTCACTGAAGCCTATTACACAACCGGCCACTACAGCATCTT  >  minE/2459581‑2459713

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: