Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 310,779 C→G 6.8% A698A (GCG→GCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported310,7790CG 6.8% 214.7 / 4.7 73A698A (GCG→GCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (28/40);  new base G (3/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTG  >  Exported/310709‑310843
                                                                      |                                                                
ccAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                  >  1:5691961/1‑71 (MQ=255)
      cGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                              >  1:4524734/1‑69 (MQ=255)
      cGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                              >  1:5275252/1‑69 (MQ=255)
      cGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                            <  1:4411969/71‑1 (MQ=255)
          tGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                          >  1:1446769/1‑69 (MQ=255)
              cGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                    <  1:3711164/71‑1 (MQ=255)
              cGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                    <  1:4757332/71‑1 (MQ=255)
               gTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                   <  1:2909431/71‑1 (MQ=255)
               gTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                   <  1:5607459/71‑1 (MQ=255)
                  cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                   <  1:4161674/68‑1 (MQ=255)
                  cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                >  1:3586773/1‑71 (MQ=255)
                  cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                >  1:3162654/1‑71 (MQ=255)
                   gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                               <  1:1106219/71‑1 (MQ=255)
                   gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                               <  1:4482940/71‑1 (MQ=255)
                   gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                               <  1:5023481/71‑1 (MQ=255)
                       gTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           <  1:295943/71‑1 (MQ=255)
                         cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCGTGAGTCAGTTTTTCGTGAcc                                           >  1:280878/1‑69 (MQ=255)
                         cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           <  1:3969444/69‑1 (MQ=255)
                         cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           >  1:4977237/1‑69 (MQ=255)
                          aaCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                              <  1:2008415/49‑1 (MQ=255)
                          aaCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                              <  1:2260947/49‑1 (MQ=255)
                           aCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCCCCt                                                              <  1:2858905/48‑1 (MQ=255)
                            cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                      <  1:83119/71‑1 (MQ=255)
                            cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                      <  1:5482814/71‑1 (MQ=255)
                              tatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                    <  1:194227/71‑1 (MQ=255)
                                  aGGGTGCTTCGCTGCTTCTGCAAGCAACTTGTTACGCGCCTGa                                                            <  1:3392592/43‑1 (MQ=255)
                                     gTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           <  1:5345261/57‑1 (MQ=255)
                                      tGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTc                                                         >  1:160094/1‑42 (MQ=255)
                                      tGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTc                                                         >  1:3843480/1‑42 (MQ=255)
                                        cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                          >  1:5678351/1‑71 (MQ=255)
                                        cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                          <  1:4248127/71‑1 (MQ=255)
                                        cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                           <  1:2166150/70‑1 (MQ=255)
                                        cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                           <  1:4569373/70‑1 (MQ=255)
                                          tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                    <  1:3167609/43‑1 (MQ=255)
                                          tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                    <  1:3740370/43‑1 (MQ=255)
                                          tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                    <  1:1571864/43‑1 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         >  1:2705732/1‑69 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         >  1:1054075/1‑69 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         >  1:1785531/1‑69 (MQ=255)
                                             ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                 <  1:4555066/43‑1 (MQ=255)
                                              tgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                    >  1:3442614/1‑71 (MQ=255)
                                               gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                               >  1:1755266/1‑43 (MQ=255)
                                               gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                    >  1:1410835/1‑70 (MQ=255)
                                                 ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                  >  1:4169904/1‑70 (MQ=255)
                                                 ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                  >  1:2339927/1‑70 (MQ=255)
                                                  tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                 >  1:1967700/1‑54 (MQ=255)
                                                   cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                            <  1:1587239/42‑1 (MQ=255)
                                                   cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                <  1:1220595/70‑1 (MQ=255)
                                                   cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                 <  1:817166/69‑1 (MQ=255)
                                                    tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                >  1:3649713/1‑37 (MQ=255)
                                                    tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                >  1:1733112/1‑37 (MQ=255)
                                                    tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           <  1:1351816/42‑1 (MQ=255)
                                                    tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           <  1:70925/42‑1 (MQ=255)
                                                    tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                           <  1:5243633/42‑1 (MQ=255)
                                                     gcaagcaaCTTGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                          >  1:1807289/1‑42 (MQ=255)
                                                     gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                          >  1:2489477/1‑42 (MQ=255)
                                                     ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc             <  1:1626494/70‑1 (MQ=255)
                                                       aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga           <  1:999838/71‑3 (MQ=255)
                                                       aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa           <  1:1283897/71‑1 (MQ=255)
                                                       aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa           <  1:2247475/71‑1 (MQ=255)
                                                        agcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga           >  1:1523076/1‑68 (MQ=255)
                                                         gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa         >  1:4088067/1‑70 (MQ=255)
                                                         gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa         >  1:1711755/1‑70 (MQ=255)
                                                         gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag         >  1:425769/1‑71 (MQ=255)
                                                          caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc        >  1:2708401/1‑71 (MQ=255)
                                                          caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc        >  1:4575506/1‑71 (MQ=255)
                                                           aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc       <  1:5341281/71‑1 (MQ=255)
                                                             cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                               <  1:5583754/45‑1 (MQ=255)
                                                             cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         >  1:436058/1‑51 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                 <  1:2018286/56‑1 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                 <  1:835662/56‑1 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg  <  1:4797019/71‑1 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg  <  1:4700272/71‑1 (MQ=255)
                                                                  tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg  >  1:4778008/1‑69 (MQ=255)
                                                                      |                                                                
CCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTG  >  Exported/310709‑310843

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: