Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,779,407 C→A 9.1% intergenic (‑141/+20) 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,4070CA 9.1% 145.4 / 4.3 51intergenic (‑141/+20)pvdD/pvdJnon‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base A (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCC  >  NC_002947/4779260‑4779552
                                                                                                                                                   |                                                                                                                                                 
cTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGTCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCa                                                                                                                                                  >  4:27219/1‑149 (MQ=255)
         gAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCggg                                                                                                                                           <  1:199642/147‑1 (MQ=255)
         gAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATCCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCggg                                                                                                                                           <  1:11005/147‑1 (MQ=255)
             aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                     <  1:127283/149‑1 (MQ=255)
             aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                     <  5:24415/149‑1 (MQ=255)
                  cTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAgc                                                                                                                                 <  8:33643/148‑1 (MQ=255)
                      cAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTATGACCGGGCCCgg                                                                                                                                      <  4:295435/139‑2 (MQ=255)
                      cAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTATGACCGGGCCCgg                                                                                                                                      >  3:295435/1‑138 (MQ=255)
                       aGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACt                                                                                                                           <  1:10985/149‑1 (MQ=255)
                            aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGt                                                                                                                       >  1:13031/1‑148 (MQ=255)
                                 tCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtc                                                                                                                 >  5:75772/1‑149 (MQ=255)
                                          gCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCg                                                                                                                                       <  1:98200/118‑1 (MQ=255)
                                          gCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCg                                                                                                                                       >  2:98200/1‑118 (MQ=255)
                                           cGGAACGTGCCAGTCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtctc                                                                                                        <  3:27219/148‑1 (MQ=255)
                                           cGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATctcc                                                                                                       >  5:85452/1‑149 (MQ=255)
                                             gAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATctcctc                                                                                                     >  6:141101/1‑149 (MQ=255)
                                              aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCctct                                                                                                    >  2:85312/1‑149 (MQ=255)
                                               aCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTg                                                                                                   <  7:194720/149‑1 (MQ=255)
                                                         ccGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACTGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACg                                                                                         <  1:85312/149‑1 (MQ=255)
                                                                   cTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg                                                                               <  5:141101/149‑1 (MQ=255)
                                                                      atgatgCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcggc                                                                             <  6:85452/148‑1 (MQ=255)
                                                                             tgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAtt                                                                     >  8:284085/1‑149 (MQ=255)
                                                                               ccGAACAAATGAAAGAGCCGGACGGCTGCTGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                   <  5:152281/149‑1 (MQ=255)
                                                                               ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAtt                                                                     >  4:111196/1‑147 (MQ=255)
                                                                                 gAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGt                                                                                        >  1:170410/1‑126 (MQ=255)
                                                                                 gAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGt                                                                                        <  2:170410/126‑1 (MQ=255)
                                                                                 gAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTatc                                                                 >  5:253249/1‑147 (MQ=255)
                                                                                     aaaTGAAAGAGCCGGACGGCTGCCGGCTATGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCACTGCAGTAAGACGTTGCAAGGCGGCCAAGGATTTAGCCaa                                                              >  7:300516/1‑148 (MQ=255)
                                                                                            agagCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCt                                                      <  3:111196/149‑1 (MQ=255)
                                                                                                         tGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTCTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgcc                                         >  8:7935/1‑149 (MQ=255)
                                                                                                                     gtgtACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCaa                             >  1:288124/1‑149 (MQ=255)
                                                                                                                       gtacgtaaGTTGACGAGCCGCGCCCGGTCATACCGGGCCCGCAGAGCTCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAaca                           <  4:190203/143‑1 (MQ=255)
                                                                                                                        tACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAaca                           <  2:13031/148‑1 (MQ=255)
                                                                                                                          cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                         >  1:80864/1‑148 (MQ=255)
                                                                                                                          cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                         <  2:80864/148‑1 (MQ=255)
                                                                                                                          cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTCTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                         <  7:7935/148‑1 (MQ=255)
                                                                                                                           tGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCCCAGAGCGCTACTGGGTCAACACGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTGATGAAGCGACGCGACCAGCCCGTGCAACACttc                       >  8:120322/1‑147 (MQ=255)
                                                                                                                             aaGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGa                     >  1:172304/1‑149 (MQ=255)
                                                                                                                                 tGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGc                 >  8:244427/1‑149 (MQ=255)
                                                                                                                                  gACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                   <  7:37026/98‑1 (MQ=255)
                                                                                                                                  gACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                   <  7:42044/98‑1 (MQ=255)
                                                                                                                                  gACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                   >  8:37026/1‑98 (MQ=255)
                                                                                                                                  gACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                   >  8:42044/1‑98 (MQ=255)
                                                                                                                                     tctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgcg                                           <  7:134898/115‑1 (MQ=255)
                                                                                                                                     tctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgcg                                           >  8:134898/5‑119 (MQ=255)
                                                                                                                                     tctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCg             >  8:49130/5‑149 (MQ=255)
                                                                                                                                      aGCCGGGCCCGGTCATACCGGTCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACTTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTTCAACACTGCGACGGCCTCGc            <  7:120322/149‑1 (MQ=255)
                                                                                                                                          gggCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGaa        <  7:244427/149‑1 (MQ=255)
                                                                                                                                                 gTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTcc  >  1:16672/1‑148 (MQ=255)
                                                                                                                                                  tCATACCGGGCCCGCAGAGCGCTACTGGTTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc                                            <  5:59083/105‑1 (MQ=255)
                                                                                                                                                  tCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc                                            >  6:59083/1‑105 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCC  >  NC_002947/4779260‑4779552

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