Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,782,608 A→G 8.5% L1569P (CTC→CCC)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,782,6080AG 8.5% 161.1 / 5.1 59L1569P (CTC→CCC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base G (3/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

CACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTGCGCCTCAG  >  NC_002947/4782465‑4782755
                                                                                                                                               |                                                                                                                                                   
ccccacggTCCCAGCAACACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGGGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                                <  4:29352/143‑1 (MQ=255)
 aCCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                                <  8:43694/148‑1 (MQ=255)
 aCCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                                <  6:75060/148‑1 (MQ=255)
   ccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCTCGAACCGGGTTCGGGCatctat                                                                                                                                             <  6:34852/149‑6 (MQ=255)
   ccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtg                                                                                                                                             >  7:289949/1‑149 (MQ=255)
   ccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtg                                                                                                                                             >  6:1037/1‑149 (MQ=255)
    cAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtgt                                                                                                                                            <  7:236150/149‑1 (MQ=255)
       gTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCa                                                                                                                                         >  3:66889/1‑149 (MQ=255)
          ccAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                                        >  1:178325/1‑147 (MQ=255)
             gCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcag                                                                                                                                     >  1:141460/1‑147 (MQ=255)
                    gAGGTGCGCAGCACGTCATGGCGCTCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGTTGTGCATCAGCAGCCCTTc                                                                                                                            <  8:78933/149‑1 (MQ=255)
                           gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtgt                                                                                                                                            <  7:311080/126‑1 (MQ=255)
                           gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtgt                                                                                                                                            >  8:311080/1‑126 (MQ=255)
                                 cGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGAc                                                                                                               <  8:289949/149‑1 (MQ=255)
                                        gcgcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                                     <  7:182012/104‑1 (MQ=255)
                                        gcgcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                                     >  8:182012/1‑104 (MQ=255)
                                         cgcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATa                                                                                                       <  7:335759/149‑1 (MQ=255)
                                                  aCCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                                >  2:279874/1‑99 (MQ=255)
                                                  aCCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                                                <  1:279874/99‑1 (MQ=255)
                                                    cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                                            >  5:222974/1‑149 (MQ=255)
                                                    cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                                            >  8:172876/1‑149 (MQ=255)
                                                           cAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCt                                                                                                                          >  6:322334/1‑112 (MQ=255)
                                                              gCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATcgcgcgc                                                                                  >  6:247536/1‑149 (MQ=255)
                                                               cTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCagcagc                                                                                                                                 <  8:323969/101‑1 (MQ=255)
                                                               cTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCagcagc                                                                                                                                 >  7:323969/1‑101 (MQ=255)
                                                               cTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATcgcgcgcg                                                                                 <  4:66889/149‑1 (MQ=255)
                                                                   gCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                                        <  1:157203/90‑1 (MQ=255)
                                                                   gCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                                        >  2:157203/1‑90 (MQ=255)
                                                                    cTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCAc                                                                            >  5:348471/1‑149 (MQ=255)
                                                                              cGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGAGCGGCCACCGGCAGGCGa                                                                  >  4:313760/1‑149 (MQ=255)
                                                                                  aTCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCa                                                              <  2:158983/149‑1 (MQ=255)
                                                                                   tCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAg                                                             <  6:348471/149‑1 (MQ=255)
                                                                                     aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGTGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                                           <  2:178325/149‑1 (MQ=255)
                                                                                     aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                                           <  7:320331/149‑1 (MQ=255)
                                                                                     aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                                           >  3:123407/1‑149 (MQ=255)
                                                                                          cGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTtgctgc                                                      <  4:123407/149‑1 (MQ=255)
                                                                                           gTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTgctgct                                                     <  5:247536/149‑1 (MQ=255)
                                                                                              gTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGc                                                  >  8:264567/1‑149 (MQ=255)
                                                                                                  aCGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAg                                              <  6:222974/149‑1 (MQ=255)
                                                                                                      ggTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCg                                          <  5:1037/149‑1 (MQ=255)
                                                                                                          gCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                                        >  4:285752/1‑147 (MQ=255)
                                                                                                            aaTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGt                                    <  3:313760/149‑1 (MQ=255)
                                                                                                             aTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGCTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTg                                   >  3:58407/1‑149 (MQ=255)
                                                                                                              tACTGCATCCGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTgg                                  <  6:193822/149‑1 (MQ=255)
                                                                                                              tACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTgg                                  <  7:172876/149‑1 (MQ=255)
                                                                                                                        gatagatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCtt                        >  3:43111/1‑149 (MQ=255)
                                                                                                                          accttGCTCGAACCGGGTTCGGGCATGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                                                <  4:63004/103‑1 (MQ=255)
                                                                                                                          accttGCTCGAACCGGGTTCGGGCATGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                                                >  3:63004/5‑107 (MQ=255)
                                                                                                                            cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCg                    >  8:141576/3‑149 (MQ=255)
                                                                                                                            cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCCATTTCGGCCACAGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCg                    >  1:294789/3‑149 (MQ=255)
                                                                                                                                 ccGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                                                >  6:341952/1‑100 (MQ=255)
                                                                                                                                 ccGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                                                <  5:341952/100‑1 (MQ=255)
                                                                                                                                     aCCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCa                                               >  1:130405/1‑113 (MQ=255)
                                                                                                                                     aCCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCa                                               <  2:130405/113‑1 (MQ=255)
                                                                                                                                     aCCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCtt           <  3:101495/149‑1 (MQ=255)
                                                                                                                                       ccGGTTCGAGCAAGGTGTGCATCAGCAGACCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCATGCGATCCATTTTCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTgc         >  2:138809/1‑149 (MQ=255)
                                                                                                                                          gTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                                        <  6:103604/115‑1 (MQ=255)
                                                                                                                                          gTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                                        >  5:103604/1‑115 (MQ=255)
                                                                                                                                          gTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                                        <  4:153744/115‑1 (MQ=255)
                                                                                                                                          gTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                                        >  3:153744/1‑115 (MQ=255)
                                                                                                                                              gAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTGCGCCTCAg  <  4:58407/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTGCGCCTCAG  >  NC_002947/4782465‑4782755

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