Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 893,980 G→A 56.2% L149L (CTG→TTG)  sapB ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE893,9800GA56.2% 15.0 / 56.4 73L149L (CTG→TTG) sapBpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (28/4);  new base A (34/7);  total (62/11)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAAC  >  minE/893912‑894038
                                                                    |                                                          
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:911870/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:3091752/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:3105682/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:2932527/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:2897472/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:283688/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:2720287/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:3289028/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:2254880/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:3433091/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:402429/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:404972/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:1479098/69‑1 (MQ=255)
ccATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:141325/69‑1 (MQ=255)
         gTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAg                                                            <  1:818282/60‑1 (MQ=255)
          tCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGc                                                  >  1:1319028/1‑69 (MQ=255)
              aaaGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGc                                                  >  1:2198649/1‑65 (MQ=255)
              aaaGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGc                                                  >  1:746709/1‑65 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGaaa                                           >  1:1056882/1‑69 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGaaa                                           >  1:2639708/1‑69 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGaaa                                           >  1:1765922/1‑69 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                           <  1:140004/69‑1 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                           <  1:3373560/69‑1 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                           <  1:426177/69‑1 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                           <  1:1212219/69‑1 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                           <  1:3141029/69‑1 (MQ=255)
                 gCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGaaa                                           <  1:2009602/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:1169620/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:1156423/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:3139734/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:2682367/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:1085407/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:2944552/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:2990936/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:147319/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:2599396/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:1263510/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:397689/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:593209/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                    >  1:1546313/1‑68 (MQ=255)
                         tcaatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGAAGATCGAAACGccc                                                >  1:3246481/1‑56 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCg                                                        >  1:2402836/1‑46 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCg                                                        >  1:620898/1‑46 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:3340581/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:368749/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:390034/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:723497/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:568268/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:60618/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:3111664/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:2730567/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:3102611/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:1952717/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:2018467/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:158146/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:234049/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:2617349/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:1351899/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:2810396/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:2919360/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:1542813/1‑47 (MQ=255)
                           aatcaaCGCAAAACCGCTAATCGGTTTCACTTCGTAGAGCAAATCGa                                                       >  1:336841/1‑47 (MQ=255)
                             tcaaCGCAACACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGc                                                  <  1:126961/50‑1 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCt                             >  1:1704022/1‑69 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATc                                                         >  1:1941974/1‑41 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:1563862/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:1566140/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:430284/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:640432/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:1998914/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:2048613/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:2569411/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:3337582/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:751245/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:3218815/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCCCAGAAACg                                         >  1:2949504/1‑57 (MQ=255)
                               aaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAAATCGAAACGCACAGAAACg                                         >  1:3377857/1‑57 (MQ=255)
                                 cGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAg                           >  1:123739/1‑69 (MQ=255)
                                 cGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAg                           >  1:294989/1‑69 (MQ=255)
                                 cGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAg                           >  1:3292564/1‑69 (MQ=255)
                                             aaTCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGCAACCAACCTAGCGTGAGTGAACaa              <  1:1548042/69‑1 (MQ=255)
                                             aaTCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGCAACCAACCTAGCGTGAGTGAACaa              <  1:519095/69‑1 (MQ=255)
                                             aaTCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGCAACCAACCTAGCGTGAGTGAACaa              <  1:1019194/69‑1 (MQ=255)
                                                           cGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAAc  >  1:1000152/1‑68 (MQ=255)
                                                                    |                                                          
CCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAAC  >  minE/893912‑894038

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: