Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,880 C→T 15.9% P331P (CCG→CCA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8800CT15.9% 231.2 / 29.3 94P331P (CCG→CCAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (32/47);  new base T (8/7);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAG  >  minE/311810‑311950
                                                                      |                                                                      
tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                                                                        <  1:4887055/71‑1 (MQ=255)
     aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                                   <  1:5056854/71‑1 (MQ=255)
     aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                                   <  1:5323838/71‑1 (MQ=255)
      ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                                  >  1:2177368/1‑71 (MQ=255)
      ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                                  >  1:1388699/1‑71 (MQ=255)
      ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                                                                  >  1:4970918/1‑71 (MQ=255)
      ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                                                                  <  1:3748263/71‑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:3123258/1‑69 (MQ=255)
        aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                                  >  1:2557573/1‑69 (MQ=255)
        aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                                  >  1:4267877/1‑69 (MQ=255)
        aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATATTCACGAACGGCgtg                                                                  >  1:3374507/1‑69 (MQ=255)
         acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                                                               >  1:4649374/1‑71 (MQ=255)
         acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                                                               <  1:4196140/71‑1 (MQ=255)
           acGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGgtgt                                                              <  1:1719512/70‑1 (MQ=255)
             gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                            >  1:2362560/1‑70 (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)
                    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)
                       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:2827703/71‑1 (MQ=255)
                             aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                           <  1:4342279/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:3057517/1‑51 (MQ=255)
                                  cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                          >  1:2735792/1‑51 (MQ=255)
                                    ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTcat                                                           <  1:1817670/48‑3 (MQ=12)
                                    ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTcat                                                           <  1:4839119/48‑3 (MQ=12)
                                    ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTcat                                                           <  1:4251544/48‑3 (MQ=12)
                                    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)
                                        aaCCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                                   <  1:4880805/36‑1 (MQ=255)
                                        aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                            <  1:834660/43‑1 (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:3998361/47‑1 (MQ=255)
                                                            ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                    <  1:1177609/47‑1 (MQ=255)
                                                            ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                    <  1:22871/47‑1 (MQ=255)
                                                            ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                    <  1:4379600/47‑1 (MQ=255)
                                                            ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAccc                                    <  1:3293815/47‑1 (MQ=255)
                                                            ttttCACGAACGGCGTGGTGTCGTATGGGCAAACAATTTTCAGAccc                                    <  1:5313097/47‑1 (MQ=255)
                                                              ttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                       >  1:5069204/1‑42 (MQ=255)
                                                              ttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                       >  1:4099146/1‑42 (MQ=255)
                                                                 aCGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                           <  1:5685646/35‑1 (MQ=255)
                                                                 aCGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                           <  1:4524394/35‑1 (MQ=255)
                                                                 aCGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                           <  1:3068489/35‑1 (MQ=255)
                                                                 aCGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAg       <  1:1618606/71‑1 (MQ=255)
                                                                   gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGtt     >  1:1724388/1‑71 (MQ=255)
                                                                   gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGtt     <  1:1909821/71‑1 (MQ=255)
                                                                   gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGtt     <  1:91683/71‑1 (MQ=255)
                                                                     aCGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGc          >  1:4850607/1‑64 (MQ=255)
                                                                      cGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAg  <  1:84959/71‑1 (MQ=255)
                                                                      |                                                                      
TAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAG  >  minE/311810‑311950

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: