Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,779,409 T→G 10.0% 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,4090TG10.0% 142.4 / 7.0 50intergenic (‑143/+18)pvdD/pvdJnon‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base G (3/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAGG  >  NC_002947/4779269‑4779557
                                                                                                                                            |                                                                                                                                                    
gAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCggg                                                                                                                                                <  1:199642/147‑1 (MQ=255)
gAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATCCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCggg                                                                                                                                                <  1:11005/147‑1 (MQ=255)
    aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                          <  5:24415/149‑1 (MQ=255)
    aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                          <  1:127283/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                              <  2:80864/148‑1 (MQ=255)
                                                                                                                 cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                              >  1:80864/1‑148 (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                                                                        >  8:37026/1‑98 (MQ=255)
                                                                                                                         gACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                        <  7:42044/98‑1 (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)
                                                                                                                                            tACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAgg  <  6:253249/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAGG  >  NC_002947/4779269‑4779557

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: