Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,782,594 C→T 9.5% G1574S (GGC→AGC)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,782,5940CT 9.5% 114.0 / 2.8 42G1574S (GGC→AGC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (2/2);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATG  >  NC_002947/4782454‑4782742
                                                                                                                                            |                                                                                                                                                    
cACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCg                                                                                                                                              >  7:76340/1‑149 (MQ=255)
  ccATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGt                                                                                                                                            <  3:327024/149‑1 (MQ=255)
       cgcTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAg                                                                                                                                       <  3:254079/149‑1 (MQ=255)
             cccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGCCCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgt                                                                                                                                 >  7:168877/1‑149 (MQ=255)
                 ggTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGc                                                                                                                             <  7:187489/149‑1 (MQ=255)
                   tCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACCGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                           >  8:74862/1‑149 (MQ=255)
                        gCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcagca                                                                                                                      <  1:142019/149‑1 (MQ=255)
                        gCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcagca                                                                                                                      <  4:4458/149‑1 (MQ=255)
                            aaCGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGccc                                                                                                                  >  3:38562/1‑149 (MQ=255)
                                    gcgcAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCa                                                                                                          >  7:327280/1‑149 (MQ=255)
                                        agcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCgg                                                                                                      <  2:38704/149‑1 (MQ=255)
                                        agcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGAATCAGCAGCCCTTCCTGCATCgg                                                                                                      >  3:326163/1‑149 (MQ=255)
                                              tCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACaa                                                                                                <  8:327280/149‑1 (MQ=255)
                                                     gcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGa                                                                                            >  3:161113/1‑146 (MQ=255)
                                                     gcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACCGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAg                                                                                         <  7:74862/149‑1 (MQ=255)
                                                       gcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGAt                                                                                       <  6:364364/149‑1 (MQ=255)
                                                        cAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTGCCTGCATCGGCGACAACGgagagaag                                                                                      >  8:302432/1‑147 (MQ=255)
                                                            gACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGa                                                                                                   >  2:86085/1‑132 (MQ=255)
                                                            gACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGCTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGa                                                                                                   <  1:86085/132‑1 (MQ=255)
                                                              ccGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTc                                                                                 >  6:77908/1‑148 (MQ=255)
                                                               cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                               >  5:359847/1‑149 (MQ=255)
                                                               cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                               >  7:142800/1‑149 (MQ=255)
                                                               cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                               >  7:137773/1‑149 (MQ=255)
                                                                              gCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGcc                                                                 <  5:77908/148‑1 (MQ=255)
                                                                                  aaCGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCAcc                                                              <  6:359847/147‑1 (MQ=255)
                                                                                                            gACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCt                                    >  8:71673/1‑147 (MQ=255)
                                                                                                              cgcgGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGc                                <  4:161113/149‑1 (MQ=255)
                                                                                                                 ggTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCg                             >  7:294289/1‑149 (MQ=255)
                                                                                                                   tAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                           >  7:48889/1‑149 (MQ=255)
                                                                                                                         tACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTgg                     >  6:66677/1‑149 (MQ=255)
                                                                                                                             gCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGt                       >  5:39407/1‑143 (MQ=255)
                                                                                                                             gCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGt                       <  6:39407/143‑1 (MQ=255)
                                                                                                                               aTCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCACAGCCCGGCCAGTGGGAAGTc               <  4:38562/149‑1 (MQ=255)
                                                                                                                                 cAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                                   >  7:342457/1‑111 (MQ=255)
                                                                                                                                   gatagatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATACAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCtt           >  3:121172/1‑149 (MQ=255)
                                                                                                                                   gatagatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCCCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCtt           >  6:161373/1‑149 (MQ=255)
                                                                                                                                    caccttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCg                                                      <  4:145399/100‑1 (MQ=255)
                                                                                                                                    caccttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCg                                                      >  3:145399/6‑105 (MQ=255)
                                                                                                                                       cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCt                                    >  6:82496/3‑120 (MQ=255)
                                                                                                                                       cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCt                                    <  5:82496/118‑1 (MQ=255)
                                                                                                                                          gCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTg          <  3:48178/143‑1 (MQ=255)
                                                                                                                                          gCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGACTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTg          >  4:48178/1‑143 (MQ=255)
                                                                                                                                            ccGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATg  >  5:120079/1‑149 (MQ=255)
                                                                                                                                            ccGAACCCGGTTCGAGCAAGGTGTGAATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATg  <  4:326163/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATG  >  NC_002947/4782454‑4782742

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