Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,558 A→C 16.4% V105V (GTT→GTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5580AC16.4% 225.5 / 52.0 116V105V (GTT→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (58/39);  new base C (9/10);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 120 total aligned reads displayed.

ACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312493‑312627
                                                                 |                                                                     
  cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                 <  1:9975/70‑1 (MQ=255)
   ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               >  1:2540559/1‑71 (MQ=255)
   ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               <  1:3724069/71‑1 (MQ=255)
      tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa                                                             >  1:2814008/1‑70 (MQ=255)
      tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa                                                             >  1:69298/1‑70 (MQ=255)
      tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa                                                             >  1:114107/1‑70 (MQ=255)
       gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa                                                             <  1:923469/69‑1 (MQ=255)
        ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                                          <  1:2940871/71‑1 (MQ=255)
        ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                                          <  1:1885283/71‑1 (MQ=255)
        ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                                           <  1:3637246/70‑1 (MQ=255)
        ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCAGCACCTGCACTTGCGCGatat                                                          <  1:2365335/71‑1 (MQ=255)
          gctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTgcgc                                                               <  1:2523939/64‑1 (MQ=255)
           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                  >  1:3640926/1‑60 (MQ=255)
           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                       >  1:734831/1‑71 (MQ=255)
           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                       >  1:1844797/1‑71 (MQ=255)
            tgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                            >  1:1130171/1‑65 (MQ=255)
            tgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATcag                                                       >  1:3561527/1‑68 (MQ=255)
             gAACTTCTTGCGGCAGCAACTGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                      >  1:846843/1‑70 (MQ=255)
             gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCa                                                     >  1:265643/1‑71 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     <  1:9925/54‑1 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                      <  1:859700/69‑1 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                      <  1:854367/69‑1 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                      <  1:849499/69‑1 (MQ=255)
              aaCCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATcagg                                                      <  1:598713/69‑4 (MQ=19)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:3351000/1‑70 (MQ=255)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:3482184/1‑70 (MQ=255)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:2828559/1‑70 (MQ=255)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATc                                                   <  1:1409466/71‑1 (MQ=255)
                cttcttGCGGCAGCAACGGTATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:2097907/1‑69 (MQ=255)
                cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCt                                                                   >  1:2098160/1‑54 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                 <  1:1706671/69‑1 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                <  1:2299812/70‑1 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                >  1:2067582/1‑70 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                               >  1:10732/1‑71 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTTCACTTGCGCGATATCAGGATCGGt                                                <  1:94435/68‑1 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                              >  1:2054891/2‑69 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                              >  1:1607135/2‑69 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                               <  1:1321715/69‑2 (MQ=16)
                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:1500614/1‑69 (MQ=255)
                       cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccgg                                           <  1:2303622/71‑3 (MQ=255)
                       cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                            <  1:1935227/70‑1 (MQ=255)
                       cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                           <  1:497714/71‑1 (MQ=255)
                       cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                           <  1:3529200/71‑1 (MQ=255)
                        ggCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                          >  1:484397/1‑71 (MQ=16)
                         gcagcaTCGGCATCGCCAGCTGCAGTTTGTTGTGCACCTGAACCTGCGCGATATCCGCa                                                     <  1:1993315/59‑1 (MQ=255)
                         gcagcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGAc                                         <  1:352053/71‑1 (MQ=16)
                         gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                         <  1:2401338/71‑1 (MQ=255)
                         gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                         <  1:503449/71‑1 (MQ=255)
                          cagcaACGGCGTGGCGAGCTGGAGTTTGTTctgcacctgcac                                                                     >  1:2626849/1‑42 (MQ=16)
                          cagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                         >  1:315622/1‑70 (MQ=255)
                          cagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                        <  1:2520448/71‑1 (MQ=255)
                           agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATc                                                   >  1:3313306/1‑59 (MQ=255)
                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTc                                       >  1:1205334/1‑70 (MQ=255)
                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:3253030/71‑1 (MQ=255)
                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:684315/1‑71 (MQ=255)
                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:689156/1‑71 (MQ=255)
                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:962620/1‑71 (MQ=255)
                                cGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAgg                                  >  1:605973/1‑71 (MQ=255)
                                cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                  >  1:1005238/1‑71 (MQ=255)
                                cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAgg                                  >  1:460897/1‑71 (MQ=255)
                                cGGCATCGCCAGCTGCAGTTTGTTATGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                  >  1:3402512/1‑71 (MQ=255)
                                 ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 >  1:651667/1‑71 (MQ=255)
                                   cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                               <  1:3712403/71‑1 (MQ=255)
                                      cGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                        >  1:598168/1‑59 (MQ=255)
                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaaa                                    >  1:3284913/1‑62 (MQ=255)
                                        cgaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGta                          <  1:2853999/69‑2 (MQ=12)
                                        ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:1102180/71‑1 (MQ=255)
                                          aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        <  1:2366295/71‑1 (MQ=255)
                                            cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                     <  1:3310785/56‑1 (MQ=255)
                                              gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:1511077/1‑56 (MQ=255)
                                              gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                     >  1:633201/1‑70 (MQ=255)
                                              gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAAACTCAAAGGTCAGGGTGATCTg                     >  1:2604417/1‑70 (MQ=255)
                                               cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                   >  1:1096764/1‑71 (MQ=255)
                                                aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                               >  1:2509207/1‑42 (MQ=255)
                                                aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                  >  1:1255720/1‑71 (MQ=255)
                                                aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                  >  1:1640108/1‑71 (MQ=255)
                                                  tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:949551/54‑1 (MQ=255)
                                                  tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:2747186/54‑1 (MQ=255)
                                                   ttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                 >  1:2051342/1‑37 (MQ=255)
                                                    tGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                  >  1:1010098/1‑35 (MQ=255)
                                                    tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:1017383/1‑71 (MQ=255)
                                                     gTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              <  1:188869/70‑1 (MQ=255)
                                                       tCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAGTCAAAGGTCAggg                            <  1:3450658/54‑1 (MQ=255)
                                                        cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGc                    >  1:1241334/1‑61 (MQ=255)
                                                        cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           >  1:3303864/1‑70 (MQ=255)
                                                        cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg          >  1:541457/1‑71 (MQ=255)
                                                         tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt         >  1:1424638/1‑71 (MQ=255)
                                                          gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg          <  1:2500555/69‑1 (MQ=255)
                                                          gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg          <  1:2697790/69‑1 (MQ=255)
                                                            aCCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                         <  1:3729467/36‑1 (MQ=255)
                                                            aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        >  1:63853/1‑69 (MQ=255)
                                                            aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        >  1:1255845/1‑69 (MQ=255)
                                                              cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:1536096/37‑1 (MQ=255)
                                                              cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:1609936/1‑61 (MQ=255)
                                                              cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt    >  1:3256021/1‑71 (MQ=255)
                                                               tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt    >  1:1926290/1‑70 (MQ=255)
                                                                gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc   >  1:2730804/1‑70 (MQ=255)
                                                                gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa  >  1:1004204/1‑71 (MQ=255)
                                                                 aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc   <  1:1145281/69‑1 (MQ=255)
                                                                 aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCACGGTGATCTGCACGGTACCCGTGGAGTCa  <  1:2806951/70‑1 (MQ=255)
                                                                 |                                                                     
ACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312493‑312627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: