Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,830,389 A→G 12.8% A282A (GCT→GCC pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,830,3890AG12.8% 116.9 / 8.4 47A282A (GCT→GCCpvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (3/3);  total (23/24)
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.45e-01

GGCCAGGCCATAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCAC  >  NC_002947/4830242‑4830529
                                                                                                                                                   |                                                                                                                                            
ggccaggccaTAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGc                                                                                                                                            >  1:1000207/1‑150 (MQ=255)
   caggccaTAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTc                                                                                                                                         <  1:1015876/150‑1 (MQ=255)
      gccaTAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGACCGGACAGCGAAgagtcgtc                                                                                                                                      >  2:508381/1‑143 (MQ=255)
           aaCTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGAccc                                                                                                                                 <  2:1000207/150‑1 (MQ=255)
            aCTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGTCAGCGAGGCTTCGCTggccgg                                                                                                                                <  1:89825/150‑6 (MQ=255)
                  gaagaaGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCcgctcgct                                                                                                                          >  1:229064/1‑150 (MQ=255)
                  gaagaaGCGCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCcgctcgct                                                                                                                          <  1:994606/150‑1 (MQ=255)
                          tctGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACCGGc                                                                                                                  >  2:497638/1‑150 (MQ=255)
                          tctGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCCGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGc                                                                                                                  <  2:329416/150‑1 (MQ=255)
                            tGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGa                                                                                                                >  2:7446/1‑150 (MQ=255)
                                    gAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGaa                                                                                                        <  2:229064/150‑1 (MQ=255)
                                      aGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGt                                                                                                      <  1:7446/150‑1 (MQ=255)
                                         cacaGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTc                                                                                                     >  1:244187/1‑148 (MQ=255)
                                         cacaGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTc                                                                                                     <  2:244187/148‑1 (MQ=255)
                                          acaGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCggg                                                                                                  <  2:890080/150‑1 (MQ=255)
                                                 tGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATa                                                                                                              >  1:175981/1‑131 (MQ=255)
                                                 tGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATa                                                                                                              <  2:175981/131‑1 (MQ=255)
                                                    aaCTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTg                                                                                        >  1:375087/1‑150 (MQ=255)
                                                             gCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGccc                                                                               >  1:379557/1‑149 (MQ=255)
                                                                                ccTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCt                                                            <  2:817869/150‑1 (MQ=255)
                                                                                  tGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGc                                                          >  2:240796/1‑150 (MQ=255)
                                                                                       gATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCa                                                     <  2:375087/150‑1 (MQ=255)
                                                                                       gATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCa                                                     <  1:240796/150‑1 (MQ=255)
                                                                                         tAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAgc                                                   <  2:861462/150‑1 (MQ=255)
                                                                                          aGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAgcg                                                  >  1:159320/1‑150 (MQ=255)
                                                                                                   gCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCg                                         >  1:973172/1‑150 (MQ=255)
                                                                                                   gCCGGAAAACTCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCg                                         <  1:47644/150‑1 (MQ=255)
                                                                                                    ccGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCgg                                        <  2:629153/150‑1 (MQ=255)
                                                                                                       gAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGcc                                     >  2:627468/1‑150 (MQ=255)
                                                                                                           aGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGa                                 <  1:469199/150‑1 (MQ=255)
                                                                                                              ccACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGt                              >  2:658927/1‑150 (MQ=255)
                                                                                                              ccACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGt                              >  2:558767/1‑150 (MQ=255)
                                                                                                                aCGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTAc                            <  2:1063440/150‑1 (MQ=255)
                                                                                                                  gcgcCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCt                          >  1:872333/1‑150 (MQ=255)
                                                                                                                     ccAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTg                                                           >  1:613254/1‑114 (MQ=255)
                                                                                                                     ccAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTg                                                           <  2:613254/114‑1 (MQ=255)
                                                                                                                      cAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGc                      <  2:872333/150‑1 (MQ=255)
                                                                                                                      cAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGc                      >  2:602062/1‑150 (MQ=255)
                                                                                                                              cAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAg              <  1:508381/150‑1 (MQ=255)
                                                                                                                                        ccGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAAc                                                       >  1:738142/1‑99 (MQ=255)
                                                                                                                                        ccGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAAc                                                       <  2:738142/99‑1 (MQ=255)
                                                                                                                                         gggtcAGCGAGGCTTCGCTGGCCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGa                                 <  1:740043/116‑1 (MQ=255)
                                                                                                                                         gggtcAGCGAGGCTTCGCTGGCCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGa                                 >  2:740043/5‑120 (MQ=255)
                                                                                                                                         gggtcAGCGAGGCTTCGCTGGCCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCa   >  1:807183/5‑150 (MQ=255)
                                                                                                                                          ggCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCAc  >  1:631865/1‑150 (MQ=255)
                                                                                                                                           gtcAGCGAGGCTTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATg                                                                                     <  1:1032839/64‑1 (MQ=255)
                                                                                                                                           gtcAGCGAGGCTTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATg                                                                                     >  2:1032839/3‑66 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GGCCAGGCCATAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCAC  >  NC_002947/4830242‑4830529

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