Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 892,843 T→G 21.6% P210P (CCA→CCC sapC ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE892,8430TG21.6% 130.3 / 38.7 88P210P (CCA→CCCsapCpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (19/50);  new base G (15/4);  total (34/54)
Fisher's exact test for biased strand distribution p-value = 9.30e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTC  >  minE/892777‑892911
                                                                  |                                                                    
ggcgATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                    <  1:381973/69‑1 (MQ=255)
      atCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                    <  1:161701/63‑1 (MQ=255)
      atCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                    <  1:2248800/63‑1 (MQ=255)
      atCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                    <  1:2271852/63‑1 (MQ=255)
      atCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                    <  1:1131722/63‑1 (MQ=255)
       tCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGGGATGTTTgg                                                                    <  1:593389/62‑1 (MQ=255)
           gaATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCa                                                          <  1:2421709/68‑1 (MQ=255)
           gaATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCa                                                          <  1:612801/68‑1 (MQ=255)
                gCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGGATCACAGCAAAc                                                       >  1:1119272/1‑66 (MQ=255)
                gCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAAc                                                       >  1:2666440/1‑66 (MQ=255)
                gCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAAc                                                       >  1:2334949/1‑66 (MQ=255)
                gCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAAc                                                       >  1:3107049/1‑66 (MQ=255)
                   aTCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCagag                                                  <  1:65169/68‑1 (MQ=255)
                   aTCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCagag                                                  <  1:2745115/68‑1 (MQ=255)
                   aTCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCagag                                                  <  1:1837187/68‑1 (MQ=255)
                    tCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGaa                                                <  1:2371846/69‑1 (MQ=255)
                    tCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGaa                                                <  1:870467/69‑1 (MQ=255)
                    tCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGaa                                                <  1:397920/69‑1 (MQ=255)
                      gACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAata                                              <  1:918883/69‑1 (MQ=255)
                      gACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAata                                              <  1:3346013/69‑1 (MQ=255)
                      gACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAata                                              <  1:2719020/69‑1 (MQ=255)
                            gcgcGGGTGATCTCGGTGCCCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGc                                        >  1:1855753/1‑69 (MQ=255)
                            gcgcGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGc                                        >  1:2906796/1‑69 (MQ=255)
                            gcgcGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGc                                        >  1:500662/1‑69 (MQ=255)
                            gcgcGGGTGATCTAGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGc                                        >  1:160338/1‑69 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                    <  1:2761912/38‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:999822/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2528077/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2856702/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2929446/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3076568/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1041056/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:328100/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:631477/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1784199/69‑1 (MQ=255)
                               cgGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1066855/69‑1 (MQ=255)
                                gggTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3289337/68‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGTTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1086176/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1156054/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3376758/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2911285/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:38194/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3367913/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:576289/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3008233/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2732297/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:961690/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3081380/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:649632/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1980841/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:79285/65‑1 (MQ=255)
                                   tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2066516/65‑1 (MQ=255)
                                    gATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1468654/64‑1 (MQ=255)
                                    gATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:2955404/64‑1 (MQ=255)
                                    gATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1675991/64‑1 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1223775/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1209525/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:3342113/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:643202/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:3337878/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:981511/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1027800/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:218724/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1575077/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:2046891/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:2973052/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1663326/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1816381/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGCGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:870727/1‑68 (MQ=255)
                                         cGGTGACCAGTCCGGCGGTGATGTTGCGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1784185/1‑68 (MQ=255)
                                          ggTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:2143992/1‑67 (MQ=255)
                                          ggTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:3405688/1‑67 (MQ=255)
                                          ggTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCa                            >  1:1779976/1‑67 (MQ=255)
                                                cAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACg                    >  1:1775711/1‑69 (MQ=255)
                                                 aGCCCGGCTGTTATGTTTTGAATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:1340988/51‑1 (MQ=37)
                                                    ccGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcg                <  1:412224/69‑1 (MQ=255)
                                                    ccGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcg                <  1:2769438/69‑1 (MQ=255)
                                                    ccGGCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcg                <  1:403530/69‑1 (MQ=255)
                                                     cggcggTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcg                <  1:2215435/68‑1 (MQ=255)
                                                       gcggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcggc              >  1:3216355/1‑68 (MQ=255)
                                                       gcggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcggc              >  1:886378/1‑68 (MQ=255)
                                                       gcggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcggc              >  1:720264/1‑68 (MQ=255)
                                                           tGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:3253746/41‑1 (MQ=255)
                                                            gATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                     <  1:12725/40‑1 (MQ=255)
                                                            gATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATaa         >  1:1938368/1‑68 (MQ=255)
                                                            gATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATaa         >  1:1790523/1‑68 (MQ=255)
                                                            gATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATaa         >  1:1149835/1‑68 (MQ=255)
                                                            gATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATaa         >  1:1290370/1‑68 (MQ=255)
                                                                  tGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATgcgc                              >  1:1727678/1‑41 (MQ=255)
                                                                  tGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTActc  >  1:2232329/1‑69 (MQ=255)
                                                                  tGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTActc  >  1:1654845/1‑69 (MQ=255)
                                                                  tGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTActc  >  1:1060915/1‑69 (MQ=255)
                                                                  tGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTActc  >  1:3248945/1‑69 (MQ=255)
                                                                  |                                                                    
GGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTC  >  minE/892777‑892911

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: