Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,523,008	G→C	9.4%	A622G (GCC→GGC)	yhjG ←	protein YhjG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,523,008	0	G	C	9.4%	162.8 / 5.5	53	A622G (GCC→GGC)	yhjG	protein YhjG
Reads supporting (aligned to +/- strand):  ref base G (21/27);  new base C (2/3);  total (23/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.96e-01

TGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGAAGA  >  NC_002947/5522868‑5523152                                                                                                                                             |                                                                                                                                                 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                                                                               <  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                                                                                                           <  5:190925/112‑1 (MQ=255)                                                                     accaccCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCgaga                                                                                                           >  6:190925/1‑112 (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                                                                                          <  7:157413/101‑1 (MQ=255)                                                                                                 cACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTc                                                                                          >  8:157413/1‑101 (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)                                                                                                                                             |                                                                                                                                                 TGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCACCCCCAGCGCAACCATCCCTGCGCCACGCAACGCCAGCGGCAATGCCTGCACACCGGCGTTCGGCTTGGCGAACGGTCCGCGCACATACAGGGGGGAACGCAGCGAGAACAACCGCAGCCCTTTCGATTCAGGGGTGATGTTCAAATCCAGTTGCTCGGAAGCGAAGTTCGCCGTGCCGTCGATGTAGATGATCGCGTTCTCGGTATCGAAGA  >  NC_002947/5522868‑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: ‑