Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,779,409 T→G 8.5% intergenic (‑143/+18) pvdD ← / ← pvdJ non‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,779,4090TG 8.5% 135.9 / 4.3 47intergenic (‑143/+18)pvdD/pvdJnon‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (17/26);  new base G (2/2);  total (19/28)
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.39e-01

GCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAGG  >  NC_002947/4779263‑4779557
                                                                                                                                                  |                                                                                                                                                    
gCACTTGAACCTCTCCTACCAGGGCACCGGTCTTGTCCTGCGGAACGGGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAt                                                                                                                                                      <  6:262446/147‑1 (MQ=255)
  aCTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCg                                                                                                                                                  <  7:198120/149‑1 (MQ=255)
    ttGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCggg                                                                                                                                                >  5:147468/1‑149 (MQ=255)
        aCATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCg                                                                                                                                            <  5:117353/149‑1 (MQ=255)
        aCATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCg                                                                                                                                            <  5:117350/149‑1 (MQ=255)
           tctcCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCg                                                                                                                                            >  7:338597/1‑146 (MQ=255)
                   cAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTAc                                                                                                                                 >  7:194204/1‑149 (MQ=255)
                   cAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTAc                                                                                                                                 >  7:199474/1‑149 (MQ=255)
                   cAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTAc                                                                                                                                 <  7:39907/149‑1 (MQ=255)
                      ggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTgg                                                                                                                              <  2:41521/149‑1 (MQ=255)
                                     cTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtc                                                                                                               >  1:9492/1‑149 (MQ=255)
                                            aCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGTGCGCTACTGGGTCAACTCGTCGATCTCctct                                                                                                         <  4:321878/148‑1 (MQ=255)
                                             cGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGc                                                                                                       <  1:239896/149‑1 (MQ=255)
                                              gTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCa                                                                                                      <  8:131932/149‑1 (MQ=255)
                                              gTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCa                                                                                                      <  2:9492/149‑1 (MQ=255)
                                                 ccAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTa                                                                                                   <  8:194204/149‑1 (MQ=255)
                                                 ccAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTa                                                                                                   <  8:199474/149‑1 (MQ=255)
                                                   aGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGt                                                                                                    >  8:66653/1‑146 (MQ=255)
                                                       cgcgACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGt                                                                                             <  7:15911/149‑1 (MQ=255)
                                                          gACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTc                                                                                          <  2:340830/149‑1 (MQ=255)
                                                            cAGGCGCATGATGCTGCCCAACAATTGACAGAGCCGGACCGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCaa                                                                                        <  4:275779/149‑1 (MQ=255)
                                                               gCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggc                                                                                     <  4:274240/149‑1 (MQ=255)
                                                               gCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggc                                                                                     >  4:338399/1‑149 (MQ=255)
                                                                 tCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTAGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg                                                                                   >  1:311325/1‑149 (MQ=255)
                                                                  cATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcggc                                                                                  >  5:208551/1‑149 (MQ=255)
                                                                       tgctgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAgg                                                                             <  3:338399/149‑1 (MQ=255)
                                                                         ctgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAt                                                                           <  3:293949/149‑1 (MQ=255)
                                                                            ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTg                                                                                                                               >  7:32697/1‑94 (MQ=255)
                                                                            ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTg                                                                                                                               <  8:32697/94‑1 (MQ=255)
                                                                             cGAACAAATGAAAGACCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCGGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAg                                                                       <  8:215055/149‑1 (MQ=255)
                                                                                      gAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGt                                                              >  5:335132/1‑149 (MQ=255)
                                                                                         cgagCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTg                                                                                                        <  3:8654/103‑1 (MQ=255)
                                                                                         agagCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCt                                                           >  4:77657/1‑149 (MQ=255)
                                                                                                gACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGa                                                    <  7:66653/149‑1 (MQ=255)
                                                                                                gACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGa                                                    >  3:92892/1‑149 (MQ=255)
                                                                                                           gCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgcc                                         <  7:316789/149‑1 (MQ=255)
                                                                                                                 ggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTg                                     <  3:77657/147‑1 (MQ=255)
                                                                                                                 ggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGc                                    <  6:335132/148‑1 (MQ=255)
                                                                                                                       cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                              <  6:282410/148‑1 (MQ=255)
                                                                                                                       cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                              >  5:282410/1‑148 (MQ=255)
                                                                                                                        tGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGc                            >  6:121808/1‑149 (MQ=255)
                                                                                                                            gTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACg                        >  8:56746/1‑149 (MQ=255)
                                                                                                                             ttGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGAAGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACgg                       >  1:15043/1‑149 (MQ=255)
                                                                                                                                    tgcGGGCCCGGTATGACCGGGCCCACAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCt                >  8:6053/3‑149 (MQ=255)
                                                                                                                                      cGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGc                            <  6:338233/135‑1 (MQ=255)
                                                                                                                                      cGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGcc                     <  8:89085/142‑1 (MQ=255)
                                                                                                                                      cGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGcc                     >  7:89085/1‑142 (MQ=255)
                                                                                                                                            cGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTc        <  6:208551/149‑1 (MQ=255)
                                                                                                                                                  tACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGt                                                              <  8:31670/89‑1 (MQ=255)
                                                                                                                                                  tACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGt                                                              >  7:31670/1‑89 (MQ=255)
                                                                                                                                                  tACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAgg  <  4:92892/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAGG  >  NC_002947/4779263‑4779557

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