Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,792,573 C→T 8.7% V438M (GTG→ATG)  pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,792,5730CT 8.7% 120.3 / 3.1 46V438M (GTG→ATG) pvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AGCGGCCAGTCACCCAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCGGGCTCACGCCC  >  NC_002947/4792425‑4792718
                                                                                                                                                    |                                                                                                                                                 
aGCGGCCAGTCACCCAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACg                                                                                                                                                  <  2:473796/150‑1 (MQ=255)
  cGGCCAGTCACCCAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGgt                                                                                                                                                <  1:676422/150‑1 (MQ=255)
        gTCACCCAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAg                                                                                                                                          <  2:386030/150‑1 (MQ=255)
         tCACCCAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAgg                                                                                                                                         >  1:22855/1‑150 (MQ=255)
             ccAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTaa                                                                                                                                     >  2:360118/1‑150 (MQ=255)
               agcgcagcgACGTCCCCGCGGCTAAGGGCGGCTAGTAGGCGCGGCCAACCACGCCCCATGGCCTCCACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCCGCAACCTCCACGGTGTCCAGGGTAAgg                                                                                                                                   <  2:242352/150‑1 (MQ=255)
                    agcgACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGGAGCAGGCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGaa                                                                                                                              <  1:748113/150‑1 (MQ=255)
                        aCGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGt                                                                                                                          >  1:1107334/1‑150 (MQ=255)
                             aTCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGt                                                                                                                     >  1:555842/1‑150 (MQ=255)
                               cgcgGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGc                                                                                                                   >  1:666851/1‑150 (MQ=255)
                                    cTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTg                                                                                                              >  2:40539/1‑150 (MQ=255)
                                            ggcCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTc                                                                                                      <  1:1128740/150‑1 (MQ=255)
                                               cagcagGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAg                                                                                                   >  1:710325/1‑150 (MQ=255)
                                                     gCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGTGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCAc                                                                                                                                                   <  1:1030907/96‑1 (MQ=255)
                                                     gCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGTGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCAc                                                                                                                                                   >  2:1030907/1‑96 (MQ=255)
                                                                    ccccATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGa                                                                              >  2:69348/1‑150 (MQ=255)
                                                                            ccTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTgc                                                                      >  1:39547/1‑150 (MQ=255)
                                                                            ccTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCACGAGGCGGTgc                                                                      <  1:99415/150‑1 (MQ=255)
                                                                             cTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTgcg                                                                     >  2:1117980/1‑150 (MQ=255)
                                                                                 aCGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGcaca                                                                 <  1:360118/150‑1 (MQ=255)
                                                                                    gtgtCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTg                                                              <  1:1064134/150‑1 (MQ=255)
                                                                                        cgcgCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGa                                                          <  1:1117980/150‑1 (MQ=255)
                                                                                         gcgcTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGAc                                                         <  2:989079/150‑1 (MQ=255)
                                                                                             tCGAACAGGTCAAGGGCAAAGGGGGAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGTAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCg                                                     >  2:855762/1‑150 (MQ=255)
                                                                                               gAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCATGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGt                                                   >  1:441400/1‑150 (MQ=255)
                                                                                               gAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCATGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGGCGGt                                                   >  1:440993/1‑150 (MQ=255)
                                                                                                      tCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGGGTCCAGGGTAAGGTCGCAGCGCCCACTGCCCTTGCGCACCCCGCGCTTTTCCAGCCCAATGCCCCTTGTGTAGCGCCCAGCGGCGGCCTCGGTGACCacc                                            >  2:68356/1‑148 (MQ=255)
                                                                                                        aaGGGCAAAGGTGAAGGCGGCGCGCAGCTGCCCGGCAACCACCACGGTGTCCAGGGGAACGGCGAGGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGa                                          >  2:51847/1‑150 (MQ=255)
                                                                                                             cAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCATGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGAt                                     <  2:441400/150‑1 (MQ=255)
                                                                                                             cAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCATGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGAt                                     <  2:440993/150‑1 (MQ=255)
                                                                                                             cAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAATGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGAt                                     >  1:631158/1‑150 (MQ=255)
                                                                                                              aaaGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATg                                    <  1:627736/150‑1 (MQ=255)
                                                                                                               aaGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATgg                                   <  1:1111571/150‑1 (MQ=255)
                                                                                                                  gTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTg                                                                       >  2:860743/1‑111 (MQ=255)
                                                                                                                  gTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTg                                                                       <  1:860743/111‑1 (MQ=255)
                                                                                                                      aGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGaaaa                            <  2:555842/150‑1 (MQ=255)
                                                                                                                          ggcgCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGGTGGTTGAAAAGCAc                        <  2:1107334/150‑1 (MQ=255)
                                                                                                                          ggcgCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCAc                        >  2:327415/1‑150 (MQ=255)
                                                                                                                                   cTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAg               <  1:1100436/150‑1 (MQ=255)
                                                                                                                                    tGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGc              <  1:1087543/150‑1 (MQ=255)
                                                                                                                                       ccGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCggg           <  1:593282/150‑1 (MQ=255)
                                                                                                                                         ggCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCGGGCt         <  2:22855/150‑1 (MQ=255)
                                                                                                                                            aaCCTCCACGGTGTCCATGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCGCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCGGGCTCAc      >  2:646802/1‑150 (MQ=255)
                                                                                                                                            aaCCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCGGGCTCAc      >  2:1037760/1‑150 (MQ=255)
                                                                                                                                               cTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGaaaa                            <  1:345512/125‑1 (MQ=255)
                                                                                                                                               cTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGaaaa                            >  2:345512/1‑125 (MQ=255)
                                                                                                                                                tCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCGGGCTCACGccc  >  1:177684/1‑150 (MQ=255)
                                                                                                                                                tCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCAAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAAGCGCGGGCTCACGccc  >  1:391849/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGCGGCCAGTCACCCAGCGCAGCGACGTCATCGCGGCTGAGGGCGGCCAGCAGGCTCAGCCAATCACGCCCCATGGCCTCGACGGTGTCGCGCTCGAACAGGTCAAGGGCAAAGGTGAAGGCGGCGCGAAGCTGCCCGGCAACCTCCACGGTGTCCAGGGTAAGGTCGAAGCGTGCAGTGCTCTTGCGCAGCTCGAGTTTTTCCAGCTCAAGGCCCGAGGCGGTGCGCACAGTGGCGACGTCGGTGACCATCGACTGATGGTTGAAAAGCACCTGGAACAGCGGGCTCACGCCC  >  NC_002947/4792425‑4792718

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: ‑