Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,047 G→T 9.6% I410I (ATC→ATA mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0470GT 9.6% 333.2 / 12.1 103I410I (ATC→ATAmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (43/50);  new base T (3/7);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 103 total aligned reads displayed.

TTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGA  >  minE/2585978‑2586114
                                                                     |                                                                   
tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                                                    >  1:2428704/1‑71 (MQ=255)
   aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                                                  <  1:3102437/70‑1 (MQ=255)
      cGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCg                                                               <  1:1850986/70‑1 (MQ=255)
      cGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCg                                                               <  1:2174545/70‑1 (MQ=255)
        tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCa                                                            >  1:3425377/1‑71 (MQ=255)
        tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCa                                                            >  1:548634/1‑71 (MQ=255)
          gCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATc                                                          >  1:2163682/1‑71 (MQ=255)
           cttcttttGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatc                                                          <  1:3797960/70‑1 (MQ=12)
              ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                       >  1:2138510/1‑70 (MQ=255)
              ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                       <  1:3623113/70‑1 (MQ=255)
              ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                       <  1:1032638/70‑1 (MQ=255)
                  cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                   <  1:2144167/70‑1 (MQ=255)
                   ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                 <  1:2388022/71‑1 (MQ=255)
                    gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                <  1:3349906/71‑1 (MQ=255)
                      gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                              <  1:626892/68‑1 (MQ=9)
                      gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                              <  1:3071193/68‑1 (MQ=9)
                       gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                                                                    <  1:661694/48‑1 (MQ=255)
                       gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                              >  1:790588/1‑70 (MQ=255)
                        tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCCCCAGTAACCCTa                                            >  1:1636305/1‑71 (MQ=255)
                        tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                             <  1:52052/70‑1 (MQ=255)
                        tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                            >  1:767643/1‑71 (MQ=255)
                        tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                            <  1:3260103/71‑1 (MQ=255)
                        gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga                                            <  1:4053857/70‑3 (MQ=255)
                         aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAg                                                    >  1:1719769/1‑62 (MQ=255)
                         aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                             >  1:4113399/1‑69 (MQ=255)
                         aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                           >  1:1670558/1‑71 (MQ=255)
                          gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCGCGATGGCGTCATcca                                                        >  1:1881326/1‑57 (MQ=255)
                          gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                                        >  1:1924630/1‑57 (MQ=255)
                          gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                                        >  1:579385/1‑57 (MQ=255)
                          gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATc                                                          >  1:1506712/1‑55 (MQ=255)
                           cTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                <  1:1695424/64‑1 (MQ=255)
                           cTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatcg                                                         <  1:3896454/55‑2 (MQ=255)
                           cTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCatcatcg                                                         <  1:3377221/55‑2 (MQ=255)
                             tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                        >  1:920585/1‑70 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCCCCACGATGGCGTCATccacca                                                     >  1:837820/1‑56 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                                                                <  1:2556031/45‑1 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                                                                <  1:3594923/45‑1 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                                                                <  1:176843/45‑1 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                   >  1:1107565/1‑58 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                        >  1:1404749/1‑69 (MQ=255)
                              aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                        >  1:899383/1‑69 (MQ=255)
                               tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                             <  1:3825497/63‑1 (MQ=255)
                                cTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGt                                                              >  1:4047199/1‑45 (MQ=255)
                                cTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                     >  1:2156977/1‑70 (MQ=255)
                                 ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATc                                                          <  1:562592/48‑1 (MQ=255)
                                 ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                   >  1:3505798/1‑71 (MQ=255)
                                   ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                     <  1:2466172/67‑1 (MQ=255)
                                   ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                     <  1:1017283/67‑1 (MQ=255)
                                   ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                   <  1:2726256/69‑1 (MQ=255)
                                      cAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCAt                               <  1:3209622/70‑1 (MQ=255)
                                      cAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATc                              <  1:2465905/71‑1 (MQ=255)
                                       aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATc                              <  1:644711/70‑1 (MQ=255)
                                        aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                             <  1:1102717/54‑1 (MQ=255)
                                        aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                <  1:1663049/67‑1 (MQ=255)
                                        aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                <  1:3473991/67‑1 (MQ=255)
                                         tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                            >  1:1348189/1‑70 (MQ=255)
                                         tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                           <  1:814906/71‑1 (MQ=255)
                                           acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                         >  1:2630334/4‑71 (MQ=2)
                                           acacGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                     >  1:2704051/1‑43 (MQ=255)
                                           acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                  >  1:3560893/1‑62 (MQ=255)
                                           acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                         >  1:501919/1‑71 (MQ=255)
                                           acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                         >  1:3412148/1‑71 (MQ=255)
                                             acGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                     >  1:4037653/1‑41 (MQ=255)
                                             acGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                     >  1:3612529/1‑41 (MQ=255)
                                             acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                       <  1:2006510/71‑1 (MQ=255)
                                             acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                       <  1:3927935/71‑1 (MQ=255)
                                                 tCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAAc                        <  1:379082/66‑1 (MQ=9)
                                                  cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                   >  1:3853016/1‑70 (MQ=2)
                                                   gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                  <  1:2081560/70‑1 (MQ=255)
                                                   gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                  <  1:1334357/70‑1 (MQ=255)
                                                   gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                  <  1:3155042/70‑1 (MQ=255)
                                                   gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                 <  1:115142/71‑1 (MQ=255)
                                                   gACGTTCTCCAACACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                  <  1:1914800/70‑1 (MQ=255)
                                                       ttCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                              <  1:1429061/38‑1 (MQ=255)
                                                       ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATa                     <  1:920801/63‑1 (MQ=255)
                                                        tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                             >  1:4138693/1‑54 (MQ=255)
                                                        tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                             >  1:932832/1‑54 (MQ=255)
                                                         ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAGCATACCGAACATAGTCAACgtgt           >  1:3743712/1‑71 (MQ=255)
                                                         ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc               >  1:2643675/1‑67 (MQ=255)
                                                         ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc               >  1:3828203/1‑67 (MQ=255)
                                                         ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg            <  1:2930515/70‑1 (MQ=255)
                                                            caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                        <  1:4077617/55‑1 (MQ=255)
                                                            caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt          <  1:2970651/69‑1 (MQ=255)
                                                            caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg         <  1:1637659/70‑1 (MQ=255)
                                                            caccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa        <  1:667836/71‑1 (MQ=255)
                                                             accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt       <  1:1880151/71‑1 (MQ=255)
                                                              ccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                       <  1:3462353/54‑1 (MQ=255)
                                                              ccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                 <  1:591886/60‑1 (MQ=255)
                                                               caccacGATTGCGTCATCCACAAGTAACCCTATCGCCAGCACCATCCCGAACa                       <  1:1477405/53‑1 (MQ=21)
                                                                accacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa        >  1:1766881/1‑67 (MQ=255)
                                                                 ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGa                          <  1:3112642/48‑1 (MQ=255)
                                                                 ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa        >  1:4011644/1‑66 (MQ=255)
                                                                 ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg   <  1:3833700/71‑1 (MQ=255)
                                                                  cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg   >  1:1881580/1‑70 (MQ=255)
                                                                  cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa  <  1:2403952/71‑1 (MQ=255)
                                                                  cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa  <  1:3836052/71‑1 (MQ=255)
                                                                  cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa  <  1:1635053/71‑1 (MQ=255)
                                                                   acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa  >  1:2046148/1‑70 (MQ=255)
                                                                   acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa  >  1:3312876/1‑70 (MQ=255)
                                                                   acGATGGCGTCATCCACCACTAACCCTATCGACAGCAAGATCCCGAACATAGTCAACGTGTAGAAGGt    >  1:924041/1‑68 (MQ=25)
                                                                     |                                                                   
TTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGA  >  minE/2585978‑2586114

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: