Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,827,099 C→A 11.5% S1379I (AGT→ATT) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,827,0990CA11.5% 151.6 / 8.4 53S1379I (AGT→ATT) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (19/27);  new base A (2/4);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AGTGCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTG  >  NC_002947/4826950‑4827242
                                                                                                                                                     |                                                                                                                                               
agggCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAc                                                                                                                                                 <  2:480791/147‑1 (MQ=255)
  tGCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACtt                                                                                                                                               <  2:1263279/150‑1 (MQ=255)
   gCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTg                                                                                                                                              >  1:135487/1‑150 (MQ=255)
           gTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGc                                                                                                                                      <  1:207273/150‑1 (MQ=255)
            tGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAAGTGCCTGcagta                                                                                                                                     <  1:1135786/150‑5 (MQ=255)
            tGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAAGTGCCTGcagta                                                                                                                                     <  1:693172/150‑5 (MQ=255)
             gACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAAGTGCCTGcagtag                                                                                                                                    <  1:1203896/150‑6 (MQ=255)
                gTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGc                                                                                                                                         <  1:992799/142‑1 (MQ=255)
                gTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGc                                                                                                                                         >  2:992799/1‑142 (MQ=255)
                 ttGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCCCCTgg                                                                                                                                >  2:116032/1‑150 (MQ=255)
                      ggCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGaa                                                                                                                           <  2:97019/150‑1 (MQ=255)
                         ggTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAacac                                                                                                                        <  1:1173281/150‑1 (MQ=255)
                           tttCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCt                                                                                                                      >  2:819707/1‑150 (MQ=255)
                            ttCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTg                                                                                                                     <  1:116032/150‑1 (MQ=255)
                             tCAGTCGCGCCATACCGGTGGTGCAGTTGCACCTGCGGCAGCACCGGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGc                                                                                                                    <  2:1096283/150‑1 (MQ=255)
                               aGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCaa                                                                                                                  <  1:485815/150‑1 (MQ=255)
                                 tCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACa                                                                                                                >  1:979466/1‑150 (MQ=255)
                                 tCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACa                                                                                                                <  1:819707/150‑1 (MQ=255)
                                   gggCCATACCGGTTGTGCAGTTGCACCTGCCGCATCACGTGCAGCACGCGGTCACGCAGCGAGTCAAACAGCGTTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAgc                                                                                                              <  1:22928/150‑1 (MQ=255)
                                      ccATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCAGGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAgcggc                                                                                                           >  1:1196878/1‑150 (MQ=255)
                                              gctgtgCAGTTGCACCTGCGGCCGCACCTGCAGCACGTGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACaaa                                                                                                   <  1:72353/148‑1 (MQ=255)
                                                     aGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGc                                                                                            <  2:979466/150‑1 (MQ=255)
                                                           aCCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGc                                                                                               >  2:484120/1‑141 (MQ=255)
                                                           aCCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGc                                                                                               <  1:484120/141‑1 (MQ=255)
                                                             cTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCAc                                                                                    <  2:1020077/150‑1 (MQ=255)
                                                               gCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCAccc                                                                                  <  1:766375/150‑1 (MQ=255)
                                                                    agcaCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATaa                                                                             >  2:200407/1‑150 (MQ=255)
                                                                         cTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAac                                                                                                                          >  2:1141243/1‑100 (MQ=255)
                                                                         cTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAac                                                                                                                          <  1:1141243/100‑1 (MQ=255)
                                                                         cTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTg                                                                        <  2:1196878/150‑1 (MQ=255)
                                                                               caCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAg                                                                  <  1:781937/150‑1 (MQ=255)
                                                                                aCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAgg                                                                 <  2:959723/150‑1 (MQ=255)
                                                                                   cgGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCa                                                              <  1:200407/150‑1 (MQ=255)
                                                                                      tCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCg                                                           >  2:157654/1‑150 (MQ=255)
                                                                                                  gCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCt                                                                                                                      <  1:280554/79‑1 (MQ=255)
                                                                                                  gCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCt                                                                                                                      >  2:280554/1‑79 (MQ=255)
                                                                                                         aCGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGAt                                        <  1:157654/150‑1 (MQ=255)
                                                                                                                gCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCAcc                                 >  1:945622/1‑150 (MQ=255)
                                                                                                                   aCCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCgg                                                >  2:849923/1‑132 (MQ=255)
                                                                                                                   aCCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCgg                                                <  1:849923/132‑1 (MQ=255)
                                                                                                                   aCCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAgg                              >  1:951891/1‑150 (MQ=255)
                                                                                                                     ccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCa                                                                                              >  2:956133/1‑84 (MQ=255)
                                                                                                                     ccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCa                                                                                              <  1:956133/84‑1 (MQ=255)
                                                                                                                          aaCACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCg                       >  1:26686/1‑150 (MQ=255)
                                                                                                                                  gcgcAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAg               >  1:583570/1‑150 (MQ=255)
                                                                                                                                     cAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCtt            <  2:26686/150‑1 (MQ=255)
                                                                                                                                       ccgttatGCCGGCAGTTTCATGCGCCTCCTTGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGc          <  2:915601/144‑1 (MQ=255)
                                                                                                                                         gagccgCAGGCAAGTGCCTGCAGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAg        >  1:180611/6‑150 (MQ=255)
                                                                                                                                            ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGa                                         >  2:193269/3‑114 (MQ=255)
                                                                                                                                            ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGa                                         <  1:193269/112‑1 (MQ=255)
                                                                                                                                              gCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGt   <  2:135487/150‑1 (MQ=255)
                                                                                                                                               cAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAgg                              <  1:996598/122‑1 (MQ=255)
                                                                                                                                               cAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAgg                              >  2:996598/1‑122 (MQ=255)
                                                                                                                                               cAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCaagatcgg  >  1:915601/1‑143 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AGTGCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTG  >  NC_002947/4826950‑4827242

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