Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 93,879 G→T 100% V287V (GTG→GTT 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
*minE93,8790GT75.7% 89.7 / 35.2 70V287V (GTG→GTTsecApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base G (12/5);  new base T (0/53);  total (12/58)
Fisher's exact test for biased strand distribution p-value = 5.82e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAAAATCTCGCCAGGTGAACCTGACCGAACGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGT  >  minE/93817‑93947
                                                              |                                                                    
gAAAAATCTCGCCAGGTGAACCTGACCGAACGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAagag                                                                <  1:287321/69‑1 (MQ=255)
gAAAAATCTCGCCAGGTGAACCTGACCGAACGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAagag                                                                <  1:2997656/69‑1 (MQ=255)
gAAAAATCTCGCCAGGTGAACCTGACCGAACGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAagag                                                                <  1:2018937/69‑1 (MQ=255)
gAAAAATCTCGCCAGGTGAACCTGACCGAACGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAagag                                                                <  1:2474609/69‑1 (MQ=255)
            cAGGTGAACCTGACCGATCGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1215017/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTTGTTAAAGAAGGcatcat                                                        <  1:2918473/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:524548/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:574536/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2149046/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:512175/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:451758/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:3355088/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2282170/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2288574/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:3189119/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:3117254/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2417649/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:3096329/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:3065536/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2619060/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2639122/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:2731262/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1981038/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1113964/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1238656/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:878543/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1468923/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1473518/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:816505/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1749396/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:712179/65‑1 (MQ=255)
            cAGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:185307/65‑1 (MQ=255)
             tggTGAACCTGACCGAACGTGGTATGGTTCTGATTGAAGAACTCCTGGTTAAAGAAGGcatctt                                                        <  1:683698/63‑3 (MQ=255)
             aGGTGAACCTGACCGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1719692/64‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:744762/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:780008/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:865838/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:739331/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:893668/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2900949/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2650475/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:918014/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2864640/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:184957/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2024265/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:1845062/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2190375/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2386611/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2491609/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2521436/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:1024700/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2650835/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2663221/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGGGAGtct                                        <  1:440775/67‑1 (MQ=255)
                          cGAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGAAGAAGGTGAGtct                                        <  1:1524197/67‑1 (MQ=255)
                           gAACGTGGTCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGCATCATGGATGAAGGTGAGtct                                        <  1:2314019/66‑1 (MQ=255)
                             aCGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTa                                   <  1:3435001/69‑1 (MQ=255)
                                   tCTGGTTCTGATTGAAGAACTGCTGGTTAAAGAAGGcatcat                                                        <  1:1315512/42‑1 (MQ=25)
                                                gaagaaCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTg                                     >  1:1541191/1‑48 (MQ=255)
                                                gaagaaCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTg                                     >  1:3294195/1‑48 (MQ=255)
                                                gaagaaCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTg                                     >  1:1813661/1‑48 (MQ=255)
                                                gaagaaCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTg                                     >  1:615093/1‑48 (MQ=255)
                                                     aCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTAct                                 >  1:2928701/1‑47 (MQ=255)
                                                     aCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGa            >  1:1697614/1‑68 (MQ=255)
                                                     aCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGa            >  1:2827515/1‑68 (MQ=255)
                                                     aCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGa            >  1:2127653/1‑68 (MQ=255)
                                                     aCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGa            >  1:2259276/1‑68 (MQ=255)
                                                     aCTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGa            >  1:3268904/1‑68 (MQ=255)
                                                      ctgctgGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATg          >  1:1565062/1‑69 (MQ=255)
                                                      ctgctgGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCAACATCATGCTGATGc         >  1:2235099/1‑69 (MQ=255)
                                                              gAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGt  <  1:3260530/69‑1 (MQ=255)
                                                              gAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGt  <  1:3206571/69‑1 (MQ=255)
                                                              gAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGt  <  1:1189179/69‑1 (MQ=255)
                                                              gAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGt  <  1:1052941/69‑1 (MQ=255)
                                                              |                                                                    
GAAAAATCTCGCCAGGTGAACCTGACCGAACGTGGTCTGGTGCTGATTGAAGAACTGCTGGTGAAAGAGGGCATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGT  >  minE/93817‑93947

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: