Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 485,614 G→T N2K (AAC→AAA acrA ← multidrug efflux pump membrane fusion lipoprotein AcrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913485,6140GT98.4% 214.2 / ‑5.4 63N2K (AAC→AAAacrAmultidrug efflux pump membrane fusion lipoprotein AcrA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (50/12);  total (51/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TTAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTGTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAC  >  NC_000913/485553‑485674
                                                             |                                                            
ttAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTc                                                           >  1:1868299/1‑65 (MQ=255)
    ggCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTACACCCTCTGTTTTTTTTCatat                                                       >  1:2905100/1‑65 (MQ=255)
    ggCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCatat                                                       <  1:1593341/65‑1 (MQ=255)
     gCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATg                                                      >  1:1357535/1‑65 (MQ=255)
      cTAAGCTGCCTGCGAGCATCAGAACGACCGCCAGAGGCGTAAACACTCTGTTTTTTTTCATATGt                                                     >  1:41521/1‑65 (MQ=255)
      cTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGt                                                     >  1:3143030/1‑65 (MQ=255)
         cgCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTaaa                                                  <  1:2651895/64‑1 (MQ=255)
         aGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTaaa                                                  <  1:3002488/65‑1 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTaacc                                                 >  1:489238/1‑63 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:245204/1‑65 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:639921/1‑65 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:255374/1‑65 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:1214788/1‑65 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:3584816/1‑65 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:1090572/1‑65 (MQ=255)
          gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGGTTTTTTTCATATGTAAAc                                                 >  1:2054579/1‑65 (MQ=255)
           cTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAcc                                                <  1:2600777/65‑1 (MQ=255)
            tGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCt                                               >  1:998465/1‑65 (MQ=255)
                    agCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTc                                       >  1:1050088/1‑65 (MQ=255)
                    agCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTc                                       >  1:2319541/1‑65 (MQ=255)
                    agCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTGTTCATATGTAAACCTCGAGTGTc                                       >  1:951874/1‑65 (MQ=255)
                       aTCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAGCCTCGAGTGTCCgg                                    <  1:432417/65‑2 (MQ=255)
                        tCAGAACGACCGCCAGAGGAGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   <  1:2659840/65‑1 (MQ=255)
                            aaCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCa                               >  1:691992/1‑65 (MQ=255)
                             aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTGATATGTAAACCTCGAGTGTCCGATTTCaa                              >  1:2874188/1‑65 (MQ=255)
                             aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCaa                              >  1:2012348/1‑65 (MQ=255)
                             aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCaa                              >  1:3299766/1‑65 (MQ=255)
                             aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGNAAACCTCGAGTGTCCGATTTCaa                              >  1:2932825/1‑65 (MQ=255)
                                       gagGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaat                    >  1:3300944/1‑65 (MQ=255)
                                       gagGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaat                    >  1:3606797/1‑65 (MQ=255)
                                       gagGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaat                    >  1:1754639/1‑65 (MQ=255)
                                       gagGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtagtcaat                    >  1:1960941/1‑65 (MQ=255)
                                         ggCGTAAACGCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggncaatgg                  <  1:2213076/65‑1 (MQ=255)
                                         ggCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatgg                  >  1:1903559/1‑65 (MQ=255)
                                         ggCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatgg                  >  1:1231499/1‑65 (MQ=255)
                                           cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:3343306/1‑65 (MQ=255)
                                           cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:677776/1‑65 (MQ=255)
                                           cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:2884491/1‑65 (MQ=255)
                                           cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:1046956/1‑65 (MQ=255)
                                              aaaCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCaaa             <  1:1017810/65‑1 (MQ=255)
                                              aaaCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCaaa             <  1:1408263/65‑1 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:395267/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1163610/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:841219/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:733697/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:691765/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1375183/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1401549/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:494191/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:490163/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1881092/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:392192/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3622264/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:30111/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3436508/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3170424/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2602244/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2460573/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGCCAAAAg           >  1:1979827/1‑65 (MQ=255)
                                                aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGACCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2868084/1‑65 (MQ=255)
                                                aCCCTCTGTTTGTTTTCANATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2772820/1‑65 (MQ=255)
                                                 cccTCTGTTTTTTTTCATCTGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGt          >  1:2787585/1‑65 (MQ=255)
                                                 cccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGt          >  1:1823229/1‑65 (MQ=255)
                                                 cccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGt          >  1:2436981/1‑65 (MQ=255)
                                                  ccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGtt         <  1:1236260/65‑1 (MQ=255)
                                                   ttctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTta        <  1:40649/64‑1 (MQ=255)
                                                   ctctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTta        <  1:319033/65‑1 (MQ=255)
                                                   ctctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTta        <  1:2971192/65‑1 (MQ=255)
                                                   ctctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTta        <  1:808149/65‑1 (MQ=255)
                                                         ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:534387/65‑1 (MQ=255)
                                                         ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:779772/65‑1 (MQ=255)
                                                         ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:2651129/65‑1 (MQ=255)
                                                         ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:3730906/65‑1 (MQ=255)
                                                         ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  >  1:2031254/1‑65 (MQ=255)
                                                         cttttttCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:2645119/64‑1 (MQ=255)
                                                             |                                                            
TTAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTGTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAC  >  NC_000913/485553‑485674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: