Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,523,003	G→C	8.9%	P624A (CCG→GCG)	yhjG ←	protein YhjG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,523,003	0	G	C	8.9%	169.1 / 6.0	56	P624A (CCG→GCG)	yhjG	protein YhjG
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (2/3);  total (26/30)
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.46e-01

GCAACAACGGTGTGCACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGA  >  NC_002947/5522856‑5523149                                                                                                                                                     |                                                                                                                                                   gcAACAACGGTGTGCACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGc                                                                                                                                                   >  7:152845/1‑149 (MQ=255)   aacaacGGTGTGCACTGGTTCGGGTCGATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCt                                                                                                                                                  >  1:197586/1‑149 (MQ=255)    acaacGGTGTGCACTGGTTCGTGTC‑ATACCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTg                                                                                                                                                <  8:84599/149‑1 (MQ=255)             tgCACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgc                                                                                                                                        <  3:294731/148‑2 (MQ=255)              gCACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgccg                                                                                                                                      >  6:204960/1‑146 (MQ=255)               cACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTgg                                                                                                                                               <  6:363612/139‑1 (MQ=255)               cACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTgg                                                                                                                                               >  5:363612/1‑139 (MQ=255)               cACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgccg                                                                                                                                      <  5:204960/148‑4 (MQ=255)               cACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCCGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTcc                                                                                                                                     >  5:28910/1‑149 (MQ=255)                 cTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCcgc                                                                                                                                   <  1:371070/149‑1 (MQ=255)                       cgtgTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCATCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGcaca                                                                                                                              <  2:173523/145‑1 (MQ=255)                          gTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATa                                                                                                                            >  3:124864/1‑147 (MQ=255)                          gTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgc                                                                                                                                        <  7:351260/135‑2 (MQ=255)                          gTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGCCAAGCCGAACgc                                                                                                                                        >  8:351260/1‑134 (MQ=255)                              aTCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAggg                                                                                                                       <  6:28910/149‑1 (MQ=255)                               tCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCcgcg                                                                                                                                  >  4:99768/1‑137 (MQ=255)                               tCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCcgcg                                                                                                                                  <  3:99768/137‑1 (MQ=255)                               tCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAgggg                                                                                                                      <  5:186417/149‑1 (MQ=255)                                        ctgGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAg                                                                                                             >  6:46793/1‑149 (MQ=255)                                           ggCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGc                                                                                                            <  4:329914/147‑1 (MQ=255)                                           ggCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGc                                                                                                            >  4:39140/1‑147 (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                                                                            <  5:388050/145‑1 (MQ=255)                                                                             cTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGt                                                                            >  6:388050/1‑145 (MQ=255)                                                                             cTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATg                                                                        >  8:174944/1‑149 (MQ=255)                                                                              tGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGACCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGAGATga                                                                       >  8:136279/1‑148 (MQ=255)                                                                                ccaccaCCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTc                                                                     <  5:135037/149‑1 (MQ=255)                                                                                ccaccaCCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTc                                                                     <  2:277772/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                                                                                        >  8:84903/1‑107 (MQ=255)                                                                                                       ccTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCttt                                                                                        <  7:84903/107‑1 (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)                                                                                                                                                     |                                                                                                                                                   GCAACAACGGTGTGCACTGGTTCGGGTC‑ATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGA  >  NC_002947/5522856‑5523149

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