Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,065 T→C 6.7% L404L (TTA→TTG) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0650TC 6.7% 297.8 / 10.7 104L404L (TTA→TTG) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (41/56);  new base C (3/4);  total (44/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 112 total aligned reads displayed.

CGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAT  >  minE/2585996‑2586133
                                                                     |                                                                    
cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                                      <  1:2144167/70‑1 (MQ=255)
  gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                                   <  1:3349906/71‑1 (MQ=255)
    gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                                 <  1:3071193/68‑1 (MQ=9)
    gggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                                 <  1:626892/68‑1 (MQ=9)
     gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                                 >  1:790588/1‑70 (MQ=255)
      tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCCCCAGTAACCCTa                                                               >  1:1636305/1‑71 (MQ=255)
      tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                                                <  1:52052/70‑1 (MQ=255)
      tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                               <  1:3260103/71‑1 (MQ=255)
      gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga                                                               <  1:4053857/70‑3 (MQ=255)
      gaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGcccga                                                               >  1:222344/2‑69 (MQ=255)
       aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                                                >  1:4113399/1‑69 (MQ=255)
       aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                              >  1:1670558/1‑71 (MQ=255)
         cTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                                   <  1:1695424/64‑1 (MQ=255)
           tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                           >  1:920585/1‑70 (MQ=255)
            aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                                      >  1:1107565/1‑58 (MQ=255)
            aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                           >  1:899383/1‑69 (MQ=255)
             tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCt                                                                <  1:3825497/63‑1 (MQ=255)
              cTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                        >  1:2156977/1‑70 (MQ=255)
               ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                      >  1:3505798/1‑71 (MQ=255)
                 ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                        <  1:1017283/67‑1 (MQ=255)
                 ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                        <  1:2466172/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:3473991/67‑1 (MQ=255)
                      aTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                                   <  1:1663049/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)
                         acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                     >  1:3560893/1‑62 (MQ=255)
                         acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                            >  1:501919/1‑71 (MQ=255)
                         acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                            >  1:3412148/1‑71 (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:3828203/1‑67 (MQ=255)
                                       ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                  >  1:2643675/1‑67 (MQ=255)
                                       ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                               <  1:2930515/70‑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)
                                               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:1635053/71‑1 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                     >  1:2046148/1‑70 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                     >  1:3312876/1‑70 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACAATCCCGAACATAGTCAACGTGTTGATGGTGa                     >  1:3819739/1‑70 (MQ=255)
                                                 acGATGGCGTCATCCACCACTAACCCTATCGACAGCAAGATCCCGAACATAGTCAACGTGTAGAAGGt                       >  1:924041/1‑68 (MQ=25)
                                                     tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                    <  1:35588/67‑1 (MQ=255)
                                                     tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                    <  1:3964587/67‑1 (MQ=255)
                                                     tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                 >  1:2920476/1‑70 (MQ=255)
                                                       gCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgac              >  1:3746943/1‑71 (MQ=255)
                                                        cGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                  >  1:4051349/1‑50 (MQ=255)
                                                        cGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgacc             >  1:1517052/1‑71 (MQ=255)
                                                        cGTCATCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgacc             <  1:3706088/70‑1 (MQ=255)
                                                          tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                             <  1:2191374/53‑1 (MQ=255)
                                                           cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc          >  1:1595167/1‑71 (MQ=255)
                                                           cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc          >  1:2909912/1‑71 (MQ=255)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccgac              >  1:338247/1‑66 (MQ=255)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc          >  1:3148884/1‑70 (MQ=255)
                                                             tCCACCAGTAACCCTATCGCTAGCACCATCCCGAACATAGTc                                     <  1:2491892/42‑1 (MQ=255)
                                                             tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                               <  1:2125051/48‑1 (MQ=255)
                                                                accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACaa     <  1:1224088/71‑1 (MQ=255)
                                                                 ccaGTAACCCTATCGCCAGCACCATCCCGAACATa                                        <  1:3934213/35‑1 (MQ=255)
                                                                 ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACGGCCGACAAg    <  1:213206/71‑1 (MQ=255)
                                                                 ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAg    <  1:3971846/71‑1 (MQ=255)
                                                                  caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa   <  1:954895/71‑1 (MQ=255)
                                                                  caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa   <  1:2940710/71‑1 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                  >  1:2026449/1‑39 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                  >  1:1037108/1‑39 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg                      >  1:2995348/1‑51 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg                      >  1:1754906/1‑51 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc       >  1:2816409/1‑66 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc       >  1:1342981/1‑66 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAt  <  1:3074936/71‑1 (MQ=255)
                                                                    gTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                           <  1:3167683/45‑1 (MQ=255)
                                                                     tAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                             <  1:3786760/42‑1 (MQ=255)
                                                                     tAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                             <  1:1538144/42‑1 (MQ=255)
                                                                     tAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTg                      >  1:3734539/1‑49 (MQ=255)
                                                                     tAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACaa     <  1:3942450/66‑1 (MQ=255)
                                                                     tAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa   >  1:3215870/1‑68 (MQ=255)
                                                                     |                                                                    
CGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAT  >  minE/2585996‑2586133

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: