Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,523,005 T→A 8.4% K623M (AAG→ATG) ‡ yhjG ← protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,523,0050TA 8.4% 155.6 / 3.1 56K623M (AAG→ATG) ‡yhjGprotein YhjG
Reads supporting (aligned to +/- strand):  ref base T (22/29);  new base A (2/3);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGAAGA  >  NC_002947/5522858‑5523152
                                                                                                                                                    |                                                                                                                                                   
aacaacGGTGTGCACTGGTTCGGGTCGATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCt                                                                                                                                                     >  1:197586/1‑149 (MQ=255)
 acaacGGTGTGCACTGGTTCGTGTCATACCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTg                                                                                                                                                   <  8:84599/149‑1 (MQ=255)
          tgCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgc                                                                                                                                           <  3:294731/148‑2 (MQ=255)
           gCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgccg                                                                                                                                         >  6:204960/1‑146 (MQ=255)
            cACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTgg                                                                                                                                                  <  6:363612/139‑1 (MQ=255)
            cACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTgg                                                                                                                                                  >  5:363612/1‑139 (MQ=255)
            cACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgccg                                                                                                                                         <  5:204960/148‑4 (MQ=255)
            cACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCCGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTcc                                                                                                                                        >  5:28910/1‑149 (MQ=255)
              cTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCcgc                                                                                                                                      <  1:371070/149‑1 (MQ=255)
                    cgtgTCATCCCCGCTGCTGGGCGCAATCAGCGCCATCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGcaca                                                                                                                                 <  2:173523/145‑1 (MQ=255)
                       gTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATa                                                                                                                               >  3:124864/1‑147 (MQ=255)
                       gTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgc                                                                                                                                           <  7:351260/135‑2 (MQ=255)
                       gTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgc                                                                                                                                           >  8:351260/1‑134 (MQ=255)
                           aTCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAggg                                                                                                                          <  6:28910/149‑1 (MQ=255)
                            tCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCcgcg                                                                                                                                     <  3:99768/137‑1 (MQ=255)
                            tCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCcgcg                                                                                                                                     >  4:99768/1‑137 (MQ=255)
                            tCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAgggg                                                                                                                         <  5:186417/149‑1 (MQ=255)
                                     ctgGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAg                                                                                                                >  6:46793/1‑149 (MQ=255)
                                        ggCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGc                                                                                                               >  4:39140/1‑147 (MQ=255)
                                        ggCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGc                                                                                                               <  4:329914/147‑1 (MQ=255)
                                        ggCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCga                                                                                                             <  2:197586/149‑1 (MQ=255)
                                         gcgcAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCgag                                                                                                            <  3:39140/149‑1 (MQ=255)
                                            cAATCAGCGCCAGCAACCCCGCCGCCGGCCCTTCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCTGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGaac                                                                                                         <  3:92751/149‑1 (MQ=255)
                                               tCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGaacaac                                                                                                      >  1:277772/1‑149 (MQ=255)
                                                 aGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAAcc                                                                                                     <  8:152845/148‑1 (MQ=255)
                                                      cAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGcc                                                                                               <  2:170003/149‑1 (MQ=255)
                                                              ccgccgccgGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGCACAACCGCCGCCCTTTCGAt                                                                                       >  5:182494/1‑149 (MQ=255)
                                                              ccgccgccgGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGAt                                                                                       >  6:135037/1‑149 (MQ=255)
                                                                 ccgccgGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGAtt                                                                                      >  2:254474/1‑147 (MQ=255)
                                                                    ccgGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAggg                                                                                 >  8:33071/1‑149 (MQ=255)
                                                                          cTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGt                                                                               >  6:388050/1‑145 (MQ=255)
                                                                          cTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGt                                                                               <  5:388050/145‑1 (MQ=255)
                                                                          cTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATg                                                                           >  8:174944/1‑149 (MQ=255)
                                                                           tGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGACCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGAGATga                                                                          >  8:136279/1‑148 (MQ=255)
                                                                             ccaccaCCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTc                                                                        <  2:277772/149‑1 (MQ=255)
                                                                             ccaccaCCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTc                                                                        <  5:135037/149‑1 (MQ=255)
                                                                               accaccCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCgaga                                                                                                           >  6:190925/1‑112 (MQ=255)
                                                                               accaccCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCgaga                                                                                                           <  5:190925/112‑1 (MQ=255)
                                                                                    ccccAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCa                                                                 <  5:46793/149‑1 (MQ=255)
                                                                                      ccAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGTGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGt                                                               <  7:136279/149‑1 (MQ=255)
                                                                                      ccAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGt                                                               >  8:189393/1‑149 (MQ=255)
                                                                                       cAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGtt                                                              <  7:189393/149‑1 (MQ=255)
                                                                                        aGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTg                                                             <  4:124864/149‑1 (MQ=255)
                                                                                                    ccTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCttt                                                                                           <  7:84903/107‑1 (MQ=255)
                                                                                                    ccTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCttt                                                                                           >  8:84903/1‑107 (MQ=255)
                                                                                                       gcgcCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGtt                                              >  6:299193/1‑149 (MQ=255)
                                                                                                           cACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTc                                                                                          >  8:157413/1‑101 (MQ=255)
                                                                                                           cACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTc                                                                                          <  7:157413/101‑1 (MQ=255)
                                                                                                           cACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATACAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCgcc                                          >  6:181996/1‑149 (MQ=255)
                                                                                                                aaCGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCGCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCgccgtgcc                                     <  7:33071/149‑1 (MQ=255)
                                                                                                                     cAGCGGCAATTCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGa                                <  4:225513/149‑1 (MQ=255)
                                                                                                                      aGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGAt                               <  1:254474/149‑1 (MQ=255)
                                                                                                                                 cTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAGCAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATcg                    >  8:305117/1‑149 (MQ=255)
                                                                                                                                        ccGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTtctc             <  5:299193/149‑1 (MQ=255)
                                                                                                                                         cGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCg            >  7:28168/1‑149 (MQ=255)
                                                                                                                                                 ggCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGa     <  7:174944/148‑1 (MQ=255)
                                                                                                                                                   cTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAGCAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGAAGa  <  7:305117/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGAAGA  >  NC_002947/5522858‑5523152

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: