Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,740,210 T→G 7.5% L218R (CTG→CGG) ‡ PP_3304 → Bcr/CflA family multidrug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,740,2100TG 7.5% 162.8 / 3.7 53L218R (CTG→CGG) ‡PP_3304Bcr/CflA family multidrug resistance transporter
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base G (2/2);  total (26/27)
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.03e-01

GTATGCTGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAGCTCGCTCCAGGTGG  >  NC_002947/3740066‑3740356
                                                                                                                                                |                                                                                                                                                  
gTATGCTGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATc                                                                                                                                                <  5:115598/149‑1 (MQ=255)
    gCTGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCa                                                                                                                                             >  5:150611/1‑148 (MQ=255)
     cTGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCAGGATCGCag                                                                                                                                            <  2:262162/148‑2 (MQ=255)
     cTGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGGTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTTCCTGCGATCCTGATCGCAt                                                                                                                                            <  6:108306/148‑1 (MQ=255)
      tGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATAc                                                                                                                                          >  5:308424/1‑149 (MQ=255)
        aaCACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACgg                                                                                                                                        >  5:23839/1‑149 (MQ=255)
             cacTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCAGGATCGCagg                                                                                                                                           <  2:140407/141‑3 (MQ=255)
             cacTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCAGGATCGCagg                                                                                                                                           >  1:140407/1‑139 (MQ=255)
               cTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTg                                                                                                                                 >  7:328674/1‑149 (MQ=255)
                  tGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGc                                                                                                                               <  6:23839/148‑1 (MQ=255)
                                 ttttCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGc                                                                                                                               <  3:155792/133‑1 (MQ=255)
                                 ttttCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGc                                                                                                                               >  4:155792/1‑133 (MQ=255)
                                 ttttCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCtt                                                                                                                >  4:262662/1‑148 (MQ=255)
                                   ttCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTc                                                                                                               <  3:262662/147‑1 (MQ=255)
                                    tCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATc                                                                                                            >  7:41710/1‑149 (MQ=255)
                                      aTAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCtt                                                                                                          >  2:393/1‑149 (MQ=255)
                                                gCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGCTGACTCAATCAGTACCTGCGATCCTGAACGCATACGGCAAATTGCCGTCCAACCGCCTCTTCATCTTGCCCGCACTc                                                                                                <  5:92507/149‑1 (MQ=255)
                                                  aTGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTcgc                                                                                              >  6:204981/1‑149 (MQ=255)
                                                     tCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGTAAATTGGCGTCCAACCGCctct                                                                                                                 <  8:181909/127‑1 (MQ=255)
                                                     tCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGTAAATTGGCGTCCAACCGCctct                                                                                                                 >  7:181909/1‑127 (MQ=255)
                                                     tCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGAc                                                                                           <  4:349646/149‑1 (MQ=255)
                                                     tCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGAc                                                                                           <  6:150611/149‑1 (MQ=255)
                                                         ataTCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGc                                                                                       <  2:81692/147‑1 (MQ=255)
                                                                 aTCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTc                                                                                                >  7:231974/1‑132 (MQ=255)
                                                                 aTCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTc                                                                                                <  8:231974/132‑1 (MQ=255)
                                                                               ccTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTc                                                                  <  1:336650/148‑1 (MQ=255)
                                                                                cTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGc                                                                <  4:273370/149‑1 (MQ=255)
                                                                                cTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGc                                                                >  6:195511/1‑149 (MQ=255)
                                                                                       ggATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTc                                                         <  5:204981/149‑1 (MQ=255)
                                                                                              tGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTg                                                  >  7:82546/1‑149 (MQ=255)
                                                                                               gAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTgc                                                 >  6:370214/1‑149 (MQ=255)
                                                                                                aaaCCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTgcg                                                >  4:331273/1‑149 (MQ=255)
                                                                                                  aCCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGcc                                              >  2:208038/1‑149 (MQ=255)
                                                                                                    ccTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGt                                            >  3:228126/1‑149 (MQ=255)
                                                                                                        ccTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCaa                                         <  2:54457/148‑1 (MQ=255)
                                                                                                            cccATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCtt                                     >  8:110734/1‑148 (MQ=255)
                                                                                                            cccATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTg                                    <  3:331273/149‑1 (MQ=255)
                                                                                                            cccATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTg                                    <  1:393/149‑1 (MQ=255)
                                                                                                             ccATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGa                                   <  8:328674/149‑1 (MQ=255)
                                                                                                             ccATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCGTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGa                                   >  4:250271/1‑149 (MQ=255)
                                                                                                             ccATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGCCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGa                                   >  1:270754/1‑149 (MQ=255)
                                                                                                              cATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGAt                                  <  4:212671/149‑1 (MQ=255)
                                                                                                                  gtcgGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGAt                                  >  6:197754/1‑145 (MQ=255)
                                                                                                                  gtcgGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGAt                                  <  5:197754/145‑1 (MQ=255)
                                                                                                                     gggTGACTCAATCAGTACCTGCGATCCTGATCTCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAgg                            >  8:229119/1‑148 (MQ=255)
                                                                                                                     gggTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAgg                            <  5:195511/148‑1 (MQ=255)
                                                                                                                        tGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATg                        <  8:41710/149‑1 (MQ=255)
                                                                                                                          aCTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATgg                       <  6:308424/148‑1 (MQ=255)
                                                                                                                                  aGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAgc               <  4:228126/148‑1 (MQ=255)
                                                                                                                                    cgtatGCGATCAGGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAgctcg            >  5:93432/5‑149 (MQ=255)
                                                                                                                                         gCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAgctcgctc         >  3:204538/1‑147 (MQ=255)
                                                                                                                                          cGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTAATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAGCTCGCTCCAg      >  2:131313/1‑149 (MQ=255)
                                                                                                                                               cTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAGCTCGCTCCAGGTgg  <  7:229119/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTATGCTGAACACCACTTTGGGGTGGCGCGCACTTTTCATAGCCGTCGGCATGTCGGCAATCCTTATCGCGCTGTTCTACCTCAGAGGGATAGGTGAAACCCTTCCTGCCCATCGTCGGGTGACTCAATCAGTACCTGCGATCCTGATCGCATACGGCAAATTGGCGTCCAACCGCCTCTTCATCTTGCCCGCACTCGCGACCAGCCTGCTGATGAGCGGCTTATTCGCCTCGTTCGCCGCTGCGCCGTCAATCTTGATGGAAGGCATGGGGCTCAGCTCGCTCCAGGTGG  >  NC_002947/3740066‑3740356

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