Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 891,127 T→C 37.3% G154G (GGA→GGG sapF ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE891,1270TC37.3% 28.9 / 36.1 59G154G (GGA→GGGsapFpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (15/22);  new base C (5/17);  total (20/39)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAT  >  minE/891061‑891189
                                                                  |                                                              
ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGccc                                                              <  1:1078189/69‑1 (MQ=255)
ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGccc                                                              <  1:1829871/69‑1 (MQ=255)
ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGccc                                                              <  1:1444968/69‑1 (MQ=255)
ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGccc                                                              <  1:2655185/69‑1 (MQ=255)
ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGccc                                                              <  1:1358428/69‑1 (MQ=255)
ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGccc                                                              <  1:1082078/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACGCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1759805/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:691947/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:528612/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:3199875/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:202248/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1009602/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1120834/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:2624491/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:2393781/69‑1 (MQ=255)
          ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACACAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:2785350/69‑1 (MQ=255)
           taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:654236/69‑1 (MQ=255)
           taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:144028/69‑1 (MQ=255)
           taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:135431/69‑1 (MQ=255)
           taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:1338967/69‑1 (MQ=255)
           taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1887612/68‑1 (MQ=255)
            aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:2446006/1‑69 (MQ=255)
            aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:2513947/1‑69 (MQ=255)
            aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:84543/1‑69 (MQ=255)
            aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:1218344/1‑69 (MQ=255)
            aatGACTTTTGGGCCCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:3319989/1‑69 (MQ=255)
             atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 <  1:1908445/69‑1 (MQ=255)
             atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 <  1:534694/69‑1 (MQ=255)
             atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 >  1:1397968/1‑69 (MQ=255)
             atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:1169067/1‑68 (MQ=255)
                   tttGGGCGCAATATCAGCTCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:3244083/69‑1 (MQ=255)
                   tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:1308071/69‑1 (MQ=255)
                   tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:749808/69‑1 (MQ=255)
                   tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:1173222/69‑1 (MQ=255)
                   tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:3207814/69‑1 (MQ=255)
                   tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:1331132/69‑1 (MQ=255)
                   tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:1992712/69‑1 (MQ=255)
                        gcgcAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                       <  1:1492844/68‑1 (MQ=255)
                        gcgcAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                       <  1:2103818/68‑1 (MQ=255)
                           cAATATCGGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACt                                   <  1:578571/69‑1 (MQ=255)
                           cAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACt                                   <  1:1224908/69‑1 (MQ=255)
                           cAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACt                                   <  1:2457419/69‑1 (MQ=255)
                             atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:3092266/1‑69 (MQ=255)
                             atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:1633100/1‑69 (MQ=255)
                             atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:514984/1‑69 (MQ=255)
                             atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:3307869/1‑69 (MQ=255)
                                tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:1841128/1‑69 (MQ=255)
                                tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:711889/1‑69 (MQ=255)
                                tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:3093324/1‑69 (MQ=255)
                                tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:2151661/1‑69 (MQ=255)
                                    cgcgcgcgcCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 <  1:41012/46‑1 (MQ=255)
                                    cgcgcgcgcCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:2956091/52‑1 (MQ=255)
                                            ccagacccagacGCTGTTTTTGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTagca                  <  1:1665104/69‑1 (MQ=255)
                                            ccagacccagacGCTGTTTTTGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTagca                  <  1:76054/69‑1 (MQ=255)
                                            ccagacccagacGCTGTTTTTGCCCGGGGGCTAACATATGCGGATAGTAACTGACGTGATCCGGTagca                  <  1:1390962/69‑1 (MQ=255)
                                             cagacccagacGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagcag                 <  1:528389/69‑1 (MQ=255)
                                             cagacccagacGCTGTTTTTGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTagca                  <  1:1880784/68‑1 (MQ=255)
                                             cagacccagacGCTGTTTTTGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTagca                  <  1:3028326/68‑1 (MQ=255)
                                                      acGCTGTTTTTGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTagca                  <  1:795635/59‑1 (MQ=255)
                                                           gTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagcag                 <  1:421891/55‑1 (MQ=255)
                                                            tttttGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCa   >  1:466426/1‑68 (MQ=255)
                                                            tttttGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt  >  1:3260738/1‑69 (MQ=255)
                                                            tttttGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt  >  1:2201290/1‑69 (MQ=255)
                                                            tttttGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt  >  1:13751/1‑69 (MQ=255)
                                                            tttttGCCAGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt  >  1:2049549/1‑69 (MQ=255)
                                                                  |                                                              
TTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAT  >  minE/891061‑891189

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: