Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,893,820 G→A 18.7% R132R (CGC→CGT mltA ← membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,893,8200GA18.7% 80.5 / 9.0 53R132R (CGC→CGTmltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base G (30/13);  new base A (3/7);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 3.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATA  >  minE/1893761‑1893879
                                                           |                                                           
gaCGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                    <  1:1075611/69‑1 (MQ=255)
gaCGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                    <  1:804421/69‑1 (MQ=255)
 aCGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                    <  1:1479982/68‑1 (MQ=255)
    gCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                    <  1:2650733/65‑1 (MQ=255)
    gCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                    <  1:797756/65‑1 (MQ=255)
     cAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                    <  1:1211785/64‑1 (MQ=255)
                      gcggcATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtg                               <  1:2360531/68‑1 (MQ=255)
                      gcggcATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtg                               <  1:1957127/68‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTTGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:82965/67‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:880593/67‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:2656369/67‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:2429606/67‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:641543/67‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:696677/67‑1 (MQ=255)
                        ggcATACGGTAAATCGGATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACAGGCgtgt                              <  1:2000558/67‑1 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtgt                              >  1:1674699/1‑57 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:640809/57‑1 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa                   >  1:2764166/1‑68 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa                   >  1:2534312/1‑68 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa                   >  1:1852046/1‑68 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa                   >  1:2025315/1‑68 (MQ=255)
                                  aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa                   >  1:2043328/1‑68 (MQ=255)
                                  aaaTAGGATACTGGAACTCCCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:1057793/57‑1 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGcc                                              >  1:554288/1‑40 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:1692877/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:2936152/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:1172879/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:2342482/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:2281342/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:2984449/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:2050350/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:40744/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            <  1:520124/42‑1 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTg                                            >  1:1729728/1‑42 (MQ=255)
                                   aaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtgt                              <  1:1037844/56‑1 (MQ=255)
                                       ggATACTGGAACTCGCCCTTACGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAGCCTGTAAACTGc             >  1:2734274/1‑69 (MQ=255)
                                       ggATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAGCCTGTAAACTGc             >  1:2575965/1‑69 (MQ=255)
                                       ggATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAGCCTGTAAACTGc             >  1:1785175/1‑69 (MQ=255)
                                           aCTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGt         >  1:2245052/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAAc                >  1:121105/1‑61 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:891671/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1006853/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1094074/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:819391/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1121640/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1494682/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1644289/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1776822/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:180677/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:1806833/1‑69 (MQ=255)
                                            cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:2457449/1‑69 (MQ=255)
                                            cTGGAACTAGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt        >  1:3346090/1‑69 (MQ=255)
                                                  cTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATa  <  1:3060947/69‑1 (MQ=255)
                                                           |                                                           
GACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATA  >  minE/1893761‑1893879

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: