Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,843,140	G→A	Q566* (CAG→TAG)	gacS ←	sensor protein GacS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,843,140	0	G	A	100.0%	84.6 / NA	27	Q566* (CAG→TAG)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/14);  total (13/14)

TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCG  >  NC_002947/1842993‑1843282                                                                                                                                                    |                                                                                                                                               tCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTAGAGCTGATGATGATAATCGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGAACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTaca                                                                                                                                              <  1:764998/150‑1 (MQ=255)    ggCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTacagca                                                                                                                                           >  2:344817/1‑150 (MQ=255)          gCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGt                                                                                                                                     <  1:68371/150‑1 (MQ=255)             cTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGt                                                                                                                                  <  1:572882/150‑1 (MQ=255)                       gCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCGCCTAAAGAAACGAGGCGTAGGAGTCGaa                                                                                                                        >  2:458344/1‑150 (MQ=255)                                  gCTCGCGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGc                                                                                                             >  1:297640/1‑150 (MQ=255)                                       aGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAgg                                                                                                        <  1:180217/150‑1 (MQ=255)                                        gctgctGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGc                                                                                                       <  1:332989/150‑1 (MQ=255)                                           gctgATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTa                                                                                                                                 >  2:178926/1‑121 (MQ=255)                                           gctgATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTa                                                                                                                                 <  1:178926/121‑1 (MQ=255)                                            ctgATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCa                                                                                                   >  2:484036/1‑150 (MQ=255)                                                tgatAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCaa                                                                                                                                          <  1:162040/107‑1 (MQ=255)                                                tgatAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCaa                                                                                                                                          >  2:162040/1‑107 (MQ=255)                                                          ccAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTc                                                                                                                                    <  1:209277/103‑1 (MQ=255)                                                          ccAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTc                                                                                                                                    >  2:209277/1‑103 (MQ=255)                                                                  tCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAAGTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGt                                                                             >  1:334482/1‑150 (MQ=255)                                                                                aCCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCa                                                                <  2:819077/149‑1 (MQ=255)                                                                                aCCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCa                                                                >  1:819077/1‑149 (MQ=255)                                                                                      gTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTg                                                         <  2:431559/150‑1 (MQ=255)                                                                                            aCCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATc                                                   >  2:412346/1‑150 (MQ=255)                                                                                                 gCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACaa                                              >  1:221008/1‑150 (MQ=255)                                                                                                     cgcgAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATggc                                          >  2:213137/1‑150 (MQ=255)                                                                                                           atCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcg                                    >  2:597687/2‑150 (MQ=255)                                                                                                              gAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcggcg                                 >  2:81983/1‑150 (MQ=255)                                                                                                                         cggccTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGACTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATg                      <  1:424577/150‑1 (MQ=255)                                                                                                                                 ctggcGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGt              <  1:344817/150‑1 (MQ=255)                                                                                                                                    gcggcCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTGCGTGGCCATCGCCAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGc           <  2:162334/150‑1 (MQ=255)                                                                                                                                            cACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTc   <  1:484036/150‑1 (MQ=255)                                                                                                                                             aCCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGTTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCg  <  2:221008/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCG  >  NC_002947/1842993‑1843282

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