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

CTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCC  >  NC_002947/4779260‑4779553
                                                                                                                                                   |                                                                                                                                                  
cTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCa                                                                                                                                                   <  3:365317/149‑1 (MQ=255)
 ttGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCa                                                                                                                                                   >  4:353754/1‑148 (MQ=255)
 ttGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCa                                                                                                                                                   >  3:1924/1‑148 (MQ=255)
       ttGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCCGGCATACCgcg                                                                                                                                            >  7:134546/1‑147 (MQ=255)
         gAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGTTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCggg                                                                                                                                            <  7:289629/147‑1 (MQ=255)
            cATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTATGACCGGGCCCgg                                                                                                                                       <  5:373501/149‑2 (MQ=255)
             aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                      <  1:330056/149‑1 (MQ=255)
             aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                      <  3:362531/149‑1 (MQ=255)
                  cTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGACAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCag                                                                                                                                     <  1:122728/145‑1 (MQ=255)
                  cTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAgcg                                                                                                                                 >  2:86944/1‑149 (MQ=255)
                     ccAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTa                                                                                                                              >  8:313172/1‑149 (MQ=255)
                        gggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTg                                                                                                                           >  3:185301/1‑149 (MQ=255)
                                gTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACt                                                                                                                   >  4:282210/1‑149 (MQ=255)
                                         tGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtct                                                                                                          >  4:304956/1‑149 (MQ=255)
                                         tGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCAGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtct                                                                                                          >  7:10269/1‑149 (MQ=255)
                                             gAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATctcctc                                                                                                      <  4:1924/149‑1 (MQ=255)
                                                   gCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCaga                                                                                                                                    <  1:147952/113‑1 (MQ=255)
                                                   gCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCaga                                                                                                                                    >  2:147952/1‑113 (MQ=255)
                                                   gCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAg                                                                                                 >  1:215061/1‑148 (MQ=255)
                                                        cccGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGAc                                                                                           <  3:282210/149‑1 (MQ=255)
                                                         ccGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTATGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACAGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACg                                                                                          >  4:210392/1‑149 (MQ=255)
                                                           gcgACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGtt                                                                                        >  3:85606/1‑149 (MQ=255)
                                                           gcgACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGtt                                                                                        >  5:282986/1‑149 (MQ=255)
                                                             gACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGttt                                                                                       >  7:188404/1‑148 (MQ=255)
                                                                               ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                    <  2:215061/149‑1 (MQ=255)
                                                                               ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa                                                                    <  4:85606/149‑1 (MQ=255)
                                                                                     aaaTGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAAc                                                              >  3:24683/1‑149 (MQ=255)
                                                                                      aaTGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACt                                                             >  8:77094/1‑149 (MQ=255)
                                                                                       aTGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTc                                                            <  7:77094/149‑1 (MQ=255)
                                                                                         gAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGt                                                          >  4:287275/1‑149 (MQ=255)
                                                                                           aaGAGCCGGAAGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCTGGCATACCGGGCCCGCAGAGCGCTACTGTGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTc                                                         <  4:259447/148‑1 (MQ=255)
                                                                                            agagCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCt                                                       <  3:210392/149‑1 (MQ=255)
                                                                                               gCCGGACGGCTGCCGGCTCTGGGGGTACTGAAGTGGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGGCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGaa                                                    <  8:10269/149‑1 (MQ=255)
                                                                                                    aCGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAc                                               <  4:286490/149‑1 (MQ=255)
                                                                                                    aCGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAc                                               <  1:86944/149‑1 (MQ=255)
                                                                                                       gCTGCCGGCTCTGGGTGTACTGATGTTGACGAGCCGGGCACGTTAATACCGGGCCCGCAGATCGATACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAATTATTTAGCCAACTCGTCCTTAATCGAcgc                                             <  4:185301/148‑1 (MQ=255)
                                                                                                        cTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgcg                                            <  8:148874/148‑1 (MQ=255)
                                                                                                        cTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgcg                                            >  7:148874/1‑148 (MQ=255)
                                                                                                          gCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGc                                                                  <  7:114774/124‑1 (MQ=255)
                                                                                                                   gggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTg                                 >  2:353433/1‑148 (MQ=255)
                                                                                                                   gggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTg                                 >  6:154169/1‑148 (MQ=255)
                                                                                                                     gtgtACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCa                               >  8:182347/1‑148 (MQ=255)
                                                                                                                        tACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAacac                           <  6:282986/149‑1 (MQ=255)
                                                                                                                             aaGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGc                        <  7:313172/147‑1 (MQ=255)
                                                                                                                               gTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGCTCTCCTCTGCAGTAAGACGTTTCAAGGCTGCCAAGGATTTAGCCACATCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGaag                    >  2:306529/1‑147 (MQ=255)
                                                                                                                                ttGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACgg                   >  4:173037/1‑149 (MQ=255)
                                                                                                                                    cGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTc               >  6:367377/1‑149 (MQ=255)
                                                                                                                                      ctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACg                    <  7:37920/139‑1 (MQ=255)
                                                                                                                                      ctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACg                    >  8:37920/4‑142 (MQ=255)
                                                                                                                                         cGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAACTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCACCTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTg           >  3:42787/1‑148 (MQ=255)
                                                                                                                                                  tCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTccc  >  4:207623/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCC  >  NC_002947/4779260‑4779553

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: