Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 95,017 T→C 8.3% L667L (TTG→CTG)  secA → preprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE95,0170TC 8.3% 343.1 / 10.7 165L667L (TTG→CTG) secApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base T (108/43);  new base C (4/10);  total (112/53)
Fisher's exact test for biased strand distribution p-value = 1.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 165 total aligned reads displayed.

GGAATATGATGACGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCATTG  >  minE/94956‑95083
                                                             |                                                                  
ggAATATGATGACGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGAt                                                               <  1:517849/67‑1 (MQ=255)
       gatgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCg                                                        >  1:1714837/1‑67 (MQ=255)
       gatgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2291133/1‑69 (MQ=255)
       gatgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2840572/1‑69 (MQ=255)
       gatgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2523682/1‑69 (MQ=255)
       gatgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1307087/1‑69 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1402329/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:3091363/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1933632/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:3249596/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1548513/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:266749/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2561963/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1419185/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1379924/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2478626/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2173761/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:1085846/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2387039/1‑68 (MQ=255)
        atgaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGAGTAACGAACTGTTGGATGTCAGCGAt                                                      >  1:2173663/1‑68 (MQ=255)
          gaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAtgtg                                                   <  1:1352258/69‑1 (MQ=255)
          gaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGAtgtg                                                   <  1:1288808/69‑1 (MQ=255)
          gaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGCTGGATGTCAGCGAtgtg                                                   <  1:2656835/69‑1 (MQ=255)
          gaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGCTGGATGTCAGCGAtgtg                                                   <  1:33798/69‑1 (MQ=255)
          gaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGCTGGATGTCAGCGAtgtg                                                   <  1:1893742/69‑1 (MQ=255)
          gaCGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGCTGGATGTCAGCGAtgtg                                                   <  1:917574/69‑1 (MQ=255)
               ggCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGCTGGATGTCAGCGATGTGAGCGa                                              >  1:3218013/1‑69 (MQ=255)
               ggCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGCTGGATGTCAGCGATGTGAGCGa                                              >  1:560996/1‑69 (MQ=255)
                         cAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAAc                                    >  1:3132873/1‑69 (MQ=255)
                                     aTTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGaa                        <  1:2876889/69‑1 (MQ=255)
                                     aTTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGaa                        <  1:2398167/69‑1 (MQ=255)
                                     aTTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGaa                        <  1:366659/69‑1 (MQ=255)
                                     aTTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGaa                        <  1:642498/69‑1 (MQ=255)
                                      tttACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGaa                        <  1:708378/68‑1 (MQ=255)
                                        tACTCCGAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:4691/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGTACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       <  1:457840/51‑1 (MQ=255)
                                        tACTCCCAGCGTAACGCACTGTTGGATGTCAGCGATGTGa                                                  >  1:1776655/1‑40 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       <  1:2588183/51‑1 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:236151/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:855489/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       <  1:833053/51‑1 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       <  1:795471/51‑1 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:3221540/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:968083/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:2958286/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:3002872/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       <  1:2273428/51‑1 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:1252371/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:1160992/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:1962659/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCAtt                                       >  1:1107182/1‑51 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGc                                 >  1:1606509/1‑57 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGa                      >  1:2696430/1‑68 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGAt                     >  1:3315418/1‑69 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGAt                     >  1:566676/1‑69 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGAt                     >  1:1118014/1‑69 (MQ=255)
                                        tACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGAt                     >  1:2555184/1‑69 (MQ=255)
                                            cccAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGAtgtg                  <  1:1424058/68‑1 (MQ=255)
                                                gCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCaa             <  1:1523000/69‑1 (MQ=255)
                                                   tAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGc          <  1:740580/69‑1 (MQ=255)
                                                   tAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGc          <  1:2951124/69‑1 (MQ=255)
                                                    aaCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         <  1:2164326/69‑1 (MQ=255)
                                                    aaCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         <  1:1768827/69‑1 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAACCCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2544120/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTaa                                     >  1:2709286/1‑40 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAgag         >  1:3085176/1‑66 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGc          >  1:2115112/1‑67 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGc          >  1:1049372/1‑67 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1004293/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1051386/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2344599/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:3435620/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:34890/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:151648/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:417320/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1454885/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:103939/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:632238/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1936462/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2339902/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2445322/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2075602/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2651490/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:19868/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1939783/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2845752/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:3286142/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:2937391/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1877088/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1871825/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:3031224/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1808957/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCg         >  1:1631623/1‑68 (MQ=255)
                                                     aCGAACTGTTGGATGTCAGCCATGTGAGCGAAACCATTAACAg                                  >  1:804171/1‑43 (MQ=255)
                                                       gAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGAc       <  1:358663/68‑1 (MQ=255)
                                                          cTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCAtt   <  1:3012897/69‑1 (MQ=255)
                                                           tGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCATTg  <  1:1463445/69‑1 (MQ=255)
                                                           tGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCATTg  <  1:275475/69‑1 (MQ=255)
                                                           tGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCATTg  <  1:2036025/69‑1 (MQ=255)
                                                           tGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCATTg  <  1:1565775/69‑1 (MQ=255)
                                                             |                                                                  
GGAATATGATGACGTGGCTAACGATCAGCGTCGCGCCATTTACTCCCAGCGTAACGAACTGTTGGATGTCAGCGATGTGAGCGAAACCATTAACAGCATTCGTGAAGATGTGTTCAAAGCGACCATTG  >  minE/94956‑95083

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: