Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,753,542	C→A	13.2%	intergenic (+47/‑319)	glnA → / → glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,753,542	0	C	A	13.2%	152.0 / 9.7	53	intergenic (+47/‑319)	glnA/glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base A (3/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCTG  >  NC_002947/5753395‑5753689                                                                                                                                                    |                                                                                                                                                    acgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACGTGATCAAAACGGCCTCCCTCgg                                                                                                                                                   >  2:1108536/1‑150 (MQ=255) acgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   <  2:889174/150‑1 (MQ=255) acgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   <  2:257763/150‑1 (MQ=255)         aTCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   >  1:548034/1‑142 (MQ=255)         aTCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   <  2:548034/142‑1 (MQ=255)         aTCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGt                                                                                                                                           >  2:465455/1‑150 (MQ=255)         aTCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGt                                                                                                                                           >  1:330830/1‑150 (MQ=255)          tCGATGCCTTCATCGAGCTTAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtt                                                                                                                                          <  1:843748/150‑1 (MQ=255)           cGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCATACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttt                                                                                                                                         <  2:1029259/150‑1 (MQ=255)               gCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGt                                                                                                                                     <  2:799927/150‑1 (MQ=255)                    cATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCTTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGAACAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTg                                                                                                                                <  1:1149265/150‑1 (MQ=255)                     aTCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTgc                                                                                                                               >  1:1137169/1‑150 (MQ=255)                              aaGAGCGAAGAAGAAATCAAGGTCCATACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGc                                                                                                                      >  2:273751/1‑150 (MQ=255)                                  gCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGaa                                                                                                                  <  2:182071/150‑1 (MQ=255)                                   cGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGACTCCTGACTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAg                                                                                                                 <  1:884665/150‑1 (MQ=255)                                       gaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACGTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgcc                                                                                                             <  1:1108536/150‑1 (MQ=255)                                       gaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgcc                                                                                                             <  2:1137169/150‑1 (MQ=255)                                          gaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCt                                                                                                          >  2:1423573/1‑150 (MQ=255)                                              tCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGAtt                                                                                                      >  1:916677/1‑150 (MQ=255)                                                           cTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   <  1:635285/92‑1 (MQ=255)                                                           cTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   >  2:635285/1‑92 (MQ=255)                                                           cTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCCCGGGAGGCCGTTTTTGTTTCTgc                                                                                                                               <  2:18959/112‑1 (MQ=255)                                                           cTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCCCGGGAGGCCGTTTTTGTTTCTgc                                                                                                                               >  1:18959/1‑112 (MQ=255)                                                                       gCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCACCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCtt                                                                             >  1:277306/1‑150 (MQ=255)                                                                               aCGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   <  2:1418892/72‑1 (MQ=255)                                                                               aCGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                   >  1:1418892/1‑72 (MQ=255)                                                                                 gAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGt                                                                                               >  2:1221835/1‑122 (MQ=255)                                                                                 gAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGt                                                                                               <  1:1221835/122‑1 (MQ=255)                                                                                  aGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAt                                                                  <  1:465455/150‑1 (MQ=255)                                                                                   gCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCACCAAtt                                                                 <  2:916677/150‑1 (MQ=255)                                                                                             aGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGc                                                       <  2:330830/150‑1 (MQ=255)                                                                                               ctgctgATCTGATTGGCGCCTCGCACGGTCTTCATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGccc                                                     >  1:1107535/1‑150 (MQ=255)                                                                                                 gctgATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTc                                                   >  2:440173/1‑150 (MQ=255)                                                                                                        tgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGc                                                                                                                                              >  1:413191/1‑52 (MQ=37)                                                                                                        tgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGc                                                                                                                                              <  2:413191/52‑1 (MQ=37)                                                                                                         gatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTc                                           <  1:273751/150‑1 (MQ=255)                                                                                                           ttGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGGCGCTGATTCGCCGGTCACGCTATGCTCTATCTTGGTGCATGAATTCACACCATGCCCTCAATTGGTccc                                         >  1:654714/1‑150 (MQ=255)                                                                                                              gcgcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                            >  2:610992/1‑112 (MQ=255)                                                                                                              gcgcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                            >  2:611048/1‑112 (MQ=255)                                                                                                              gcgcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                            <  1:611048/112‑1 (MQ=255)                                                                                                              gcgcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                            <  1:610992/112‑1 (MQ=255)                                                                                                                         cGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCAt                                                                      >  1:460086/1‑107 (MQ=255)                                                                                                                         cGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCAt                                                                      <  2:460086/107‑1 (MQ=255)                                                                                                                         cGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAgg                           >  2:990852/1‑150 (MQ=255)                                                                                                                               cATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAAc                     <  1:440173/150‑1 (MQ=255)                                                                                                                                aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACg                    >  2:1083694/1‑150 (MQ=255)                                                                                                                                  gATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGtt                  >  2:1183943/1‑150 (MQ=255)                                                                                                                                   aTCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGttt                 >  2:1082920/1‑150 (MQ=255)                                                                                                                                      aaaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTACACGTTTCAc              >  2:619109/1‑150 (MQ=255)                                                                                                                                        aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctct                                                                                  <  2:1244250/80‑1 (MQ=38)                                                                                                                                        aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctct                                                                                  >  1:1244250/1‑80 (MQ=38)                                                                                                                                        aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGt                                                                          <  2:1048161/88‑1 (MQ=255)                                                                                                                                        aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGt                                                                          >  1:1048161/1‑88 (MQ=255)                                                                                                                                              ctccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTg                                  >  2:498443/1‑122 (MQ=255)                                                                                                                                              ctccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTg                                  <  1:498443/122‑1 (MQ=255)                                                                                                                                               cccctcGGGAGGCCGTTTTTGTTTCTGCGCAACTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTCTGCTCTATCTTGCTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACttt     <  2:654714/149‑1 (MQ=255)                                                                                                                                                ccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAACAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc    <  1:1143636/150‑1 (MQ=255)                                                                                                                                                 ccttGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCt   >  1:1035273/1‑150 (MQ=255)                                                                                                                                                  ctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGGCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCTg  >  1:968149/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCTG  >  NC_002947/5753395‑5753689

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑