Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 311,760 G→C 11.1% A371A (GCC→GCG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported311,7600GC11.1% 265.8 / 17.5 90A371A (GCC→GCGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base C (6/4);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATG  >  Exported/311692‑311829
                                                                    |                                                                     
gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg                                                                      >  1:1724965/1‑70 (MQ=255)
  tAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa                                                                   >  1:196454/1‑71 (MQ=255)
        gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCg                                                             >  1:3594483/1‑71 (MQ=255)
         tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCgg                                                            <  1:494426/71‑1 (MQ=255)
            aTCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt                                                         >  1:2996679/1‑71 (MQ=255)
             tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATc                                                        <  1:2303522/71‑1 (MQ=255)
                agaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa                                                       <  1:208985/69‑1 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCATCGGTACGGCaa                                                                   >  1:893220/1‑56 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                    <  1:113161/71‑1 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                    <  1:2329970/71‑1 (MQ=255)
                  aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAAc                                                     >  1:2228641/1‑69 (MQ=255)
                  aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGt                                                   >  1:2290489/1‑71 (MQ=255)
                  aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAAACCACCGGTACGGCAATGGTCGGAATCaa                                                      >  1:2715058/1‑68 (MQ=255)
                      cAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcg                                               >  1:2071402/1‑71 (MQ=255)
                         aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa                                                                   <  1:3393648/48‑1 (MQ=255)
                         aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa                                                                   <  1:3895571/48‑1 (MQ=255)
                         aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg                                             >  1:4846186/1‑70 (MQ=255)
                         aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                            <  1:3663923/71‑1 (MQ=255)
                          ggcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGa                                           <  1:2655864/71‑1 (MQ=255)
                           gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATg                                                                 >  1:2210611/1‑48 (MQ=255)
                           gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATg                                                                 >  1:778627/1‑48 (MQ=255)
                           gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCACCGTCGCGCgg                                            >  1:3497152/1‑69 (MQ=255)
                           gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                            >  1:991608/1‑69 (MQ=255)
                               cAAGGACGGCAAAGGGCCCGAGCAATACCACCGGTACGGCAATGGTCGGa                                                           <  1:4347077/50‑1 (MQ=255)
                                aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                     <  1:1388741/71‑1 (MQ=255)
                                 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                    >  1:2082518/1‑71 (MQ=255)
                                 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                    <  1:1954891/71‑1 (MQ=255)
                                 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                    >  1:2930044/1‑71 (MQ=255)
                                 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGGCGGAATCAACGTCGCGCGGAAGTTCTg                                    <  1:1104062/71‑1 (MQ=255)
                                 aGGACGGCAAAGGTCCCGAGAAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                    >  1:3748501/1‑71 (MQ=255)
                                    aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                 >  1:612465/1‑71 (MQ=255)
                                    aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                 >  1:3481233/1‑71 (MQ=255)
                                      ggCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                               <  1:3631523/71‑1 (MQ=255)
                                        cAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                             <  1:5176462/71‑1 (MQ=255)
                                          aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                           >  1:2637157/1‑71 (MQ=11)
                                          aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                           >  1:5237249/1‑71 (MQ=11)
                                          aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                           >  1:5215938/1‑71 (MQ=11)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                         <  1:2961535/56‑1 (MQ=255)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                         <  1:2566969/56‑1 (MQ=255)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                         <  1:4773430/56‑1 (MQ=255)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                         <  1:4473877/56‑1 (MQ=255)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                         <  1:13541/56‑1 (MQ=255)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                         <  1:2496179/56‑1 (MQ=255)
                                            ggTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGa                          >  1:1236659/1‑70 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                          >  1:4441074/1‑50 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                            <  1:122985/64‑1 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                      <  1:3948426/70‑1 (MQ=255)
                                                 cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAcatacata                    <  1:128063/71‑1 (MQ=255)
                                                 cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                    <  1:2566942/71‑1 (MQ=255)
                                                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                  >  1:4922309/3‑71 (MQ=2)
                                                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGAGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                  >  1:5222363/3‑71 (MQ=2)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                              >  1:2609659/1‑59 (MQ=255)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                              >  1:2727410/1‑59 (MQ=255)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                              >  1:2780819/1‑59 (MQ=255)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                  <  1:2970212/71‑1 (MQ=255)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                  <  1:966593/71‑1 (MQ=255)
                                                    gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAcc                 <  1:1008852/71‑1 (MQ=255)
                                                    acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                  >  1:1975921/2‑70 (MQ=2)
                                                    acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                  >  1:1661494/2‑70 (MQ=2)
                                                      aacaccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg               >  1:4596896/4‑71 (MQ=2)
                                                      aaTACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                     <  1:511256/65‑1 (MQ=255)
                                                       acaccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAgg              <  1:4670260/69‑1 (MQ=9)
                                                         accacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAgg              <  1:1106224/69‑1 (MQ=9)
                                                         accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                          >  1:295208/1‑41 (MQ=255)
                                                         accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                          >  1:2786853/1‑41 (MQ=255)
                                                         accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                          >  1:822193/1‑41 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa            >  1:4568543/1‑70 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                          >  1:1083101/1‑40 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                          >  1:1165736/1‑40 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                     >  1:4393492/1‑61 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa            >  1:4377344/1‑70 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAac           <  1:403968/71‑1 (MQ=255)
                                                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGAAGGAACAg                           >  1:2907685/1‑55 (MQ=255)
                                                           caccGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                            >  1:5433506/1‑37 (MQ=255)
                                                           caccGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                            >  1:4014668/1‑37 (MQ=255)
                                                           caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                             <  1:3971975/52‑1 (MQ=255)
                                                           caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                             <  1:5398554/52‑1 (MQ=255)
                                                            acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg               >  1:3124164/4‑65 (MQ=2)
                                                            acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAacac         <  1:5493422/68‑1 (MQ=9)
                                                            accGGTACGGCAATGGTCGGAATCAACTTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac         <  1:4566776/71‑1 (MQ=255)
                                                            accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg              <  1:4882378/66‑1 (MQ=255)
                                                            accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac         <  1:3584368/71‑1 (MQ=255)
                                                            accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac         <  1:297754/71‑1 (MQ=255)
                                                            accGGCACCGCGATCGTCGGGATCAGCGTCGCCCGCCAATTCTGCAGGAACAGGAACATCACCAGGAagac         <  1:4785712/71‑1 (MQ=2)
                                                             ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAg                           >  1:266931/1‑52 (MQ=255)
                                                             ccGGTACGGCAATGGTAGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                            >  1:4796944/1‑51 (MQ=255)
                                                               ggTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg      >  1:1448485/1‑71 (MQ=255)
                                                                 tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGga    >  1:2878751/1‑71 (MQ=255)
                                                                  aCGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgat   >  1:4248622/1‑71 (MQ=255)
                                                                   cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg  <  1:4147432/71‑1 (MQ=255)
                                                                   cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg  <  1:109107/71‑1 (MQ=255)
                                                                    |                                                                     
GTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATG  >  Exported/311692‑311829

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