Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,873 T→G 16.0% K334Q (AAA→CAA) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8730TG16.0% 201.7 / 36.5 94K334Q (AAA→CAA) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (36/43);  new base G (8/7);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

ATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCT  >  minE/311805‑311938
                                                                    |                                                                 
aTACATAACGAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                                                 >  1:1882947/1‑71 (MQ=255)
  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                                                >  1:4986943/1‑70 (MQ=255)
   cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                                                              <  1:2122378/71‑1 (MQ=255)
   cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                                                              <  1:1791893/71‑1 (MQ=255)
   cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                                                              <  1:4763334/71‑1 (MQ=255)
    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                                                 >  1:4590069/1‑67 (MQ=255)
    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTACGTGAATAGAGATTTTCa                                                                 >  1:4825532/1‑67 (MQ=255)
     tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                                                            <  1:4887055/71‑1 (MQ=255)
          aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                       <  1:5323838/71‑1 (MQ=255)
          aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                       <  1:5056854/71‑1 (MQ=255)
           ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                      >  1:1388699/1‑71 (MQ=255)
           ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                      >  1:2177368/1‑71 (MQ=255)
           ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                                                      >  1:4970918/1‑71 (MQ=255)
           ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                                                      <  1:3748263/71‑1 (MQ=255)
            gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                                                   <  1:326925/57‑1 (MQ=255)
            gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                                                   <  1:1019905/57‑1 (MQ=255)
            gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg                                                     >  1:4428081/1‑71 (MQ=255)
            gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg                                                     >  1:5289143/1‑71 (MQ=255)
             aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                                    >  1:3048319/1‑71 (MQ=255)
             aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                                    <  1:2953439/71‑1 (MQ=255)
             aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                      >  1:4267877/1‑69 (MQ=255)
             aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                      >  1:3123258/1‑69 (MQ=255)
             aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                      >  1:2557573/1‑69 (MQ=255)
             aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATATTCACGAACGGCgtg                                                      >  1:3374507/1‑69 (MQ=255)
              acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                                                   <  1:4196140/71‑1 (MQ=255)
              acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                                                   >  1:4649374/1‑71 (MQ=255)
                acGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGgtgt                                                  <  1:1719512/70‑1 (MQ=255)
                 cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                                                   >  1:1144417/1‑52 (MQ=255)
                 cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                                                                  >  1:4707099/1‑53 (MQ=255)
                 cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                                                 >  1:1715554/1‑54 (MQ=255)
                  gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                >  1:2362560/1‑70 (MQ=255)
                   aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                                                 >  1:1313830/1‑52 (MQ=255)
                   aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                              <  1:3480806/71‑1 (MQ=255)
                   aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGAGTGGTGTCGTa                                              >  1:2278972/1‑71 (MQ=255)
                       tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                                                           <  1:2902771/54‑1 (MQ=255)
                       tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                       >  1:4893259/1‑58 (MQ=255)
                       tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg                                          >  1:3992616/1‑71 (MQ=255)
                        gatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                                               <  1:2354159/49‑1 (MQ=255)
                         atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg                                          >  1:3060770/4‑69 (MQ=12)
                         atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg                                          >  1:3219919/4‑69 (MQ=12)
                         atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa                                        <  1:1614961/68‑1 (MQ=16)
                         atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACTGCGTGGTGTCGTATGGGTa                                        <  1:2824688/71‑1 (MQ=255)
                         atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa                                        <  1:3785131/71‑1 (MQ=255)
                          tggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg                                          >  1:132925/3‑68 (MQ=16)
                          tCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                                             >  1:2947451/1‑49 (MQ=255)
                            gCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                                                      >  1:614230/1‑54 (MQ=255)
                            gCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa                                        >  1:4444151/1‑68 (MQ=16)
                             cTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAAc                                     >  1:3778978/1‑70 (MQ=255)
                             cTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACa                                    <  1:5563301/71‑1 (MQ=255)
                              ttCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                   >  1:3547085/1‑71 (MQ=255)
                               tCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAt                                  >  1:3574996/1‑71 (MQ=255)
                               tCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAt                                  >  1:4448835/1‑71 (MQ=255)
                                cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                   <  1:2911483/69‑1 (MQ=255)
                                  aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                               <  1:286767/71‑1 (MQ=255)
                                  aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                               <  1:4342279/71‑1 (MQ=255)
                                  aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                               <  1:2827703/71‑1 (MQ=255)
                                  aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                >  1:5557755/1‑70 (MQ=255)
                                  aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGAAAACAAtttt                               <  1:1681204/71‑1 (MQ=255)
                                    cAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                                               <  1:5513063/37‑1 (MQ=255)
                                    cAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCa                             <  1:3497376/71‑1 (MQ=255)
                                     aGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:3461784/51‑1 (MQ=255)
                                      gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                           <  1:1235295/71‑1 (MQ=255)
                                      gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                           <  1:5143234/71‑1 (MQ=255)
                                       cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                              >  1:3057517/1‑51 (MQ=255)
                                       cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                              >  1:2735792/1‑51 (MQ=255)
                                         ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTcat                                               <  1:4251544/48‑3 (MQ=12)
                                         ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTcat                                               <  1:4839119/48‑3 (MQ=12)
                                         ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTcat                                               <  1:1817670/48‑3 (MQ=12)
                                         ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                <  1:2709341/63‑1 (MQ=255)
                                         ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                <  1:742067/63‑1 (MQ=255)
                                         ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtt                                 <  1:3855210/62‑1 (MQ=255)
                                         ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                        <  1:2905741/71‑1 (MQ=255)
                                           ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                             >  1:1164868/1‑48 (MQ=255)
                                           ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                             >  1:4541384/1‑48 (MQ=255)
                                             aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:668516/43‑1 (MQ=255)
                                             aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACg                    >  1:3010050/1‑71 (MQ=255)
                                              aCCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                       <  1:5637196/35‑1 (MQ=255)
                                               ccACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:2750151/41‑1 (MQ=255)
                                                cACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:439659/40‑1 (MQ=255)
                                                  cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:1360000/38‑1 (MQ=255)
                                                  cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:1152826/38‑1 (MQ=255)
                                                  cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTc                              >  1:3240829/1‑56 (MQ=255)
                                                  cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag               <  1:4687316/71‑1 (MQ=255)
                                                  cTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTAGAGAACCTTCATTCCCTGCGGgaag               <  1:1628984/71‑1 (MQ=255)
                                                   ttCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                <  1:3032935/37‑1 (MQ=255)
                                                     cGTGAATAGAGATTTTCACGAACGGCGTGGTGACGTATgg                                           >  1:2802816/1‑40 (MQ=255)
                                                      gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                              >  1:2875209/1‑36 (MQ=255)
                                                      gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                             >  1:587779/1‑37 (MQ=255)
                                                      gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                             >  1:1221363/1‑37 (MQ=255)
                                                           tAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTcc      <  1:1808024/71‑1 (MQ=255)
                                                             gagaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCAt    >  1:743894/1‑71 (MQ=255)
                                                               gaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCt  <  1:5619689/71‑1 (MQ=255)
                                                                 ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                        <  1:3293815/47‑1 (MQ=255)
                                                                 ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                        <  1:22871/47‑1 (MQ=255)
                                                                 ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                        <  1:3998361/47‑1 (MQ=255)
                                                                 ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                        <  1:4379600/47‑1 (MQ=255)
                                                                 ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                        <  1:1177609/47‑1 (MQ=255)
                                                                 ttttCACGAACGGCGTGGTGTCGTATGGGCAAACAATTTTCAGAccc                        <  1:5313097/47‑1 (MQ=255)
                                                                   ttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                           >  1:4099146/1‑42 (MQ=255)
                                                                   ttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                           >  1:5069204/1‑42 (MQ=255)
                                                                    |                                                                 
ATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCT  >  minE/311805‑311938

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: