Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,785,772 G→A 8.0% R514R (CGC→CGT pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,785,7720GA 8.0% 137.4 / 2.4 50R514R (CGC→CGTpvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CGGTCACCCATTCGATGGAGCGCTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAGGTGCCGCGACCATGTCCTGC  >  NC_002947/4785624‑4785916
                                                                                                                                                     |                                                                                                                                                
cGGTCACCCATTCGATGGAGCGCTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCgggg                                                                                                                                                  <  4:145984/149‑1 (MQ=255)
     aCCCATTCGATGGAGCGCTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTc                                                                                                                                             >  4:103906/1‑149 (MQ=255)
             gATGGAGCGCTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCAt                                                                                                                                     >  8:120006/1‑149 (MQ=255)
             gATGGAGCGCTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCAt                                                                                                                                     >  2:158546/1‑149 (MQ=255)
                     gcTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAgc                                                                                                                             <  8:222812/149‑1 (MQ=255)
                         ttGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGACGCTCGCCGGCATTGGCCAGCGccc                                                                                                                         >  7:55799/1‑149 (MQ=255)
                         ttGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGACGCTCGCCGGCATTGGCCAGCGccc                                                                                                                         >  7:271476/1‑149 (MQ=255)
                              acaAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCcggcg                                                                                                                     >  3:155946/1‑148 (MQ=255)
                              acaAGGCAACCACTGCGCCCGGCTTCACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGCCCAGTGCCcggcg                                                                                                                     >  3:50131/1‑148 (MQ=255)
                                aaGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGaa                                                                                                                  <  4:148250/149‑1 (MQ=255)
                                 aGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAAc                                                                                                                 <  3:232457/149‑1 (MQ=255)
                                    cAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAgg                                                                                                              <  1:233506/149‑1 (MQ=255)
                                     aaCCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGAGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGc                                                                                                             >  5:295716/1‑149 (MQ=255)
                                      aCCACTGCGCCCGGGTTGACCCCCATTGCGAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGCTCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACTGGCt                                                                                                            <  8:288590/149‑1 (MQ=255)
                                      aCCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTTGATAGAGCTGGGCATAGCTCAACTGCTATGCGCCCGCCCGCAGCGCGGGGCGCTAGCCGGCATTGGCCAGCGCCCGGCGGAACAGGc                                                                                                             >  1:242638/1‑149 (MQ=255)
                                        cACTGCGCCCGGGTTGACGGCCATTGCCTGCAGGTGGTGCGCCAACTGGGTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCTTTGGCCAGTGCCCGGGGGAACAGGCTga                                                                                                          <  4:50131/149‑1 (MQ=255)
                                                ccGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTcc                                                                                                  >  4:161487/1‑149 (MQ=255)
                                                 cGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCg                                                                                                 <  4:155946/149‑1 (MQ=255)
                                                  ggCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGt                                                                                                >  8:177329/1‑149 (MQ=255)
                                                   gCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTg                                                                                               <  3:2319/149‑1 (MQ=255)
                                                       gACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCaa                                                                                            <  2:242638/148‑1 (MQ=255)
                                                        aCGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAg                                                                                           <  3:103906/148‑1 (MQ=255)
                                                            ccATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAg                                                                                      >  1:216539/1‑149 (MQ=255)
                                                                   cagcagGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCgg                                                                               >  2:93207/1‑149 (MQ=255)
                                                                         gtggtgCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCaa                                                                                            <  5:50362/130‑1 (MQ=255)
                                                                         gtggtgCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCaa                                                                                            >  6:50362/1‑130 (MQ=255)
                                                                          tggtgCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCAcc                                                                        <  3:161487/149‑1 (MQ=255)
                                                                            gtgCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGACGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGt                                                                      <  8:55799/149‑1 (MQ=255)
                                                                            gtgCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGACGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGt                                                                      <  8:271476/149‑1 (MQ=255)
                                                                                      tGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCg                                                             >  3:267791/1‑148 (MQ=255)
                                                                                        gTTCGAACGCTGATCAAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCAt                                                          >  3:110287/1‑149 (MQ=255)
                                                                                         ttCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATc                                                         <  1:158546/149‑1 (MQ=255)
                                                                                           cGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGt                                                                      >  5:237338/1‑134 (MQ=255)
                                                                                           cGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGt                                                                      <  6:237338/134‑1 (MQ=255)
                                                                                             aaCGCTGATCGAGCTTGGCATAGCTCAACTGCTCTGGGCCCGCAGGCAGGGGGGGGGGCTCGCAGGCATTGGCCAGCGCCCGGCGGa                                                                                                                   <  2:300363/87‑1 (MQ=255)
                                                                                             aaCGCTGATCGAGCTGGGCATAGCTCAACTGCTCAGCGCCCGCACGCAGCGCGCGGCGCTCGCCCGCATTGGCCAGCGCCCGGCGGa                                                                                                                   >  1:300363/1‑87 (MQ=255)
                                                                                               cGCTGATCAAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATc                                                   <  4:110287/149‑1 (MQ=255)
                                                                                                gCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCa                                                  <  7:177329/149‑1 (MQ=255)
                                                                                                gCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCa                                                  >  7:214531/1‑149 (MQ=255)
                                                                                                              ggCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTc                                       <  8:214531/147‑1 (MQ=255)
                                                                                                              ggCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTc                                       <  1:93207/147‑1 (MQ=255)
                                                                                                                 aTAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCcag                                    >  7:50785/1‑147 (MQ=255)
                                                                                                                 aTAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCcag                                    >  1:29743/1‑147 (MQ=255)
                                                                                                                             ctctGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAg                      >  4:12960/1‑149 (MQ=255)
                                                                                                                              tctGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAgg                     >  4:102579/1‑149 (MQ=255)
                                                                                                                                      cgcacgcaGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAGGTGCCGCGa             >  8:8993/1‑149 (MQ=255)
                                                                                                                                           gcaGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCATGTGCCGCGACCATg        <  6:295716/149‑1 (MQ=255)
                                                                                                                                           gcaGCGCGGGGCGCTCGCCCGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCCCATCGGCATCACCGACACAGTCCAGCAGCCGGCAGGCAGGTGCCGCGACCAt         >  8:156309/1‑148 (MQ=255)
                                                                                                                                            caGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAGGTGCCGCGACCATGt       >  3:291861/1‑149 (MQ=255)
                                                                                                                                              gcgcGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAGGTGCCGCGACCATGTcc     <  8:50785/149‑1 (MQ=255)
                                                                                                                                                 cgGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAGGTGCCGCGACCATGTCCTGc  <  7:120006/149‑1 (MQ=255)
                                                                                                                                                     gcgcTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCcagc                                   <  2:209052/112‑1 (MQ=255)
                                                                                                                                                     gcgcTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCcagc                                   >  1:209052/1‑112 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGGTCACCCATTCGATGGAGCGCTCTTGGCACAAGGCAACCACTGCGCCCGGCTTGACGCCCATTGCCAGCAGGTGGTGCGCCAACTGGTTCGAACGCTGATCGAGCTGGGCATAGCTCAACTGCTCTGCGCCCGCACGCAGCGCGGGGCGCTCGCCGGCATTGGCCAGCGCCCGGCGGAACAGGCTGAGAAAGTCCGTGCAAGGGAAGCTCGCGGGCTCACCGTGCAAAACCGCATCGGCATCACCGACACAGTCCAGCAGCCGGCTGGCAGGTGCCGCGACCATGTCCTGC  >  NC_002947/4785624‑4785916

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: