Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,824,317 G→C 11.1% S2306R (AGC→AGG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,824,3170GC11.1% 104.3 / 3.9 36S2306R (AGC→AGG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base C (2/2);  total (18/18)
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.93e-01

ACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGA  >  NC_002947/4824174‑4824465
                                                                                                                                               |                                                                                                                                                    
aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcg                                                                                                                                                >  2:549405/1‑150 (MQ=255)
  cGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcggc                                                                                                                                              >  2:267509/1‑150 (MQ=255)
  cGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggc                                                                                                                                              <  1:549405/150‑1 (MQ=255)
    ccTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCAGGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCAt                                                                                                                                            >  1:146740/1‑150 (MQ=255)
        gCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggc                                                                                                                                                 >  2:591551/1‑141 (MQ=255)
        gCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggc                                                                                                                                                 <  1:591551/141‑1 (MQ=255)
             gCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCa                                                                                                                                   <  2:148592/150‑1 (MQ=255)
             gCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCa                                                                                                                                   <  1:551192/150‑1 (MQ=255)
             ccATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCa                                                                                                                                   >  2:502544/2‑150 (MQ=255)
              cATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGGTACGTAGTTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAg                                                                                                                                  <  2:478885/150‑1 (MQ=255)
              cATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAg                                                                                                                                  >  2:607983/1‑150 (MQ=255)
                  gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTc                                                                                                                              <  1:406358/150‑1 (MQ=255)
                       ttCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCt                                                                                                                         <  1:502544/149‑1 (MQ=255)
                        tCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTg                                                                                                                        >  2:481594/1‑150 (MQ=255)
                         cGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTgg                                                                                                                       <  2:427279/150‑1 (MQ=255)
                          gCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGGGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGc                                                                                                                      >  1:567838/1‑150 (MQ=255)
                                   aaccaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTc                                                                                                             <  2:601655/149‑1 (MQ=255)
                                      caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAgg                                                                                                          >  1:446603/1‑150 (MQ=255)
                                      caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctacag                                                                                                                                        >  1:10043/1‑114 (MQ=255)
                                      caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctacag                                                                                                                                        <  2:10043/120‑7 (MQ=255)
                                          acGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCTGGGAGGTCGGGTTCTCCGTCAAACAAAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGc                                                                                                      >  2:26451/1‑150 (MQ=255)
                                                           ggTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCCCCGTTTCAGGCAGCTCGCCCAGGGTGTgcag                                                                                     <  1:267509/150‑1 (MQ=255)
                                                                 gCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAgcgc                                                                               >  1:179805/1‑150 (MQ=255)
                                                                         tgtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggc                                                                       <  1:481594/150‑1 (MQ=255)
                                                                           tAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggctg                                                                     >  1:184199/1‑150 (MQ=255)
                                                                                        ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgc                                                                           >  2:192153/1‑131 (MQ=255)
                                                                                        ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgc                                                                           <  1:192153/131‑1 (MQ=255)
                                                                                        ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGa                                                        <  2:202691/150‑1 (MQ=255)
                                                                                        ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGAACGCCAGGGTTACAGCATCGTCCTGGAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGa                                                        >  1:76172/1‑150 (MQ=255)
                                                                                         cTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGAt                                                       >  1:328848/1‑150 (MQ=255)
                                                                                         cTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGAt                                                       <  2:125063/150‑1 (MQ=255)
                                                                                          tGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCCGCAGCGCACGCTGGCTGAGCAGCAAGCCGATg                                                      <  2:567838/150‑1 (MQ=255)
                                                                                             tgctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCgg                                                   <  2:328848/150‑1 (MQ=255)
                                                                                                             ggTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGt                                   >  2:29765/1‑150 (MQ=255)
                                                                                                                         gTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGt                       <  2:146740/150‑1 (MQ=255)
                                                                                                                                         cgccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCggg       >  2:279868/1‑150 (MQ=255)
                                                                                                                                              ggCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGa  <  2:179805/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTCGA  >  NC_002947/4824174‑4824465

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