Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,516,244 C→T 77.3% intergenic (‑229/+30) 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,2440CT77.3% 36.1 / 7.3 35intergenic (‑229/+30)cvpA/dedDmembrane protein required for colicin V production/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base T (13/14);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGCCCATTACCACGCCACTTAA  >  minE/1516178‑1516309
                                                                  |                                                                 
cGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGAcgc                                                                 >  1:389857/1‑69 (MQ=255)
cGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGAcgc                                                                 >  1:434621/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:3168003/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:2873095/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:3326974/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:374672/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:1885099/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:1830008/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:561681/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:1651708/1‑69 (MQ=21)
                     gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg                                            >  1:612660/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:3170951/69‑1 (MQ=21)
                                         cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat                        <  1:663647/69‑1 (MQ=21)
                                         cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat                        <  1:961474/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)
                                                     aGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat                        <  1:1814456/57‑1 (MQ=21)
                                                               aGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGCCCATTACCACGCCACTTaa  >  1:1279962/1‑69 (MQ=255)
                                                                  |                                                                 
CGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGCCCATTACCACGCCACTTAA  >  minE/1516178‑1516309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: