Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,612,881	C→A	C124F (TGC→TTC)	glsB ←	glutaminase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,612,881	0	C	A	100.0%	56.2 / NA	18	C124F (TGC→TTC)	glsB	glutaminase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/7);  total (11/7)

CAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGCAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGTTGCCAGATTTCCTCTTCGGAGTAATGACGCATGGCGACAACGAGACTCAGCACTTTGGAAATAGACTGAATGGAAAAACGTTCTTGCGCGTCTCCGGCC  >  NC_000913/1612637‑1613090                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                  cAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGAAACATATCGAAGacca                                                                                                                                                                                                             <  2:137368/251‑1 (MQ=255)           caGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCAcc                                                                                                                                                                                                   >  1:89539/1‑251 (MQ=255)                       agagCGCAGTAATGCAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGCGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATgc                                                                                                                                                                                        >  2:36474/1‑249 (MQ=255)                        gagCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAAc                                                                                                                                                                                      >  1:179598/1‑251 (MQ=255)                        gagCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTAGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGaa                                                                                                                                                                                       >  1:95705/1‑250 (MQ=255)                            gcAGTAATGAAAGTAGGTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCTATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCATCGCACCGGCATTAATGAACGGAt                                                                                                                                                                                  <  2:68899/251‑1 (MQ=255)                                  aTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTAcgcg                                                                                                                                                                            <  1:137778/251‑1 (MQ=255)                                                                       tgGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAAtt                                                                                                                                       >  1:17090/1‑251 (MQ=255)                                                                         gAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTg                                                                                                                                      >  1:152891/1‑250 (MQ=255)                                                                         gAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGAAACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTg                                                                                                                                      >  1:152879/1‑250 (MQ=255)                                                                               gCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACgggg                                                                                                                               >  2:157739/1‑249 (MQ=255)                                                                                       aCTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGAc                                                                                                                       <  1:219483/251‑1 (MQ=255)                                                                                         ttCATCAGCCCGGCGATAGCCCCCTTTCGCGCGGCATGTCCCAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATACGACACACCGCTTCAGCCGCGCCCGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCCACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGGATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGCGATCCAGACgg                                                                                                                     >  2:101288/1‑251 (MQ=255)                                                                                                                 tttCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGTTGcc                                                                                                 >  1:18655/1‑247 (MQ=255)                                                                                                                                   aaTTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGTTGCCAGATTTCCGATTCGGAGTAATg                                                                           <  1:69510/251‑1 (MQ=255)                                                                                                                                                               aGGAAATATCAGACACTCCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGGTGCCAGATTTCCTCTGCGGAGTAATGACTCATGGCGACAACGAGACTCAGCACt                                               <  1:30252/251‑1 (MQ=255)                                                                                                                                                                       tCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGTTGCCAGATTTCCTCTTCGGAGTAATGACGCATGGCGACAACGCGACTCAGCACTTTGGAAAt                                       >  1:56392/1‑251 (MQ=255)                                                                                                                                                                                                             aCGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGAAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGTTGCCAGATTTCCTCTTCGGAGTAATGACGCATGGCGACAACGAGACTCAGCACTTTGGAAATAGACTGAATGGAAAAACGTTCTTGCGCGTCTCCGGcc  <  1:21076/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                  CAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTAATCGCCCTTGCAACATATCGCAGACCACCAGCGCACCGGCATTAATGAACGGATTACGCGGTATACCCTGCTCCATTTCCAGTTGCACTAAGGAATTGAACGGTGATCCAGACGGATCTTTGCCGACGCGTTGCCAGATTTCCTCTTCGGAGTAATGACGCATGGCGACAACGAGACTCAGCACTTTGGAAATAGACTGAATGGAAAAACGTTCTTGCGCGTCTCCGGCC  >  NC_000913/1612637‑1613090

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑