Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,782,594 C→T 9.1% 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.1% 148.3 / 4.5 55G1574S (GGC→AGC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (23/27);  new base T (3/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

CTGCACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCG  >  NC_002947/4782451‑4782737
                                                                                                                                               |                                                                                                                                               
cTGCACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAAc                                                                                                                                            >  3:184165/1‑149 (MQ=255)
   cACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCg                                                                                                                                         <  6:72525/149‑1 (MQ=255)
   cACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCg                                                                                                                                         <  4:312396/149‑1 (MQ=255)
   cACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCg                                                                                                                                         <  4:184165/149‑1 (MQ=255)
       aTGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACGGAATCAGATAGATGCCCGCACCCGGTTc                                                                                                                                     >  3:29352/1‑149 (MQ=255)
              ccccacggTCCCAGCAACACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGGGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                              <  4:29352/143‑1 (MQ=255)
               aCCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                              <  8:43694/148‑1 (MQ=255)
               aCCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                              <  6:75060/148‑1 (MQ=255)
                 ccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCTCGAACCGGGTTCGGGCatctat                                                                                                                           <  6:34852/149‑6 (MQ=255)
                 ccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtg                                                                                                                           >  6:1037/1‑149 (MQ=255)
                 ccAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtg                                                                                                                           >  7:289949/1‑149 (MQ=255)
                  cAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtgt                                                                                                                          <  7:236150/149‑1 (MQ=255)
                     gTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCa                                                                                                                       >  3:66889/1‑149 (MQ=255)
                        ccAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                      >  1:178325/1‑147 (MQ=255)
                           gCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATcag                                                                                                                   >  1:141460/1‑147 (MQ=255)
                                  gAGGTGCGCAGCACGTCATGGCGCTCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGTTGTGCATCAGCAGCCCTTc                                                                                                          <  8:78933/149‑1 (MQ=255)
                                         gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtgt                                                                                                                          <  7:311080/126‑1 (MQ=255)
                                         gcagcaCGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGgtgt                                                                                                                          >  8:311080/1‑126 (MQ=255)
                                               cGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGAc                                                                                             <  8:289949/149‑1 (MQ=255)
                                                      gcgcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                   <  7:182012/104‑1 (MQ=255)
                                                      gcgcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGa                                                                                                                                   >  8:182012/1‑104 (MQ=255)
                                                       cgcgcAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATa                                                                                     <  7:335759/149‑1 (MQ=255)
                                                                aCCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                              >  2:279874/1‑99 (MQ=255)
                                                                aCCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAg                                                                                                                              <  1:279874/99‑1 (MQ=255)
                                                                  cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                          >  8:172876/1‑149 (MQ=255)
                                                                  cGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGa                                                                          >  5:222974/1‑149 (MQ=255)
                                                                         cAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCt                                                                                                        >  6:322334/1‑112 (MQ=255)
                                                                            gCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATcgcgcgc                                                                >  6:247536/1‑149 (MQ=255)
                                                                             cTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCagcagc                                                                                                               <  8:323969/101‑1 (MQ=255)
                                                                             cTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCagcagc                                                                                                               >  7:323969/1‑101 (MQ=255)
                                                                             cTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATcgcgcgcg                                                               <  4:66889/149‑1 (MQ=255)
                                                                                 gCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                      <  1:157203/90‑1 (MQ=255)
                                                                                 gCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCAt                                                                                                                      >  2:157203/1‑90 (MQ=255)
                                                                                  cTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCAc                                                          >  5:348471/1‑149 (MQ=255)
                                                                                            cGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGAGCGGCCACCGGCAGGCGa                                                >  4:313760/1‑149 (MQ=255)
                                                                                                aTCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCa                                            <  2:158983/149‑1 (MQ=255)
                                                                                                 tCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAg                                           <  6:348471/149‑1 (MQ=255)
                                                                                                   aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGTGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                         <  2:178325/149‑1 (MQ=255)
                                                                                                   aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                         <  7:320331/149‑1 (MQ=255)
                                                                                                   aaTGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGtt                                         >  3:123407/1‑149 (MQ=255)
                                                                                                        cGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTtgctgc                                    <  4:123407/149‑1 (MQ=255)
                                                                                                         gTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTgctgct                                   <  5:247536/149‑1 (MQ=255)
                                                                                                            gTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGc                                >  8:264567/1‑149 (MQ=255)
                                                                                                                aCGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAg                            <  6:222974/149‑1 (MQ=255)
                                                                                                                    ggTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCg                        <  5:1037/149‑1 (MQ=255)
                                                                                                                        gCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGc                      >  4:285752/1‑147 (MQ=255)
                                                                                                                          aaTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGt                  <  3:313760/149‑1 (MQ=255)
                                                                                                                           aTACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGCTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTg                 >  3:58407/1‑149 (MQ=255)
                                                                                                                            tACTGCATCCGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTgg                <  6:193822/149‑1 (MQ=255)
                                                                                                                            tACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTgg                <  7:172876/149‑1 (MQ=255)
                                                                                                                                      gatagatGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCtt      >  3:43111/1‑149 (MQ=255)
                                                                                                                                        accttGCTCGAACCGGGTTCGGGCATGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                              >  3:63004/5‑107 (MQ=255)
                                                                                                                                        accttGCTCGAACCGGGTTCGGGCATGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                              <  4:63004/103‑1 (MQ=255)
                                                                                                                                          cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCg  >  8:141576/3‑149 (MQ=255)
                                                                                                                                          cttGCTCGAACCGGGTTCGGGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCCATTTCGGCCACAGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCg  >  1:294789/3‑149 (MQ=255)
                                                                                                                                               ccGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                              <  5:341952/100‑1 (MQ=255)
                                                                                                                                               ccGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATc                                              >  6:341952/1‑100 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CTGCACCATGCGCTCACCCAGGTCCCAGCAAAACGAGGTGCGCAGCACGTCATGGCGCGCAATGACCGCTTTCCAGGCTTCGCTGAACGCGCCGACATCAATGGCGTGGTCGACGCGGTAGCAATACTGCATCAGATAGATGCCCGAACCCGGTTCGAGCAAGGTGTGCATCAGCAGCCCTTCCTGCATCGGCGACAACGGATAGATGTCTTCGATATCGCGCGCGGCCACCGGCAGGCGATCCAGTTGCTGCTGGCTCAGCCCGGCCAGTGGGAAGTCGCTTGGCG  >  NC_002947/4782451‑4782737

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