Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,817,649 T→C 34.6% intergenic (‑138/+7) ygaD ← / ← mltB conserved hypothetical protein/membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,817,6490TC34.6% 52.1 / 42.2 75intergenic (‑138/+7)ygaD/mltBconserved hypothetical protein/membrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (18/31);  new base C (13/13);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGT  >  minE/1817582‑1817711
                                                                   |                                                              
tGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               <  1:2163415/69‑1 (MQ=255)
tGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               <  1:366070/69‑1 (MQ=255)
tGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               <  1:2045004/69‑1 (MQ=255)
tGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               <  1:1186550/69‑1 (MQ=255)
 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               <  1:1335878/68‑1 (MQ=255)
 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               >  1:768470/1‑68 (MQ=255)
 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTa                                                               <  1:1487174/68‑1 (MQ=255)
    aTTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAAt                                                            <  1:1446823/68‑1 (MQ=255)
    aTTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAAt                                                            <  1:2331331/68‑1 (MQ=255)
    aTTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAAt                                                            <  1:1931906/68‑1 (MQ=255)
                   ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACt                                                >  1:857203/1‑65 (MQ=255)
                   ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                             >  1:853197/1‑68 (MQ=255)
                   ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                             >  1:2371044/1‑68 (MQ=255)
                   ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                             >  1:750286/1‑68 (MQ=255)
                   ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                             >  1:2755764/1‑68 (MQ=255)
                   ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                             >  1:2761288/1‑68 (MQ=255)
                     cAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCgcgcc                                          >  1:510357/1‑69 (MQ=255)
                          aaaaGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCgcgccagcgc                                     <  1:177441/69‑1 (MQ=255)
                          aaaaGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCgcgccagcgc                                     <  1:3356617/69‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCATGGCtt                             <  1:1899802/67‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCtt                             <  1:1818623/67‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCtt                             <  1:2616369/67‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCtt                             <  1:1477686/67‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCtt                             <  1:1340276/67‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCtt                             <  1:1116261/67‑1 (MQ=255)
                                    tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTAGCGCCAGCGCCACGGCtt                             <  1:959621/67‑1 (MQ=255)
                                     aGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCgcgccagcgc                                     <  1:528069/58‑1 (MQ=255)
                                             ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                  <  1:390289/69‑1 (MQ=255)
                                             ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                  <  1:236306/69‑1 (MQ=255)
                                             ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                  <  1:1800881/69‑1 (MQ=255)
                                             ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                  <  1:214517/69‑1 (MQ=255)
                                             ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                  <  1:1298740/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:1067302/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:3263305/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:3062034/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:1404003/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:294974/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:2807587/69‑1 (MQ=255)
                                              gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:2580614/69‑1 (MQ=255)
                                               aaGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCgcgccagcgc                                     <  1:1908473/48‑1 (MQ=255)
                                               aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:2188396/68‑1 (MQ=255)
                                               aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:2657551/68‑1 (MQ=255)
                                               aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:2371927/68‑1 (MQ=255)
                                               aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:572902/68‑1 (MQ=255)
                                               aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:224114/68‑1 (MQ=255)
                                                aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:502465/1‑69 (MQ=255)
                                                aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:1702214/1‑69 (MQ=255)
                                                aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:2862365/1‑69 (MQ=255)
                                                aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:3211251/1‑69 (MQ=255)
                                                 gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:2609378/1‑68 (MQ=255)
                                                 gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:1510038/1‑68 (MQ=255)
                                                 gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg               >  1:1335522/1‑68 (MQ=255)
                                                 gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGa              >  1:2049657/1‑69 (MQ=255)
                                                 gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGa              >  1:387798/1‑69 (MQ=255)
                                                 gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGa              >  1:172528/1‑69 (MQ=255)
                                                    cACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                 <  1:2367887/63‑1 (MQ=255)
                                                       gCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGc          <  1:3037622/67‑1 (MQ=255)
                                                             ctctTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt  <  1:1556406/69‑1 (MQ=255)
                                                             ctctTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt  <  1:1450976/69‑1 (MQ=255)
                                                             ctctTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt  <  1:3313987/69‑1 (MQ=255)
                                                             ctctTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt  <  1:695141/69‑1 (MQ=255)
                                                             ctctTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt  <  1:70044/69‑1 (MQ=255)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTcc                         >  1:2631916/1‑44 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:661538/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:2154800/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:2137442/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:233472/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:2346581/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:721043/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:1030809/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:203348/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:1843265/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:2814402/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:2146073/1‑45 (MQ=25)
                                                               ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                        >  1:3298121/1‑45 (MQ=25)
                                                                   |                                                              
TGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGT  >  minE/1817582‑1817711

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: