Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,817,662 T→C 35.5% V360V (GTA→GTG mltB ← membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,817,6620TC35.5% 48.7 / 44.4 76V360V (GTA→GTGmltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base C (13/14);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGT  >  minE/1817601‑1817729
                                                             |                                                                   
ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACt                                                                  >  1:857203/1‑65 (MQ=255)
ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                                               >  1:2371044/1‑68 (MQ=255)
ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                                               >  1:2755764/1‑68 (MQ=255)
ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                                               >  1:853197/1‑68 (MQ=255)
ttCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTcgc                                                               >  1:750286/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:236306/69‑1 (MQ=255)
                          ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                                    <  1:214517/69‑1 (MQ=255)
                          ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                                    <  1:1800881/69‑1 (MQ=255)
                          ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                                    <  1:390289/69‑1 (MQ=255)
                          ggAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                                    <  1:1298740/69‑1 (MQ=255)
                           gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:1067302/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)
                           gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:3062034/69‑1 (MQ=255)
                           gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:1404003/69‑1 (MQ=255)
                           gAAGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:3263305/69‑1 (MQ=255)
                            aaGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCgcgccagcgc                                                       <  1:1908473/48‑1 (MQ=255)
                            aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:2188396/68‑1 (MQ=255)
                            aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:224114/68‑1 (MQ=255)
                            aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:2657551/68‑1 (MQ=255)
                            aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:572902/68‑1 (MQ=255)
                            aaGCTCACGCCCCCCTCTTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:2371927/68‑1 (MQ=255)
                             aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:3211251/1‑69 (MQ=255)
                             aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:2862365/1‑69 (MQ=255)
                             aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:1702214/1‑69 (MQ=255)
                             aGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:502465/1‑69 (MQ=255)
                              gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:1510038/1‑68 (MQ=255)
                              gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:2609378/1‑68 (MQ=255)
                              gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                 >  1:1335522/1‑68 (MQ=255)
                              gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGa                                >  1:387798/1‑69 (MQ=255)
                              gCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGa                                >  1:2049657/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:2814402/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:661538/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:3298121/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:721043/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:1030809/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:2346581/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:2137442/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:203348/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:2146073/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:2154800/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:233472/1‑45 (MQ=25)
                                            ctTGCAAATAGTTACTGCACGCGCGCCAGCGCCACGGCTTGTCCt                                          >  1:1843265/1‑45 (MQ=25)
                                                     aGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtg         >  1:225082/1‑69 (MQ=255)
                                                     aGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtg         >  1:2739348/1‑69 (MQ=255)
                                                      gttACTGCACGCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                   <  1:319774/40‑1 (MQ=37)
                                                        tACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg       <  1:2505156/68‑1 (MQ=255)
                                                        tACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg       <  1:1477113/68‑1 (MQ=255)
                                                        tACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg       <  1:2713450/68‑1 (MQ=255)
                                                        tACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg       <  1:2809156/68‑1 (MQ=255)
                                                        tACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg       <  1:1529/68‑1 (MQ=255)
                                                        tACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg       <  1:2939767/68‑1 (MQ=255)
                                                            gTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGt  <  1:725250/69‑1 (MQ=255)
                                                            gTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGt  <  1:1375544/69‑1 (MQ=255)
                                                             |                                                                   
TTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGT  >  minE/1817601‑1817729

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: