Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,515,612 G→A 10.9% S135L (TCA→TTA)  cvpA ← membrane protein required for colicin V production

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,515,6120GA10.9% 188.6 / 11.9 64S135L (TCA→TTA) cvpAmembrane protein required for colicin V production
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (4/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTTGACGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAAC  >  minE/1515543‑1515679
                                                                     |                                                                   
aaaCTTGACGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGa                                                                    <  1:1313561/71‑1 (MQ=255)
aaaCTTGACGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGa                                                                    <  1:2121717/71‑1 (MQ=255)
    ttGACGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTg                                                                >  1:4488989/1‑71 (MQ=255)
    ttGACGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTg                                                                >  1:3181876/1‑71 (MQ=255)
       aCGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTc                                                             <  1:370359/71‑1 (MQ=255)
       aCGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTc                                                             <  1:3199788/71‑1 (MQ=255)
        cGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTc                                                             <  1:1876323/70‑1 (MQ=255)
               tGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTc                                                     >  1:3931873/1‑71 (MQ=255)
                 cAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGc                                                   >  1:4004905/1‑71 (MQ=255)
                 cAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGc                                                   >  1:1739177/1‑71 (MQ=255)
                   gATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCtt                                                 <  1:4338920/71‑1 (MQ=255)
                   gATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCtt                                                 <  1:2404010/71‑1 (MQ=255)
                    aTAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCttt                                                >  1:1979979/1‑71 (MQ=255)
                     tAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCttt                                                <  1:457021/70‑1 (MQ=255)
                      aaTCAAAAAACCATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:2805215/1‑71 (MQ=255)
                      aaTCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCtttt                                               >  1:4619063/1‑70 (MQ=255)
                      aaTCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:4318263/1‑71 (MQ=255)
                      aaTCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:4139084/1‑71 (MQ=255)
                      aaTCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:2597708/1‑71 (MQ=255)
                       aTCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              <  1:4569285/70‑1 (MQ=255)
                       aTCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              <  1:2291759/70‑1 (MQ=255)
                        tCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGa                                                                    <  1:2026969/47‑1 (MQ=255)
                        tCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGa                                                                    <  1:2268818/47‑1 (MQ=255)
                        tCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGa                                                                    <  1:3874826/47‑1 (MQ=255)
                        tCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGa                                                                    <  1:3007757/47‑1 (MQ=255)
                        tCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGt                                                         >  1:2566718/1‑58 (MQ=255)
                          aaaaaaaCATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTg                                                                >  1:1520839/1‑49 (MQ=255)
                          aaaaaaaCATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTg                                                                >  1:4736554/1‑49 (MQ=255)
                                aCATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCtttt                                               >  1:2006008/1‑60 (MQ=255)
                                aCATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCtttt                                               >  1:2773489/1‑60 (MQ=255)
                                     tgatgaTAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGt                                                         >  1:3486812/1‑45 (MQ=255)
                                        tgatAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACAcccc                                      <  1:2022925/61‑1 (MQ=255)
                                          atAAAACTGAATTGCGGGATCAGCTGTGATTTGCt                                                              >  1:1756991/1‑35 (MQ=255)
                                          atAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTc                          <  1:2965375/71‑1 (MQ=255)
                                          atAAAACTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTCtt                                                      <  1:2308435/43‑1 (MQ=255)
                                          atAAAACTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTCtt                                                      <  1:3331424/43‑1 (MQ=255)
                                          atAAAACTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTAtt                                                      <  1:3261911/43‑1 (MQ=38)
                                           tAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTc                          <  1:104820/70‑1 (MQ=255)
                                           tAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTc                          <  1:1334322/70‑1 (MQ=255)
                                            aaaaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGt                                                         >  1:4035569/1‑38 (MQ=255)
                                            aaaaCTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTCTTCg                                                    >  1:1255848/1‑43 (MQ=255)
                                            aaaaCTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTCTTCg                                                    >  1:384815/1‑43 (MQ=255)
                                            aaaaCTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTCTTCg                                                    >  1:1250386/1‑43 (MQ=255)
                                            aaaaCTGAATTGCGGGATCAGCTGTAATTTGCTCCAGTCTTCg                                                    >  1:4362712/1‑43 (MQ=255)
                                             aaaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCgag                       <  1:1598720/71‑1 (MQ=255)
                                              aaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGt                                                         >  1:1212509/1‑36 (MQ=255)
                                              aaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCg                                             >  1:4054237/1‑48 (MQ=255)
                                              aaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCgag                       >  1:3473187/1‑70 (MQ=255)
                                              aaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCgag                       >  1:3317677/1‑70 (MQ=255)
                                              aaCTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCgag                       >  1:3833195/1‑70 (MQ=255)
                                                 tGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:4079017/1‑44 (MQ=255)
                                                 tGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:1897216/1‑44 (MQ=255)
                                                 tGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:1809296/1‑44 (MQ=255)
                                                 tGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTc                                              >  1:1902773/1‑44 (MQ=255)
                                                       gCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAagag              <  1:3594685/70‑1 (MQ=255)
                                                       gCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAagag              <  1:327547/70‑1 (MQ=255)
                                                       gCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAagag              <  1:4117062/70‑1 (MQ=255)
                                                       gCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAagag              <  1:2448822/70‑1 (MQ=255)
                                                          ggATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAATGAAGAGAAt           <  1:4328983/70‑1 (MQ=255)
                                                            aTCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGgc        <  1:733929/71‑1 (MQ=255)
                                                              cAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGagaaagaa                 <  1:2328901/60‑1 (MQ=255)
                                                               aGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGagaaagaa                 <  1:1121443/59‑1 (MQ=255)
                                                                 cTGTGATTTGGTCCTGTCTTCGGTTTTCGACAGCCCGGTAATGGAGTCGagaaag                   <  1:3460499/55‑1 (MQ=255)
                                                                 cTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGagaaag                   <  1:4033889/55‑1 (MQ=255)
                                                                  tgtgATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCg                         <  1:3780605/48‑1 (MQ=255)
                                                                  tgtgATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCg                         <  1:3185365/48‑1 (MQ=255)
                                                                   gtgATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGcaac  >  1:1299505/1‑70 (MQ=255)
                                                                     |                                                                   
AAACTTGACGAGCTTTGCAGATAATCAAAAAAACATCTGATGATAAAACTGAATTGCGGGATCAGCTGTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAAC  >  minE/1515543‑1515679

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: