Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,892 G→A 10.0% Y327Y (TAC→TAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8920GA10.0% 199.7 / 9.0 70Y327Y (TAC→TATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base A (4/3);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCG  >  minE/311823‑311962
                                                                     |                                                                      
gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                        >  1:2362560/1‑70 (MQ=255)
 aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                      <  1:3480806/71‑1 (MQ=255)
 aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGAGTGGTGTCGTa                                                                      >  1:2278972/1‑71 (MQ=255)
     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:4448835/1‑71 (MQ=255)
             tCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAt                                                          >  1:3574996/1‑71 (MQ=255)
              cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                                           <  1:2911483/69‑1 (MQ=255)
                aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:4342279/71‑1 (MQ=255)
                aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:2827703/71‑1 (MQ=255)
                aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:286767/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)
                   aGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                        <  1:3461784/51‑1 (MQ=255)
                    gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                                   <  1:5143234/71‑1 (MQ=255)
                    gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                                   <  1:1235295/71‑1 (MQ=255)
                     cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                      >  1:2735792/1‑51 (MQ=255)
                     cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                      >  1:3057517/1‑51 (MQ=255)
                       ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATaag                                                                    <  1:2968612/51‑3 (MQ=9)
                       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:4541384/1‑48 (MQ=255)
                         ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                                                     >  1:1164868/1‑48 (MQ=255)
                           aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                        <  1:668516/43‑1 (MQ=255)
                           aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                        <  1:834660/43‑1 (MQ=255)
                           aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACg                                            >  1:3010050/1‑71 (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: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:3998361/47‑1 (MQ=255)
                                               ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:4379600/47‑1 (MQ=255)
                                               ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:22871/47‑1 (MQ=255)
                                               ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:1177609/47‑1 (MQ=255)
                                               ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                                <  1:3293815/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)
                                                                     |                                                                      
GAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCG  >  minE/311823‑311962

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: