Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 258,220 C→T 7.4% P83S (CCG→TCG)  secF → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE258,2200CT 7.4% 212.3 / 4.0 67P83S (CCG→TCG) secFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base T (2/3);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGA  >  minE/258153‑258284
                                                                   |                                                                
taagCTCGAAAAACCGGCTGAAATTGACGTAATGTGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGa                                                               >  1:1592079/4‑71 (MQ=255)
tACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGa                                                               >  1:686854/1‑71 (MQ=255)
        aaaaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaa                                                       <  1:3595544/71‑1 (MQ=255)
         aaaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaa                                                       <  1:1791895/70‑1 (MQ=255)
          aaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaaaa                                                     <  1:3332424/71‑1 (MQ=255)
              cGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACttt                                                 <  1:2867943/71‑1 (MQ=255)
              cGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACttt                                                 <  1:500513/71‑1 (MQ=255)
              cGGCTGAAATTGACGTAATGCGGGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACttt                                                 <  1:2735763/71‑1 (MQ=255)
                  tGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTgg                                               >  1:2227712/1‑69 (MQ=255)
                         gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGTCGATGCTGCAAAACTTTGGTAGCAGCCa                                      >  1:391447/1‑71 (MQ=255)
                         gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGTCGATGCTGCAAAACTTTGGTAGCAGCCa                                      >  1:2215164/1‑71 (MQ=255)
                          aCGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCAt                                     <  1:1374942/71‑1 (MQ=255)
                               aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTag                                            >  1:2371277/1‑59 (MQ=255)
                               aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTag                                            >  1:1273833/1‑59 (MQ=255)
                               aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTag                                            >  1:3501673/1‑59 (MQ=255)
                                tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACtt                                                  <  1:1680632/52‑1 (MQ=255)
                                tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACtt                                                  <  1:3393813/52‑1 (MQ=255)
                                tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAcatc                               <  1:377504/71‑1 (MQ=255)
                                tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAcatc                               <  1:3237348/71‑1 (MQ=255)
                                tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAcatc                               <  1:3924109/71‑1 (MQ=255)
                                 gCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACtt                                                  <  1:932087/51‑1 (MQ=255)
                                  cGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGa                                                               >  1:3138500/1‑37 (MQ=255)
                                      aTGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTc                         >  1:1715939/1‑71 (MQ=255)
                                        gCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTcc                        >  1:1098539/1‑70 (MQ=255)
                                        gCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       >  1:999022/1‑71 (MQ=255)
                                         cATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:3670694/70‑1 (MQ=255)
                                         cATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:2597386/70‑1 (MQ=255)
                                         cATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:1908419/70‑1 (MQ=255)
                                         cATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:2156252/70‑1 (MQ=255)
                                         cATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCAGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:3160955/70‑1 (MQ=255)
                                          aTTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:3900390/69‑1 (MQ=255)
                                          aTTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:762545/69‑1 (MQ=255)
                                           ttGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTgg                                               <  1:830682/44‑1 (MQ=255)
                                           ttGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCg                       <  1:2278350/68‑1 (MQ=255)
                                           ttGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTAt                    >  1:2219667/1‑71 (MQ=255)
                                           ttGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTAt                    >  1:857660/1‑71 (MQ=255)
                                            tGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaa                                                       <  1:3349313/35‑1 (MQ=255)
                                            tGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaa                                                       <  1:3468944/35‑1 (MQ=255)
                                            tGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaa                                                       <  1:3830610/35‑1 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaaaa                                                     >  1:1578900/1‑36 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCaaaa                                                     >  1:380827/1‑36 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACtt                                                  >  1:2951083/1‑39 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTa                                             >  1:304150/1‑44 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGt                      >  1:4083332/1‑67 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGt                      >  1:1123697/1‑67 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGt                      >  1:1724235/1‑67 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTa                     >  1:525064/1‑68 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgc                  >  1:3120641/1‑71 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgc                  >  1:748147/1‑71 (MQ=255)
                                             gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgc                  >  1:3108921/1‑71 (MQ=255)
                                              cAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTa                     >  1:1632913/1‑67 (MQ=255)
                                                 aaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAcat                                >  1:3490105/1‑53 (MQ=255)
                                                 aaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAcat                                >  1:3014672/1‑53 (MQ=255)
                                                 aaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgccgcc              <  1:2367175/71‑1 (MQ=255)
                                                      cGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAc                                  <  1:2052457/46‑1 (MQ=255)
                                                       ggTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAc                                  <  1:568582/45‑1 (MQ=255)
                                                        gTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGa                                   >  1:3396300/1‑43 (MQ=255)
                                                           ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCGTGGTCCGTa                     >  1:788201/1‑54 (MQ=255)
                                                           ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATgg                           >  1:3928320/1‑48 (MQ=255)
                                                           ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTc                         >  1:2065419/1‑50 (MQ=255)
                                                           ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTa                     >  1:3557475/1‑54 (MQ=255)
                                                           ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTa                     >  1:2667062/1‑54 (MQ=255)
                                                           ttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGc    >  1:270146/1‑71 (MQ=255)
                                                            ttaaGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGc    >  1:843041/3‑70 (MQ=255)
                                                              aaGAGTCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGa  <  1:1465799/70‑1 (MQ=255)
                                                              aaGAGTCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGa  <  1:3479408/70‑1 (MQ=255)
                                                              aaGAGTCGATGCTGCAAAACTTTGCTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGa  <  1:1260753/70‑1 (MQ=255)
                                                                   |                                                                
TACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGA  >  minE/258153‑258284

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: