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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported311,6460GC 6.5% 244.5 / 4.6 77A409A (GCC→GCGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (41/31);  new base C (2/3);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCA  >  Exported/311585‑311715
                                                             |                                                                     
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:4820727/70‑1 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                          <  1:1766026/70‑1 (MQ=255)
                     cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                          <  1:1815962/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:458097/1‑41 (MQ=255)
                              tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                              >  1:1459547/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:252856/71‑1 (MQ=255)
                                acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                              >  1:2018696/1‑71 (MQ=255)
                                acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                              >  1:5234260/1‑71 (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:1579717/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                        >  1:5357395/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                        >  1:4643062/1‑71 (MQ=255)
                                      tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                        >  1:4315453/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)
                                            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)
                                                             |                                                                     
CTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCA  >  Exported/311585‑311715

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: