Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,420,124 A→G 72.4% D131D (GAT→GAC ccmG ← periplasmic thioredoxin of cytochrome c‑type biogenesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,420,1240AG72.4% 57.6 / 19.3 54D131D (GAT→GACccmGperiplasmic thioredoxin of cytochrome c‑type biogenesis
Reads supporting (aligned to +/- strand):  ref base A (5/10);  new base G (18/21);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

TAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGATTGCCTTCTGGCG  >  minE/1420072‑1420192
                                                    |                                                                    
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1027989/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:913070/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:871622/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:757087/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:339265/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:327924/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:2631420/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:2364156/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:2323622/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1826908/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1565670/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1536397/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1223780/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1078358/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1220519/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1132193/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCGTCGCCATCAAACAgg                                                       >  1:1324782/1‑68 (MQ=255)
tAAGGAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTACATGCCGTCGCCATCAAACAgg                                                       >  1:1398522/1‑67 (MQ=255)
     aaCGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAgc                                                 <  1:574282/69‑1 (MQ=255)
     aaCGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAgc                                                 <  1:1792931/69‑1 (MQ=255)
             tGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAg                                          <  1:1068302/68‑1 (MQ=255)
             tGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAg                                          <  1:1782038/68‑1 (MQ=255)
             tGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAg                                          <  1:482207/68‑1 (MQ=255)
                         gACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGctc                             >  1:1786987/1‑69 (MQ=255)
                         gACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGctc                             >  1:1630574/1‑69 (MQ=255)
                         gACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGctc                             >  1:3025377/1‑69 (MQ=255)
                         gACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGctc                             >  1:2576921/1‑69 (MQ=255)
                         gACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGctc                             >  1:13447/1‑69 (MQ=255)
                              cGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTca                        <  1:1155528/69‑1 (MQ=255)
                              cGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTca                        <  1:137312/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:3258438/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:1009585/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:873073/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:134801/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:760476/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:1346070/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:580593/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:1670311/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:56104/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:1746504/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:1098613/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:3307903/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:1948450/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:3204997/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:3185789/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:2038689/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:2974653/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:2865868/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:2653350/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:2587078/69‑1 (MQ=255)
                                       gCCCTAACATGCCGTCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGa               <  1:2489891/69‑1 (MQ=255)
                                       gCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTg                <  1:2302505/68‑1 (MQ=255)
                                       gCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTg                <  1:2978906/68‑1 (MQ=255)
                                       gCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTg                <  1:2144943/68‑1 (MQ=255)
                                                    atcgccatcAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGATTGCCTTCTGGCg  <  1:2829911/69‑1 (MQ=255)
                                                    atcgccatcAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGATTGCCTTCTGGCg  <  1:2085309/69‑1 (MQ=255)
                                                    |                                                                    
TAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCCTAACATGCCATCGCCATCAAACAGGCTTAGCGCGTAAGGATTGCCCAGCTCTTTCAGCCAGCTGATTGCCTTCTGGCG  >  minE/1420072‑1420192

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: