Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,782,601 C→G 9.1% P1571P (CCG→CCC)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,782,6010CG 9.1% 131.0 / 3.6 44P1571P (CCG→CCC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (23/17);  new base G (2/2);  total (25/19)
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.35e-01

CACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTG  >  NC_002947/4782454‑4782747
                                                                                                                                                   |                                                                                                                                                  
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:137773/1‑149 (MQ=255)
                                                               cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                                    >  7:142800/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                                         <  5:82496/118‑1 (MQ=255)
                                                                                                                                       cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCt                                         >  6:82496/3‑120 (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)
                                                                                                                                                 cccGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTg  <  8:137773/149‑1 (MQ=255)
                                                                                                                                                 cccGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTg  <  8:142800/149‑1 (MQ=255)
                                                                                                                                                  ccGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGt                     >  2:227461/1‑129 (MQ=255)
                                                                                                                                                  ccGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGt                     <  1:227461/129‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCGTAATGCCTTG  >  NC_002947/4782454‑4782747

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