Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,843,996	T→G	E266D (GAA→GAC)	gacS ←	sensor protein GacS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,843,996	0	T	G	97.1%	113.0 / NA	35	E266D (GAA→GAC)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (19/15);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTTTGCAGCAGGTGGGTAAAGCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGG  >  NC_002947/1843853‑1844135                                                                                                                                                |                                                                                                                                            cTTTTTTGCAGCAGGTGGGTAAAGCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGAt                                                                                                                                        <  1:233889/149‑1 (MQ=255)  gtgtttGCAGCAGGTGGGTAAAGCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATg                                                                                                                                       <  2:272741/146‑1 (MQ=255)         cagcagGTGGGTAAAGCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCa                                                                                                                                <  1:79767/149‑1 (MQ=255)            cagGTGGGTAAAGCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGt                                                                                                                             >  2:71160/1‑149 (MQ=255)                        gCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTATGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCCATGGTTTCCAGGATCTGGCGCACga                                                                                                                 >  3:311405/1‑148 (MQ=255)                         ccGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTc                                                                                                                >  7:297706/1‑149 (MQ=255)                                               tGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCtg                                                                                          <  1:321859/149‑1 (MQ=255)                                               tGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCtg                                                                                          >  6:118700/1‑149 (MQ=255)                                               tGCGGATTTCGTGGATCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCAAGGTTATGGCGCAAGTCTTCGGTGGCCTGGTCGATGCtg                                                                                          >  2:237341/1‑149 (MQ=255)                                                   gaattCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGtt                                                                                      <  8:273263/146‑1 (MQ=255)                                                   gATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGtt                                                                                      >  7:2440/1‑149 (MQ=255)                                                                aTGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGt                                                                          <  8:2440/148‑1 (MQ=255)                                                                aTGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTg                                                                         <  5:330239/149‑1 (MQ=255)                                                                 tGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTgg                                                                        <  7:246379/149‑1 (MQ=255)                                                                 tGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTgg                                                                        <  8:297706/149‑1 (MQ=255)                                                                  gTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTgg                                                                        <  2:315318/148‑1 (MQ=255)                                                                  gTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTgg                                                                        >  1:315318/1‑148 (MQ=255)                                                                   ttGGCCAGGAACTCCGACTTGATGCGGCTTGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGc                                                                      <  1:237341/149‑1 (MQ=255)                                                                    tGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGgcg                                                                                                                      >  3:40485/1‑100 (MQ=255)                                                                      gCCAGGAACTACGACTTGATGCGGCTGGCCTCAAGTGACTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGACGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGACGATGATGTGTTGCAGTTCTTCCTGGGCAtt                                                                   >  3:145219/1‑149 (MQ=255)                                                                               tCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGt                                                          >  1:24609/1‑149 (MQ=255)                                                                                   aCTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCg                                                      <  2:24609/149‑1 (MQ=255)                                                                                     ttGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGc                                                    >  4:295881/1‑149 (MQ=255)                                                                                                  ccTCATGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCTTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATg                                       <  7:270051/149‑1 (MQ=255)                                                                                                        gTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTgg                                 >  1:206908/1‑149 (MQ=255)                                                                                                            ctctTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCa                             <  2:206908/149‑1 (MQ=255)                                                                                                            ctctTTGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGAAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCa                             >  5:25663/1‑149 (MQ=255)                                                                                                               tttGCGCGCCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGtt                          >  4:80524/1‑149 (MQ=255)                                                                                                                    gcgcCATGTCCAGCTAGATGTTCTGGATGTGGTTGGTTTGCTGGTTCTGGCGCACGTCTTCGTTGGCATTGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGAAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTcc                     <  4:145219/149‑1 (MQ=255)                                                                                                                      gcCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGc                                                               >  7:252622/1‑105 (MQ=255)                                                                                                                      gcCATGTCCAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGc                                                               <  8:252622/105‑1 (MQ=255)                                                                                                                              cAGCTCGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGc           >  6:141527/1‑149 (MQ=255)                                                                                                                                   cGATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGcc       >  3:298961/1‑148 (MQ=255)                                                                                                                                    gATGTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCa     >  8:41278/1‑149 (MQ=255)                                                                                                                                       gTTCTGGATGTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATgg  >  7:27857/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                            CTTTTTTGCAGCAGGTGGGTAAAGCCGAGGATGCCATTGAGCGGGGTGCGGATTTCGTGGCTCATGTTGGCCAGGAACTCCGACTTGATGCGGCTGGCCTCAAGTGCCTCTTTGCGCGCCATGTCCAGCTCGATGTTCTGGATTTCGATGGTTTCCAGGTTCTGGCGCACGTCTTCGGTGGCCTGGTCGATGCTGTGTTGCAGTTCTTCGTGGGCATTCTGCAGGGTTTCGGCCATGCGGTTGATGCCGTGGGCCAGTTCGTCCAGCTCATGGCTGCCCATGG  >  NC_002947/1843853‑1844135

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