Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,516,230 T→C 59.2% intergenic (‑215/+44) cvpA ← / ← dedD membrane protein required for colicin V production/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,516,2300TC59.2% 10.7 / 27.3 44intergenic (‑215/+44)cvpA/dedDmembrane protein required for colicin V production/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/7);  new base C (13/13);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCGGGCAATGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGC  >  minE/1516165‑1516290
                                                                 |                                                            
gCATCGGGCAATGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTatg                                                           <  1:2872324/69‑3 (MQ=255)
gCATCGGGCAATGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAgg                                                           <  1:1923785/69‑1 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:2079937/1‑69 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:3069683/1‑69 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:2334916/1‑69 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:1542502/1‑69 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:165644/1‑69 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:2514363/1‑69 (MQ=255)
          aTGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGa                                                 >  1:2473092/1‑69 (MQ=255)
             cGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGAcgc                                              >  1:434621/1‑69 (MQ=255)
             cGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGAcgc                                              >  1:389857/1‑69 (MQ=255)
                      tGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGt                                      >  1:2731005/1‑68 (MQ=255)
                      tGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGt                                      >  1:2585682/1‑68 (MQ=255)
                          gcgTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                  <  1:1205899/68‑1 (MQ=255)
                          gcgTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                  <  1:1731852/68‑1 (MQ=255)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:374672/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:2873095/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:3326974/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:3168003/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:1885099/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:561681/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:1830008/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:612660/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:1651708/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:147323/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:1457175/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:1443978/1‑69 (MQ=21)
                                  gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                         >  1:1279956/1‑69 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:51339/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:663647/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:961474/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:3170951/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:1254729/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:3050464/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:3017230/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:257830/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:2397845/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:1872313/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:184224/69‑1 (MQ=21)
                                                      cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:1713975/69‑1 (MQ=21)
                                                      cTACACCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat     <  1:1341058/69‑1 (MQ=21)
                                                           ccTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCg        <  1:886418/61‑1 (MQ=255)
                                                           ccTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGc  <  1:2503991/67‑1 (MQ=255)
                                                           ccTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGc  <  1:1612176/67‑1 (MQ=255)
                                                                 |                                                            
GCATCGGGCAATGCGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGC  >  minE/1516165‑1516290

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: