Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,895 T→A 9.4% P326P (CCA→CCT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8950TA 9.4% 174.0 / 9.0 64P326P (CCA→CCTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (20/38);  new base A (4/2);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCG  >  minE/311828‑311962
                                                                   |                                                                   
tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg                                                                  >  1:3992616/1‑71 (MQ=255)
  atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg                                                                  >  1:3219919/4‑69 (MQ=12)
  atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg                                                                  >  1:3060770/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)
     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:4342279/71‑1 (MQ=255)
           aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:286767/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)
             cAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCa                                                     <  1:3497376/71‑1 (MQ=255)
               gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                                   <  1:5143234/71‑1 (MQ=255)
               gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                                   <  1:1235295/71‑1 (MQ=255)
                  ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                                        <  1:742067/63‑1 (MQ=255)
                  ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                                        <  1:2709341/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)
                      aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACg                                            >  1:3010050/1‑71 (MQ=255)
                           cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTc                                                      >  1:3240829/1‑56 (MQ=255)
                           cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                                       <  1:4687316/71‑1 (MQ=255)
                           cTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTAGAGAACCTTCATTCCCTGCGGgaag                                       <  1:1628984/71‑1 (MQ=255)
                              cGTGAATAGAGATTTTCACGAACGGCGTGGTGACGTATgg                                                                   >  1:2802816/1‑40 (MQ=255)
                               gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                                                     >  1:1221363/1‑37 (MQ=255)
                               gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                                                     >  1:587779/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:4379600/47‑1 (MQ=255)
                                          ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:3998361/47‑1 (MQ=255)
                                          ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:1177609/47‑1 (MQ=255)
                                          ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:3293815/47‑1 (MQ=255)
                                          ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:22871/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)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:5685646/35‑1 (MQ=255)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:3068489/35‑1 (MQ=255)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:4524394/35‑1 (MQ=255)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAg                   <  1:1618606/71‑1 (MQ=255)
                                                 gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGtt                 <  1:91683/71‑1 (MQ=255)
                                                 gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGtt                 >  1:1724388/1‑71 (MQ=255)
                                                 gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGtt                 <  1:1909821/71‑1 (MQ=255)
                                                   aCGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGc                      >  1:4850607/1‑64 (MQ=255)
                                                    cGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAg              <  1:84959/71‑1 (MQ=255)
                                                      gCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAAc                                    >  1:1707413/1‑47 (MQ=255)
                                                      gCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAAc                                    >  1:4015709/1‑47 (MQ=255)
                                                        gtggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTc                               >  1:2605170/1‑50 (MQ=255)
                                                           gtgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGt                  >  1:3636944/1‑60 (MQ=255)
                                                           gtgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGAt       <  1:2725965/71‑1 (MQ=255)
                                                            tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGAtt      <  1:4266013/71‑1 (MQ=255)
                                                            tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGAtt      <  1:2673216/71‑1 (MQ=255)
                                                             gtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCtt                         >  1:1919033/1‑51 (MQ=255)
                                                              tCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGc    <  1:172392/71‑1 (MQ=255)
                                                                gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCAc           >  1:713368/1‑62 (MQ=255)
                                                                gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCg  <  1:2286793/71‑1 (MQ=255)
                                                                gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCg  <  1:489665/71‑1 (MQ=255)
                                                                 tATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCg  <  1:3289451/70‑1 (MQ=255)
                                                                   tGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAg              <  1:3183231/56‑1 (MQ=255)
                                                                   tGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACg          <  1:4887912/60‑1 (MQ=255)
                                                                   |                                                                   
TGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCG  >  minE/311828‑311962

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: