Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,176,179	T→C	N48S (AAC→AGC)	gatZ ←	tagatose‑1,6‑bisphosphate aldolase 2 subunit GatZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,176,179	0	T	C	100.0%	79.8 / NA	23	N48S (AAC→AGC)	gatZ	tagatose‑1,6‑bisphosphate aldolase 2 subunit GatZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/13);  total (10/13)

AGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGG  >  NC_000913/2176042‑2176305                                                                                                                                          |                                                                                                                               agcagTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGctt                                                                                                                              >  2:711185/1‑140 (MQ=255) agcagTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGctt                                                                                                                              >  1:793998/1‑140 (MQ=255) agcagTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGctt                                                                                                                              >  1:138178/1‑140 (MQ=255)      tttGGTCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGt                                                                                                                         >  2:813430/1‑140 (MQ=255)                             gagaATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGcgcg                                                                                                  <  1:711185/140‑1 (MQ=255)                               gaATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTg                                                                                                <  1:813430/140‑1 (MQ=255)                                          tCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATc                                                                                     >  2:371532/1‑140 (MQ=255)                                          tCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATc                                                                                     >  2:905739/1‑140 (MQ=255)                                          tCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATc                                                                                     >  1:563460/1‑140 (MQ=255)                                                    aCCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAgc                                                                           >  1:217273/1‑140 (MQ=255)                                                              aTCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGa                                                                 <  2:836661/140‑1 (MQ=255)                                                              aTCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGa                                                                 >  1:805335/1‑140 (MQ=255)                                                                     aTCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCaa                                                          <  2:563460/140‑1 (MQ=255)                                                                     aTCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCaa                                                          >  1:549912/1‑140 (MQ=255)                                                                        gTAAAAACAATTCCGCGAAAGTCTGCCGGTGTCATCCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACgg                                                       <  2:948739/140‑1 (MQ=255)                                                                                 aaTTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGa                                              <  1:26851/140‑1 (MQ=255)                                                                                        gAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACt                                       <  2:549912/140‑1 (MQ=255)                                                                                            gTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAAc                                   <  1:905739/140‑1 (MQ=255)                                                                                                 ccGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACa                                           <  2:139413/127‑1 (MQ=255)                                                                                                 ccGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACa                                           >  1:139413/1‑127 (MQ=255)                                                                                                      gtCATTCCGGTATAACCGCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGat                         <  2:28270/140‑1 (MQ=255)                                                                                                                       gCCAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTAt        >  2:424470/1‑140 (MQ=255)                                                                                                                         cAAATTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGc      <  1:371532/140‑1 (MQ=255)                                                                                                                          aaaTTGATTGACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGcc     >  1:570876/1‑140 (MQ=255)                                                                                                                            attgattgACCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCgg   <  2:793998/140‑1 (MQ=255)                                                                                                                             ttgattgaCCTGGCTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCggg  <  2:805335/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                               AGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCCGCTTCGATAACCAACGGATGGGCAGAACAGACTGAACATATGCCGATATGTTCACCAGCTTTATGCCGGG  >  NC_000913/2176042‑2176305

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