Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,573 C→A 14.0% A100A (GCG→GCT) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5730CA14.0% 356.7 / 25.2 120A100A (GCG→GCT) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (51/52);  new base A (10/7);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

CTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGT  >  minE/312504‑312643
                                                                     |                                                                      
ctgAACTTCTTGCGGCAGCACCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGCACCTGCGCGATATCCg                                                                       >  1:2945349/1‑71 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                    >  1:515958/1‑70 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                    >  1:4642006/1‑70 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                    >  1:448469/1‑70 (MQ=255)
        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                 <  1:4296053/69‑1 (MQ=255)
          tGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                              >  1:4264842/1‑70 (MQ=255)
           gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                                               <  1:870112/68‑2 (MQ=16)
           gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                                            >  1:233171/1‑70 (MQ=16)
           gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                              <  1:5343728/69‑1 (MQ=255)
           gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                              >  1:2560355/1‑69 (MQ=255)
            cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccgg                                                           <  1:4162183/71‑3 (MQ=16)
            cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                                            <  1:4544287/70‑2 (MQ=16)
            cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                            <  1:3411071/70‑1 (MQ=255)
                agcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                       <  1:4679745/71‑2 (MQ=16)
                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                      <  1:4628002/71‑1 (MQ=255)
                  caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                                          <  1:3765029/66‑1 (MQ=255)
                  caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                     <  1:1095740/71‑1 (MQ=255)
                   aaCGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtgg                                                      >  1:4200181/1‑67 (MQ=12)
                   aaCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                     <  1:3928987/70‑1 (MQ=255)
                     cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                                  >  1:5432914/1‑71 (MQ=255)
                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                   >  1:3842056/1‑69 (MQ=255)
                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                 <  1:1776993/71‑1 (MQ=255)
                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                 >  1:1006632/1‑71 (MQ=255)
                         aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                              >  1:4784177/1‑71 (MQ=255)
                          tCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                                             >  1:4466570/1‑71 (MQ=255)
                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                           >  1:5243350/1‑71 (MQ=255)
                               aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                         <  1:3901193/70‑1 (MQ=255)
                               aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                        <  1:3479724/71‑1 (MQ=255)
                                gCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                                               >  1:5688232/1‑46 (MQ=16)
                                   ggaGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGAc                                                         >  1:1986088/3‑50 (MQ=12)
                                   gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                                  >  1:3125337/1‑41 (MQ=255)
                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGc                                    <  1:4687738/70‑1 (MQ=255)
                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                   <  1:3328822/71‑1 (MQ=255)
                                     aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                              >  1:444222/1‑59 (MQ=255)
                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                   <  1:165212/52‑1 (MQ=255)
                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                                            <  1:208285/59‑1 (MQ=255)
                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                                            <  1:159424/59‑1 (MQ=255)
                                         tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                                            >  1:1841353/1‑40 (MQ=12)
                                         tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                                            >  1:5196927/1‑40 (MQ=12)
                                         tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTccc                                                             >  1:2778589/1‑40 (MQ=12)
                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                      <  1:809661/47‑1 (MQ=255)
                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                              >  1:1367844/1‑71 (MQ=255)
                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                              >  1:2388103/1‑71 (MQ=255)
                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                              >  1:4119164/1‑71 (MQ=255)
                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                              >  1:2858834/1‑71 (MQ=255)
                                           ttCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGTAATTGTCACGCTAcc                            >  1:183813/1‑71 (MQ=255)
                                            tCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGTAATTGTCACGCTAccg                           >  1:3060230/1‑70 (MQ=255)
                                            tCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                                <  1:1061429/66‑1 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                           <  1:3835512/70‑1 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                           <  1:5346171/70‑1 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                           >  1:2931164/1‑70 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                          >  1:515794/1‑71 (MQ=255)
                                              tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                                               <  1:2099292/49‑1 (MQ=255)
                                               gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                        >  1:2156103/1‑71 (MQ=255)
                                                cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                                             >  1:1186796/1‑49 (MQ=255)
                                                cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                                             >  1:1963865/1‑49 (MQ=255)
                                                cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                        <  1:4528965/70‑1 (MQ=255)
                                                  ccTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                                       <  1:4420204/53‑1 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                                       >  1:4096186/1‑51 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                                       >  1:1548633/1‑51 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                                       >  1:342486/1‑51 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                                       >  1:5405150/1‑51 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                                       >  1:2754127/1‑51 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                    >  1:4291963/1‑70 (MQ=255)
                                                     gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                  >  1:5319846/1‑71 (MQ=255)
                                                       aCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                          <  1:5263897/45‑1 (MQ=255)
                                                       aCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                          <  1:2689198/45‑1 (MQ=255)
                                                         cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg                       <  1:2092730/62‑1 (MQ=255)
                                                         cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg                       <  1:3187097/62‑1 (MQ=255)
                                                         cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg                       <  1:3717452/62‑1 (MQ=255)
                                                         cTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg                       <  1:2477612/62‑1 (MQ=255)
                                                          tGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                                               <  1:2324496/37‑1 (MQ=255)
                                                          tGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt              >  1:3784515/1‑70 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTc                                                >  1:5035813/1‑35 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                   >  1:5584548/1‑48 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt              >  1:1688633/1‑69 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt              >  1:5504038/1‑69 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt              >  1:887521/1‑69 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt              >  1:2221081/1‑69 (MQ=255)
                                                            cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt              >  1:3921548/1‑68 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                              <  1:4701024/47‑1 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTCCCCGTgg                       >  1:3660772/1‑54 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                    <  1:3748658/57‑1 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                    <  1:66138/57‑1 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                    <  1:4944176/57‑1 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag           >  1:5381166/1‑66 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa       <  1:4340702/70‑1 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa       >  1:2985424/1‑70 (MQ=255)
                                                                  atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGtt             <  1:68859/63‑1 (MQ=255)
                                                                   tCCTCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                     <  1:4604137/38‑1 (MQ=255)
                                                                   tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                     <  1:2836004/38‑1 (MQ=255)
                                                                   tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                     <  1:5204235/38‑1 (MQ=255)
                                                                   tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGc                                    <  1:3000640/39‑1 (MQ=255)
                                                                   tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAcc                            >  1:2715166/1‑47 (MQ=255)
                                                                   tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACAt    <  1:4856884/71‑1 (MQ=255)
                                                                     cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg   <  1:2920438/70‑1 (MQ=255)
                                                                     cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg   <  1:1534703/70‑1 (MQ=255)
                                                                     cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg   <  1:4864189/70‑1 (MQ=255)
                                                                     cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg   <  1:5139095/70‑1 (MQ=255)
                                                                     cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGt  <  1:4151029/71‑1 (MQ=255)
                                                                     |                                                                      
CTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGT  >  minE/312504‑312643

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: