Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 311,664 G→C 9.5% G403G (GGC→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported311,6640GC 9.5% 318.9 / 16.4 105G403G (GGC→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (52/43);  new base C (3/7);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 110 total aligned reads displayed.

GCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCA  >  Exported/311594‑311733
                                                                      |                                                                     
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)
           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:1815962/70‑1 (MQ=255)
            cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                            <  1:1766026/70‑1 (MQ=255)
            cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                            <  1:4820727/70‑1 (MQ=255)
              tcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                              <  1:19042/66‑1 (MQ=255)
              tcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                              <  1:59304/66‑1 (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:4990330/1‑71 (MQ=255)
                    aTAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                                   >  1:1644918/1‑71 (MQ=255)
                     tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                                  <  1:3669603/71‑1 (MQ=255)
                     tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                                  <  1:4198635/71‑1 (MQ=255)
                      aaCACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                                  >  1:3062090/1‑70 (MQ=255)
                      aaCACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                                  >  1:3579065/1‑70 (MQ=255)
                       acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                <  1:252856/71‑1 (MQ=255)
                       acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                >  1:5234260/1‑71 (MQ=255)
                       acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                >  1:2018696/1‑71 (MQ=255)
                         acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                              >  1:4959289/1‑71 (MQ=255)
                         acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                              >  1:2951861/1‑71 (MQ=255)
                         acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                              >  1:3007397/1‑71 (MQ=255)
                           gCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                                            >  1:4596310/1‑71 (MQ=255)
                             tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                                          >  1:4315453/1‑71 (MQ=255)
                             tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                                          >  1:4643062/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:2560903/1‑44 (MQ=255)
                                 cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                                 >  1:2791053/1‑44 (MQ=255)
                                   ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                                   <  1:4518397/40‑1 (MQ=255)
                                   ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                                   <  1:1845934/40‑1 (MQ=255)
                                   ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGa                                                                   <  1:1807354/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)
                                            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)
                                                     gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                   >  1:1607906/4‑70 (MQ=9)
                                                     gCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                    >  1:3622721/1‑69 (MQ=255)
                                                       gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                >  1:5396240/1‑71 (MQ=255)
                                                       gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                <  1:2435204/71‑1 (MQ=255)
                                                       gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                >  1:5049028/1‑71 (MQ=255)
                                                        tCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                                          <  1:1496414/44‑1 (MQ=9)
                                                        tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc               <  1:1420939/71‑1 (MQ=255)
                                                          aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg             >  1:539611/1‑71 (MQ=255)
                                                           tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                   <  1:2642586/64‑1 (MQ=255)
                                                           tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                   <  1:5550556/64‑1 (MQ=255)
                                                           tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc            <  1:5510265/71‑1 (MQ=255)
                                                           tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg             >  1:4525407/1‑70 (MQ=255)
                                                           tCCACCAACAGGCAGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                   <  1:4851666/64‑1 (MQ=255)
                                                             caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGtt                                         >  1:4206916/1‑40 (MQ=255)
                                                             caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTGTCGAGAAGCCAAAg                 >  1:4383075/1‑64 (MQ=255)
                                                             caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                              >  1:2497605/1‑51 (MQ=255)
                                                             caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                           >  1:2975313/1‑54 (MQ=255)
                                                             caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc            >  1:65230/1‑69 (MQ=255)
                                                             caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                 >  1:2174257/1‑64 (MQ=255)
                                                              accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg         <  1:2964428/71‑1 (MQ=255)
                                                              accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg         <  1:3755918/71‑1 (MQ=255)
                                                              accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg         <  1:3421285/71‑1 (MQ=255)
                                                               ccaACAGGCCGATGGCGAGCACCATCCCGAACATTg                                           <  1:3975930/36‑1 (MQ=255)
                                                               ccaACAGGCCGATGGCGAGCACCATCCCGAACATTg                                           <  1:5381640/36‑1 (MQ=255)
                                                                gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagc               <  1:600496/60‑1 (MQ=9)
                                                                caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                              <  1:5616912/48‑1 (MQ=255)
                                                                caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                              <  1:540908/48‑1 (MQ=255)
                                                                caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                              <  1:4071180/48‑1 (MQ=255)
                                                                caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                           <  1:510460/51‑1 (MQ=255)
                                                                caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa       >  1:2970675/1‑71 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                      >  1:2718452/1‑55 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                  >  1:497563/1‑59 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                  >  1:1893672/1‑59 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                  >  1:793641/1‑59 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                    >  1:560383/1‑57 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                    >  1:5318445/1‑57 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc      >  1:5483622/1‑71 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc      >  1:469348/1‑71 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc      >  1:3244469/1‑71 (MQ=255)
                                                                 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc      >  1:4240022/1‑71 (MQ=255)
                                                                  gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAgtgt                                    <  1:161080/39‑1 (MQ=18)
                                                                  aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg         <  1:5028315/67‑1 (MQ=255)
                                                                    aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc   <  1:841289/71‑1 (MQ=255)
                                                                     ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa  <  1:572928/71‑1 (MQ=255)
                                                                     ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa  <  1:4191323/71‑1 (MQ=255)
                                                                      |                                                                     
GCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCA  >  Exported/311594‑311733

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: