Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 891,136 C→T 37.2% L151L (TTG→TTA sapF ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE891,1360CT37.2% 34.2 / 33.8 59L151L (TTG→TTAsapFpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (15/22);  new base T (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 = 9.66e-01

ATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAT  >  minE/891071‑891189
                                                                 |                                                     
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACGCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1759805/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:3199875/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:528612/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:202248/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:691947/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:2624491/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:2393781/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1009602/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1120834/69‑1 (MQ=255)
ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACACAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:2785350/69‑1 (MQ=255)
 taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:1338967/69‑1 (MQ=255)
 taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:135431/69‑1 (MQ=255)
 taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:654236/69‑1 (MQ=255)
 taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACa                                                   <  1:144028/69‑1 (MQ=255)
 taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGCCCGGGTGCTAAc                                                    <  1:1887612/68‑1 (MQ=255)
  aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:2446006/1‑69 (MQ=255)
  aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:84543/1‑69 (MQ=255)
  aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:2513947/1‑69 (MQ=255)
  aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:1218344/1‑69 (MQ=255)
  aatGACTTTTGGGCCCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:3319989/1‑69 (MQ=255)
   atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 >  1:1397968/1‑69 (MQ=255)
   atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 <  1:1908445/69‑1 (MQ=255)
   atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACata                                                 <  1:534694/69‑1 (MQ=255)
   atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACat                                                  >  1:1169067/1‑68 (MQ=255)
         tttGGGCGCAATATCAGCTCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:3244083/69‑1 (MQ=255)
         tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:3207814/69‑1 (MQ=255)
         tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:1992712/69‑1 (MQ=255)
         tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGa                                           <  1:1331132/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)
              gcgcAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                       <  1:1492844/68‑1 (MQ=255)
              gcgcAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                       <  1:2103818/68‑1 (MQ=255)
                 cAATATCGGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACt                                   <  1:578571/69‑1 (MQ=255)
                 cAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACt                                   <  1:2457419/69‑1 (MQ=255)
                 cAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACt                                   <  1:1224908/69‑1 (MQ=255)
                   atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:514984/1‑69 (MQ=255)
                   atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:3092266/1‑69 (MQ=255)
                   atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:1633100/1‑69 (MQ=255)
                   atatCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGa                                 >  1:3307869/1‑69 (MQ=255)
                      tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:2151661/1‑69 (MQ=255)
                      tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:3093324/1‑69 (MQ=255)
                      tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:1841128/1‑69 (MQ=255)
                      tCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                              >  1:711889/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:13751/1‑69 (MQ=255)
                                                  tttttGCCCGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt  >  1:2201290/1‑69 (MQ=255)
                                                  tttttGCCAGGGTGCTAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt  >  1:2049549/1‑69 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:1097939/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:1200661/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:651541/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:1769934/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:3051009/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:2770912/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:2766718/44‑1 (MQ=255)
                                                                ccAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa             <  1:2719149/44‑1 (MQ=255)
                                                                 |                                                     
ATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAT  >  minE/891071‑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: