Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 893,993 T→A 36.5% S144S (TCA→TCT sapB ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE893,9930TA36.5% 37.8 / 42.4 52S144S (TCA→TCTsapBpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (29/4);  new base A (12/7);  total (41/11)
Fisher's exact test for biased strand distribution p-value = 7.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACA  >  minE/893929‑894057
                                                                |                                                                
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGaaa                                                              >  1:1056882/1‑69 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGaaa                                                              >  1:2639708/1‑69 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGaaa                                                              >  1:1765922/1‑69 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:2009602/69‑1 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:3373560/69‑1 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:3141029/69‑1 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:426177/69‑1 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:140004/69‑1 (MQ=255)
gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:1212219/69‑1 (MQ=255)
  cAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACg                                                            >  1:2435216/1‑69 (MQ=255)
   aGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCTTAGAGCAAATCGAAACGCCCAGaaa                                                              <  1:1512895/66‑1 (MQ=255)
        tcaatcaaCGGAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:2001492/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGGAACGCCCTGAAACGGGcaa                                                       >  1:2536368/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:3139734/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:1156423/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:2599396/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:2682367/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:1169620/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:397689/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:2944552/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:2990936/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:593209/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:1085407/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:1546313/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:147319/1‑68 (MQ=255)
        tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                       >  1:1263510/1‑68 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCt                                                >  1:1704022/1‑69 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:751245/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:2949504/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:3218815/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:3337582/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:2569411/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:430284/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:2048613/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:1998914/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:1566140/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:1563862/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                                            >  1:640432/1‑57 (MQ=255)
              aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCACAGAAACg                                                            >  1:3377857/1‑57 (MQ=255)
                cGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAg                                              >  1:294989/1‑69 (MQ=255)
                cGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAg                                              >  1:3292564/1‑69 (MQ=255)
                cGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAg                                              >  1:123739/1‑69 (MQ=255)
                            aaTCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGCAACCAACCTAGCGTGAGTGAACaa                                 <  1:1019194/69‑1 (MQ=255)
                            aaTCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGCAACCAACCTAGCGTGAGTGAACaa                                 <  1:519095/69‑1 (MQ=255)
                            aaTCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGCAACCAACCTAGCGTGAGTGAACaa                                 <  1:1548042/69‑1 (MQ=255)
                                          cGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAAc                     >  1:1000152/1‑68 (MQ=255)
                                                      cGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTc                        <  1:2284879/53‑1 (MQ=255)
                                                            gCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca  >  1:274569/1‑69 (MQ=255)
                                                            gCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca  >  1:2737526/1‑69 (MQ=255)
                                                            gCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca  >  1:2509614/1‑69 (MQ=255)
                                                            gCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca  >  1:1816892/1‑69 (MQ=255)
                                                            gCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca  >  1:1692744/1‑69 (MQ=255)
                                                                |                                                                
GCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACA  >  minE/893929‑894057

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: