Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 311,649 G→A 6.2% D408D (GAC→GAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported311,6490GA 6.2% 251.5 / 4.1 80D408D (GAC→GATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base A (2/3);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGG  >  Exported/311585‑311719
                                                                |                                                                      
cttcttTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATc                                                                   <  1:2086679/70‑1 (MQ=255)
cttcttTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATc                                                                   <  1:1499404/70‑1 (MQ=255)
    ttttGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacc                                                               >  1:1251386/1‑70 (MQ=255)
       tGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                             >  1:1358828/1‑69 (MQ=255)
         gcggcAAACCTTCTTGCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacc                                                               >  1:5323789/1‑65 (MQ=255)
         gcggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacc                                                               >  1:1743876/1‑65 (MQ=255)
         gcggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAgg                                                         >  1:951509/1‑71 (MQ=255)
         gcggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAgg                                                         >  1:4890893/1‑71 (MQ=255)
            gcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGcc                                                       >  1:1329531/1‑70 (MQ=255)
             cAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                     <  1:406510/71‑1 (MQ=255)
              aaaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCa                                                                     <  1:819344/54‑1 (MQ=255)
                    tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              >  1:109510/1‑71 (MQ=255)
                    tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              >  1:3077261/1‑71 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCATAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              >  1:5422651/1‑70 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGat                                              <  1:463952/70‑2 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              <  1:1766026/70‑1 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              <  1:1815962/70‑1 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              <  1:4820727/70‑1 (MQ=255)
                      ttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcca                                                                 >  1:4262650/1‑50 (MQ=255)
                       tcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                <  1:19042/66‑1 (MQ=255)
                       tcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                <  1:59304/66‑1 (MQ=255)
                         cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATc                                                                   >  1:2785726/1‑45 (MQ=255)
                         cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCTATGGCGAGCACCa                                         <  1:538705/71‑1 (MQ=255)
                          gCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCAt                                        <  1:1433740/71‑1 (MQ=255)
                           ccATAACACGCTCAACGTTATCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              <  1:3258236/64‑1 (MQ=255)
                             aTAACATGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                     >  1:1056061/1‑71 (MQ=255)
                             aTAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                     >  1:1644918/1‑71 (MQ=255)
                             aTAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                     >  1:4990330/1‑71 (MQ=255)
                             aTAACACGCTCAACGTTTTCTACAACAACGATGGCGTc                                                                      >  1:4867762/1‑38 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                                  >  1:1459547/1‑41 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                                  >  1:458097/1‑41 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                                  >  1:2231522/1‑41 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                                  >  1:516200/1‑41 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                    <  1:4198635/71‑1 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                    <  1:3669603/71‑1 (MQ=255)
                               aaCACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                    >  1:3579065/1‑70 (MQ=255)
                               aaCACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                    >  1:3062090/1‑70 (MQ=255)
                                acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                  >  1:5234260/1‑71 (MQ=255)
                                acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                  >  1:2018696/1‑71 (MQ=255)
                                acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                  <  1:252856/71‑1 (MQ=255)
                                  acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                >  1:3007397/1‑71 (MQ=255)
                                  acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                >  1:2951861/1‑71 (MQ=255)
                                  acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                >  1:4959289/1‑71 (MQ=255)
                                    gCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                              >  1:4596310/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                              <  1:3480826/37‑1 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                              <  1:2937573/37‑1 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                              <  1:3967913/37‑1 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                              <  1:4555398/37‑1 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                              <  1:2701007/53‑1 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                            >  1:5357395/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                            >  1:4643062/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                            >  1:4315453/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                            >  1:1579717/1‑71 (MQ=255)
                                         aCGTTTTCTACCACAACGATGGCGTCATCCACCAACAgg                                                         <  1:1445885/39‑1 (MQ=255)
                                          cGTTTTCTACCACAACGATGGCGTCATCCCCCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                        >  1:2161141/1‑71 (MQ=255)
                                          cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                   >  1:920907/1‑44 (MQ=255)
                                          cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                   >  1:2791053/1‑44 (MQ=255)
                                          cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                   >  1:2560903/1‑44 (MQ=255)
                                            ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                     <  1:1807354/40‑1 (MQ=255)
                                            ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                     <  1:4518397/40‑1 (MQ=255)
                                            ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                     <  1:1845934/40‑1 (MQ=255)
                                            ttttCTACCACAACGATGGCGTCATACACCAACAGGCCGa                                                     <  1:4354943/40‑1 (MQ=255)
                                              ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                    >  1:1397721/1‑71 (MQ=255)
                                              ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                    >  1:5433983/1‑71 (MQ=255)
                                              ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtg                     <  1:4225786/70‑1 (MQ=255)
                                               tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                    <  1:3700339/70‑1 (MQ=255)
                                                cTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                               <  1:1941279/58‑1 (MQ=255)
                                                 caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                 >  1:4149736/2‑69 (MQ=9)
                                                 caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                  <  1:4994313/69‑2 (MQ=9)
                                                   ccaccaCTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAtg               <  1:477367/71‑2 (MQ=2)
                                                   ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc               >  1:1276196/1‑71 (MQ=255)
                                                    caccacTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                  >  1:3822069/1‑51 (MQ=2)
                                                    caccacTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAtg               <  1:5562986/70‑2 (MQ=2)
                                                     acaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga             <  1:2423164/71‑1 (MQ=255)
                                                      caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCa                                            >  1:2347169/1‑39 (MQ=255)
                                                       aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                >  1:3396943/1‑50 (MQ=255)
                                                            tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa      <  1:1071670/71‑1 (MQ=255)
                                                             ggCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa      >  1:2725324/1‑70 (MQ=255)
                                                             ggCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa      >  1:958283/1‑70 (MQ=255)
                                                              gCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa      >  1:3622721/1‑69 (MQ=255)
                                                                gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg  >  1:5396240/1‑71 (MQ=255)
                                                                gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg  >  1:5049028/1‑71 (MQ=255)
                                                                gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg  <  1:2435204/71‑1 (MQ=255)
                                                                |                                                                      
CTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGG  >  Exported/311585‑311719

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: