Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,053,497	T→C	E894E (GAA→GAG)	torS ←	sensory histidine kinase TorS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,053,497	0	T	C	100.0%	86.6 / NA	26	E894E (GAA→GAG)	torS	sensory histidine kinase TorS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/12);  total (14/12)

TGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTCTCG  >  NC_000913/1053368‑1053623                                                                                                                                  |                                                                                                                               tGCCAGTCCGTTAAGATGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATctct                                                                                                                             <  1:243013/133‑1 (MQ=255)  gCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATctctt                                                                                                                            >  1:66452/1‑133 (MQ=255)       tCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACAATGCTTCCAGAGCGGCTACACTGCGAGTAACCTCTTCgtgt                                                                                                                       >  2:263288/1‑133 (MQ=255)                 tGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAgggggg                                                                                                             <  1:481603/133‑1 (MQ=255)                     taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCATGTAATCTCTTCGTCTGGAAGAGGGGCAGa                                                                                                         >  1:172611/1‑133 (MQ=255)                               cACTATTCTACTTGTCCTTAACTGTCTGCTATCGATTCAAAGCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGgctgc                                                                                               <  2:172611/133‑1 (MQ=255)                                         gttgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGtt                                                                                     >  2:480367/1‑133 (MQ=255)                                              tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                                                >  2:69451/1‑133 (MQ=255)                                              tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                                                >  2:426897/1‑133 (MQ=255)                                               cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                                               >  1:156265/1‑133 (MQ=255)                                                       gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctgg                                                                       >  1:369582/1‑133 (MQ=255)                                                            tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                                  >  2:509841/1‑133 (MQ=255)                                                                  ttCAAATCGCGTTCAGGTCTTTCTTATGTAACGATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                                            <  2:206862/133‑1 (MQ=255)                                                                      aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                                        >  2:586740/1‑133 (MQ=255)                                                                                      ttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAAc                                        <  2:587591/133‑1 (MQ=255)                                                                                      ttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAAc                                        <  2:369582/133‑1 (MQ=255)                                                                                           aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtt                                   <  1:69451/133‑1 (MQ=255)                                                                                            tGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                                  >  1:495959/1‑133 (MQ=255)                                                                                                 ccATGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc                             >  1:573283/1‑133 (MQ=255)                                                                                                   aTGCTTCCAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg                           <  1:426897/133‑1 (MQ=255)                                                                                                         ccAGAGCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcgc                     >  2:8821/1‑133 (MQ=255)                                                                                                            gagCGGCAACACTGCGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGc                  <  2:66452/133‑1 (MQ=255)                                                                                                              gCGGCAACGCTGCGTGTACTCTATTCGTGTCGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCAc                <  1:263288/133‑1 (MQ=255)                                                                                                                       cTGCGTGTCATCTCTGCGCGTGCAAGGGGGTCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAttt       <  1:167953/133‑1 (MQ=255)                                                                                                                          cGTGTAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGGGCAACTGCGTTTTAGCTGATGCGCGGCACGTTTTATTTtct    >  2:167953/1‑133 (MQ=255)                                                                                                                            tgtAATCTCTTCGTGTGGGAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTCTCg  <  1:268036/133‑1 (MQ=255)                                                                                                                            tgtAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTCTCg  <  1:480367/133‑1 (MQ=255)                                                                                                                            tgtAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTCTCg  <  1:136128/133‑1 (MQ=255)                                                                                                                            tgtAATCTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTCTCg  >  1:12333/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                               TGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTCTCG  >  NC_000913/1053368‑1053623

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