Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,538,974	T→C	S204G (AGC→GGC)	cysM ←	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,538,974	0	T	C	100.0%	98.9 / NA	29	S204G (AGC→GGC)	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/15);  total (14/15)

CCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGGAGACAAAATGAGTGATGCG  >  NC_000913/2538842‑2539091                                                                                                                                     |                                                                                                                      cccGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCctct                                                                                                                >  1:241058/1‑140 (MQ=255)     cACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCg                                                                                                            >  1:260184/1‑140 (MQ=255)                  cgcATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGAc                                                                                               >  1:12485/1‑140 (MQ=255)                  cgcATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGAc                                                                                               >  1:242218/1‑140 (MQ=255)                          gtTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCa                                                                                       <  2:242218/140‑1 (MQ=255)                          gtTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCa                                                                                       <  2:241058/140‑1 (MQ=255)                                cGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCa                                                                                       >  2:82588/1‑134 (MQ=255)                                cGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCa                                                                                       <  1:82588/134‑1 (MQ=255)                                   aTCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTgg                                                                              >  1:161368/1‑140 (MQ=255)                                       cgcTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTg                                                                          >  2:233173/1‑140 (MQ=255)                                             tgaATAGCCAGCGCCTCATCCACCAGAGATGCGTTGAAAAGCCCCTGCAGATATGCCGTATGCCATCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGTCCGACAATTGACACCGGTTTTGATTGTTcgcg                                                                    <  2:181447/140‑1 (MQ=255)                                                   tCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATaaa                                                              >  1:86644/1‑140 (MQ=255)                                                      aGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTcg                                                                      <  1:7062/129‑1 (MQ=255)                                                      aGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTcg                                                                      >  2:7062/1‑129 (MQ=255)                                                      aGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                                                           >  1:99173/1‑140 (MQ=255)                                                          ccTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTgag                                                       <  2:260184/140‑1 (MQ=255)                                                           cTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTgaga                                                      >  1:257814/1‑140 (MQ=255)                                                            tCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGAc                                                     <  2:12485/140‑1 (MQ=255)                                                            tCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGAc                                                     <  2:86644/140‑1 (MQ=255)                                                             cATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACg                                                    >  2:207895/1‑140 (MQ=255)                                                             cATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACg                                                    >  1:369428/1‑140 (MQ=255)                                                                ccaccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCg                                                 <  1:32777/140‑1 (MQ=255)                                                                          aGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGc                                       <  2:99173/140‑1 (MQ=255)                                                                          aGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGc                                       <  2:257814/140‑1 (MQ=255)                                                                              ttGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgt                                   <  1:154881/140‑1 (MQ=255)                                                                                             aGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGga                    <  2:307231/140‑1 (MQ=255)                                                                                               atatTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGgaga                  <  1:141359/140‑1 (MQ=255)                                                                                                     cGTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGGAGACAAAAt            >  1:334388/1‑140 (MQ=255)                                                                                                      gTAGGCCAGCGGCGAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGGAGACAAAATg           <  2:369428/140‑1 (MQ=255)                                                                                                               cggcgAATGCCGGGAATGCTGCCGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGGAGACAAAATGAGTGATgcg  <  2:334388/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                      CCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCCCATGCTGGAGACAAAATGAGTGATGCG  >  NC_000913/2538842‑2539091

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