Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,074 C→G 9.2% A401A (GCG→GCC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0740CG 9.2% 358.3 / 14.9 118A401A (GCG→GCCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (44/63);  new base G (6/5);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTC  >  minE/2586004‑2586143
                                                                      |                                                                     
   tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                                     >  1:3234193/1‑70 (MQ=255)
   tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                                     >  1:1274685/1‑70 (MQ=255)
    aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCg                                                                      <  1:2119080/68‑1 (MQ=255)
    aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                                                    <  1:3554745/70‑1 (MQ=255)
     tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCa                                                                   <  1:2884450/70‑1 (MQ=255)
       ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATACACCAGTAACCCTATCGCCAGCa                                                                >  1:3488866/1‑71 (MQ=255)
        tCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                               >  1:521560/1‑71 (MQ=255)
         ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                                <  1:1871226/69‑1 (MQ=255)
         ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                                <  1:3080823/69‑1 (MQ=255)
          cGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                                              <  1:2019901/70‑1 (MQ=255)
          cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                                              <  1:59734/70‑1 (MQ=255)
          cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                                             <  1:3776617/71‑1 (MQ=255)
             aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                          >  1:3364156/1‑71 (MQ=255)
             aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                          >  1:253318/1‑71 (MQ=255)
               tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                                       >  1:1467828/1‑56 (MQ=255)
               tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                                       >  1:3748009/1‑56 (MQ=255)
               tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                                         >  1:416224/1‑70 (MQ=255)
                gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                               <  1:786780/63‑1 (MQ=255)
                 acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATGCCGaa                                                      >  1:2101319/4‑71 (MQ=2)
                 acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                                      >  1:1833409/4‑71 (MQ=2)
                 acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                                      >  1:3790593/4‑71 (MQ=2)
                    cGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACa                                                    >  1:1706856/1‑70 (MQ=2)
                    cGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                   <  1:3795494/71‑1 (MQ=255)
                     gCTCGACGTTCTCCACCACCGCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                                        >  1:3353477/1‑65 (MQ=255)
                     gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                                <  1:148099/57‑1 (MQ=255)
                     gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                                <  1:1200652/57‑1 (MQ=255)
                     gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                                <  1:3711481/57‑1 (MQ=255)
                     gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                    >  1:2023796/1‑69 (MQ=255)
                     gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                    >  1:944151/1‑69 (MQ=255)
                     gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                    >  1:1872219/1‑69 (MQ=255)
                        cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                                     <  1:928251/65‑1 (MQ=255)
                         gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                                <  1:2651753/69‑1 (MQ=255)
                         gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                               <  1:2405175/70‑1 (MQ=255)
                          aCGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAg                                                 >  1:2130414/1‑67 (MQ=255)
                           cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                                     >  1:2417964/1‑62 (MQ=255)
                           cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                            <  1:3045170/71‑1 (MQ=255)
                           cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                            >  1:417058/1‑71 (MQ=255)
                            gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                           >  1:1140588/1‑71 (MQ=255)
                            gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                           <  1:2441673/71‑1 (MQ=255)
                            gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                           <  1:2210775/71‑1 (MQ=255)
                            gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                           <  1:2033255/71‑1 (MQ=255)
                             ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                            <  1:1396246/69‑1 (MQ=255)
                             ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                          >  1:1738287/1‑71 (MQ=255)
                              tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                                         >  1:142341/1‑71 (MQ=255)
                              tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                          <  1:526951/70‑1 (MQ=255)
                                   accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                   <  1:836523/56‑1 (MQ=255)
                                   accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                   <  1:1956328/56‑1 (MQ=255)
                                   accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                    >  1:3481322/1‑71 (MQ=255)
                                   accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                    <  1:916688/71‑1 (MQ=255)
                                   accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                    >  1:3757272/1‑71 (MQ=255)
                                   accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                    >  1:2685255/1‑71 (MQ=255)
                                     caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                          <  1:1421996/63‑1 (MQ=255)
                                        cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                               <  1:2792971/71‑1 (MQ=255)
                                        cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                               <  1:3618446/71‑1 (MQ=255)
                                         acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              <  1:40117/71‑1 (MQ=255)
                                         acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              <  1:1450063/71‑1 (MQ=255)
                                         acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              <  1:3900925/71‑1 (MQ=255)
                                         acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              <  1:2680586/71‑1 (MQ=255)
                                          cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                             <  1:2915695/71‑1 (MQ=255)
                                          cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                             <  1:484202/71‑1 (MQ=255)
                                           gATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                               >  1:124831/1‑52 (MQ=255)
                                            aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                           <  1:3099023/71‑1 (MQ=255)
                                            aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                           <  1:3423260/71‑1 (MQ=255)
                                             tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccg                          <  1:3864657/71‑1 (MQ=255)
                                                 gTCATCCACCAGTAACCCTATCGCCAGCACCATCTCGAACATAGTCAAc                                            <  1:3740599/49‑1 (MQ=255)
                                                  tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATg                                   <  1:2577487/57‑1 (MQ=255)
                                                  tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc                     <  1:3925589/71‑1 (MQ=255)
                                                    aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                   >  1:178616/1‑39 (MQ=255)
                                                    aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAac                            >  1:212543/1‑62 (MQ=255)
                                                    aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                    >  1:1729638/1‑70 (MQ=255)
                                                    aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                    <  1:1155688/70‑1 (MQ=255)
                                                    aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                    >  1:1610756/1‑70 (MQ=255)
                                                     tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGa                  >  1:3519300/1‑71 (MQ=255)
                                                       caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc                 <  1:1942710/70‑1 (MQ=255)
                                                       caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc                 <  1:923957/70‑1 (MQ=255)
                                                        accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc                     <  1:625019/65‑1 (MQ=255)
                                                        accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACaa               <  1:673907/71‑1 (MQ=255)
                                                         ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                               <  1:2271318/38‑1 (MQ=255)
                                                         ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc                 >  1:2024760/1‑68 (MQ=255)
                                                         ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAg              <  1:694486/71‑1 (MQ=255)
                                                          caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                                     <  1:738470/47‑1 (MQ=255)
                                                          caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                             >  1:417796/1‑55 (MQ=255)
                                                          caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa             <  1:852538/71‑1 (MQ=255)
                                                          caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa             <  1:2857097/71‑1 (MQ=255)
                                                          caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa             <  1:3929530/71‑1 (MQ=255)
                                                          caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa             <  1:1782093/71‑1 (MQ=255)
                                                           aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCaa                                             >  1:2509047/1‑38 (MQ=255)
                                                           aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                                      >  1:146125/1‑45 (MQ=255)
                                                           aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                                      >  1:664354/1‑45 (MQ=255)
                                                               aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCa        <  1:3838294/71‑1 (MQ=255)
                                                                 ccTATCGCCTGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCaa       <  1:2559852/70‑1 (MQ=255)
                                                                  cTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                               <  1:1406112/45‑1 (MQ=255)
                                                                   tATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACg    >  1:1318356/1‑71 (MQ=255)
                                                                    aTCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGt   >  1:1824727/1‑71 (MQ=255)
                                                                     tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATTGTGAAACCGACCGCCGACAAGATCGCAAACGt   <  1:1570086/70‑1 (MQ=255)
                                                                     tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              >  1:1935721/1‑43 (MQ=255)
                                                                     tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              >  1:2267649/1‑43 (MQ=255)
                                                                     tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                              >  1:3305288/1‑43 (MQ=255)
                                                                     tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                               <  1:1491114/42‑1 (MQ=255)
                                                                     tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGt   <  1:2014738/70‑1 (MQ=255)
                                                                      |                                                                     
GCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTC  >  minE/2586004‑2586143

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: