Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 706,680 C→T 11.5% R244K (AGG→AAG)  lpxL ← lauryl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE706,6800CT11.5% 175.0 / 11.4 61R244K (AGG→AAG) lpxLlauryl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/21);  new base T (4/3);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

TCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCGTGGT  >  minE/706611‑706750
                                                                     |                                                                      
tCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGcc                                                                        <  1:235109/70‑1 (MQ=255)
        atcCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCg                                                                 <  1:3346128/69‑1 (MQ=255)
        atcCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCg                                                                 <  1:3149342/69‑1 (MQ=255)
        atcCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTg                                                               >  1:923800/1‑71 (MQ=255)
        atcCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTg                                                               >  1:2788814/1‑71 (MQ=255)
          ccAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTgg                                                              >  1:414684/1‑70 (MQ=255)
            aGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAAc                                                           >  1:2891093/1‑71 (MQ=255)
            aGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAAc                                                           >  1:1994110/1‑71 (MQ=255)
             gTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAAc                                                           >  1:2686141/1‑70 (MQ=255)
             gTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACg                                                          <  1:1378853/71‑1 (MQ=255)
             gTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACg                                                          <  1:2644544/71‑1 (MQ=255)
             gTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACg                                                          <  1:4006075/71‑1 (MQ=255)
              tGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTa                                                                      >  1:840960/1‑58 (MQ=255)
              tGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGa                                                         <  1:1013544/71‑1 (MQ=255)
              tGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGa                                                         <  1:1167868/71‑1 (MQ=255)
              tGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGa                                                         <  1:87849/71‑1 (MQ=255)
                gCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGag                                                        >  1:2528440/1‑69 (MQ=255)
                     gaACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGa                                                             >  1:3640075/1‑60 (MQ=255)
                      aaCACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGaa                                                        <  1:369782/64‑1 (MQ=255)
                       acacTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCa                                                <  1:2265718/71‑1 (MQ=255)
                        cacTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCAc                                                  >  1:2014142/1‑68 (MQ=255)
                        cacTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCAc                                                  >  1:2958471/1‑68 (MQ=255)
                        cacTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAg                                               <  1:1776823/71‑1 (MQ=255)
                          ctctGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGa                                                         <  1:2297900/59‑1 (MQ=255)
                          ctctGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGa                                              <  1:2157698/70‑1 (MQ=255)
                            ctGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGaa                                                        >  1:1776233/1‑58 (MQ=255)
                             tGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATg                                          >  1:1621795/1‑71 (MQ=255)
                                   cagcaTAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCATACAt                                           >  1:1992379/1‑64 (MQ=255)
                                   cagcaTAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAg                                                       <  1:3758644/52‑1 (MQ=255)
                                        tAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACAt                                >  1:4014513/1‑70 (MQ=255)
                                        tAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATc                               <  1:3875647/71‑1 (MQ=255)
                                        tAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATc                               <  1:157499/71‑1 (MQ=255)
                                          aTCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGAc                                             >  1:247862/1‑55 (MQ=255)
                                          aTCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCg                             >  1:3632104/1‑71 (MQ=255)
                                           tCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGaa                                                        <  1:2031499/43‑1 (MQ=255)
                                           tCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGt                            <  1:1229802/71‑1 (MQ=255)
                                            cAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGaa                                                            >  1:1845877/1‑38 (MQ=255)
                                            cAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGa                                              <  1:1121280/52‑1 (MQ=255)
                                            cAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCg                             >  1:558273/1‑69 (MQ=255)
                                             aaTTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGc                          >  1:657242/1‑71 (MQ=255)
                                             aaTTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACa                                            >  1:128316/1‑53 (MQ=255)
                                             aaTTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACa                                            >  1:224247/1‑53 (MQ=255)
                                             aaTTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACa                                            >  1:831138/1‑53 (MQ=255)
                                             aaTTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTg                           <  1:1144405/70‑1 (MQ=255)
                                                 gATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAg                       >  1:1307840/1‑70 (MQ=255)
                                                   tACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCAt                    <  1:2325480/71‑1 (MQ=255)
                                                    aCCCTTTGCCATCTTGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc                   >  1:2845442/1‑71 (MQ=255)
                                                    aCCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCAcc                                                 <  1:757077/41‑1 (MQ=255)
                                                     cccTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTg                           >  1:2314176/1‑62 (MQ=255)
                                                     cccTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc                   >  1:838207/1‑70 (MQ=255)
                                                      ccTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc                   >  1:1037803/1‑69 (MQ=255)
                                                      ccTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc                   <  1:3842795/69‑1 (MQ=255)
                                                      ccTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc                   >  1:3558424/1‑69 (MQ=255)
                                                      ccTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc                   >  1:3553429/1‑69 (MQ=255)
                                                        tttGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGa                                              >  1:1288545/1‑40 (MQ=255)
                                                              aTCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCgg         >  1:2490935/1‑71 (MQ=255)
                                                                 tGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACAt                                >  1:3616668/1‑45 (MQ=255)
                                                                 tGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCg      >  1:3863958/1‑71 (MQ=255)
                                                                 tGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCg      >  1:917545/1‑71 (MQ=255)
                                                                 tGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCg      >  1:2700048/1‑71 (MQ=255)
                                                                 tGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGCTTCCGGTCg      >  1:1910805/1‑71 (MQ=255)
                                                                   gCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCg      <  1:3878132/69‑1 (MQ=255)
                                                                     cTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCGTGGt  <  1:3001650/71‑1 (MQ=255)
                                                                     |                                                                      
TCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCCTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATCCAGGTTCCGGTCGTGGT  >  minE/706611‑706750

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: