Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,770,559 G→T 10.3% T2918K (ACG→AAG)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,770,5590GT10.3% 107.9 / 4.0 39T2918K (ACG→AAG) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base T (2/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CCAGGTAGAGTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGCCG  >  NC_002947/4770419‑4770699
                                                                                                                                            |                                                                                                                                            
ccAGGTAGAGTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTaa                                                                                                                                     >  1:755705/1‑150 (MQ=255)
 cAGGTAGAGTTCCCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAAc                                                                                                                                    <  1:272117/150‑1 (MQ=255)
     tAGAGTTCCCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACg                                                                                                                                <  1:457455/150‑1 (MQ=255)
      agagTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACgg                                                                                                                               <  1:286362/150‑1 (MQ=255)
         gTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGttt                                                                                                                            >  1:612103/1‑150 (MQ=255)
                      acCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATa                                                                                                               >  1:152526/1‑150 (MQ=255)
                           gcggcAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCACCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCAt                                                                                                          <  1:292023/150‑1 (MQ=255)
                           gcggcAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCAt                                                                                                          <  2:363874/150‑1 (MQ=255)
                                  tGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTgatgat                                                                                                   <  2:152526/150‑1 (MQ=255)
                                         aGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTtcg                                                                                                                          <  2:431371/120‑1 (MQ=255)
                                         aGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTtcg                                                                                                                          >  1:431371/1‑120 (MQ=255)
                                            tCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGc                                                                                         >  1:1121786/1‑150 (MQ=255)
                                              aaCGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGTTGAGCCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTg                                                                                       >  1:618545/1‑150 (MQ=255)
                                                   ccATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGt                                                                                  <  1:1038640/150‑1 (MQ=255)
                                                    cATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTa                                                                                 >  1:1002770/1‑150 (MQ=255)
                                                    cATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTa                                                                                 <  2:174604/150‑1 (MQ=255)
                                                          aTAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATg                                                                           >  1:1120433/1‑150 (MQ=255)
                                                            aGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATgcg                                                                         >  1:843202/1‑150 (MQ=255)
                                                                   aCGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTc                                                                  <  1:92198/150‑1 (MQ=255)
                                                                     gCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCAt                                                                >  1:1119134/1‑150 (MQ=255)
                                                                      ccgCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTGGAAAGAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATa                                                               <  1:248584/148‑1 (MQ=255)
                                                                          ccTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTg                                                           >  1:452383/1‑150 (MQ=255)
                                                                                aaCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGccc                                                                                                                  >  1:1109709/1‑89 (MQ=255)
                                                                                aaCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGccc                                                                                                                  <  2:1109709/89‑1 (MQ=255)
                                                                                aaCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGc                                                     >  2:463633/1‑150 (MQ=255)
                                                                                             gcgcATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCAc                                        <  2:755705/150‑1 (MQ=255)
                                                                                             gcgcATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCAc                                        <  2:1119134/150‑1 (MQ=255)
                                                                                                    gACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCAAGATCGCCGGGGCCACCACGATTTCGGTCGGCCCGTAGCCATTGATGATGTACTGCGGCTGCGAGTAGGCATGCCCCTGTGCATAGCTGGCCTGCGGCACCGCATCCCCGCCCaa                                 >  1:615876/1‑150 (MQ=255)
                                                                                                     aCCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAg                                <  2:1002770/150‑1 (MQ=255)
                                                                                                        gcgcCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCaa                             >  1:1032279/1‑150 (MQ=255)
                                                                                                        gcgcCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCaa                             >  1:1033753/1‑150 (MQ=255)
                                                                                                             gcAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGAt                        <  2:1120433/150‑1 (MQ=255)
                                                                                                                 ggCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGc                                                                                      <  1:751165/84‑1 (MQ=255)
                                                                                                                 ggCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGc                                                                                      >  2:751165/1‑84 (MQ=255)
                                                                                                                 ggCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATgcgc                    <  1:463633/150‑1 (MQ=255)
                                                                                                                             gcgGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCATCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTggg        >  2:1044544/1‑150 (MQ=255)
                                                                                                                                  tttCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGcc   <  2:1033753/150‑1 (MQ=255)
                                                                                                                                  tttCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGcc   <  2:1032279/150‑1 (MQ=255)
                                                                                                                                   ttCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCGCCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGCCg  >  2:987584/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                            
CCAGGTAGAGTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGCCG  >  NC_002947/4770419‑4770699

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