Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,579,168	T→C	E339G (GAA→GGA)	ggt ←	gamma‑glutamyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,579,168	0	T	C	100.0%	71.3 / NA	23	E339G (GAA→GGA)	ggt	gamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/15);  total (8/15)

AATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTC  >  NZ_CP009273/3579083‑3579248                                                                                      |                                                                                 aaTTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGa                                                                              <  2:241307/90‑1 (MQ=255)   tttGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGc                                                                            <  1:99950/90‑1 (MQ=255)      gATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGt                                                                         >  1:463550/1‑90 (MQ=255)       aTCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc                                                                        <  2:471132/90‑1 (MQ=255)       aTCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc                                                                        <  2:450908/90‑1 (MQ=255)       aTCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc                                                                        <  2:281270/90‑1 (MQ=255)                     tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  2:57983/1‑90 (MQ=255)                          tAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTcc                                                     >  1:43035/1‑90 (MQ=255)                           aGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCg                                                    >  1:270958/1‑90 (MQ=255)                           aGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCg                                                    >  1:460178/1‑90 (MQ=255)                               tttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCCAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc                                                <  1:502273/90‑1 (MQ=255)                                   ttGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGcc                                            <  2:112452/90‑1 (MQ=255)                                    tGGTCAGCGCCTGCCACGGTACTTTGGCAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCa                                           <  1:36379/90‑1 (MQ=255)                                    tGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCa                                           <  2:411714/90‑1 (MQ=255)                                              cTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcat                                 <  2:195759/90‑1 (MQ=255)                                               tcccACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatc                                >  1:411633/3‑90 (MQ=255)                                               tGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatc                                <  1:594726/90‑1 (MQ=255)                                               tGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatc                                >  2:594726/1‑90 (MQ=255)                                                    cGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatcgcatc                           <  2:494657/90‑1 (MQ=255)                                                        aCTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg                       <  1:57983/90‑1 (MQ=255)                                                           ttGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGACGGCGTAGGCGAATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTg                    >  2:564738/1‑90 (MQ=255)                                                               cAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCaa                <  2:586532/90‑1 (MQ=255)                                                                             ccAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCGATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATttcttc  <  1:548052/90‑1 (MQ=255)                                                                                      |                                                                                 AATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTC  >  NZ_CP009273/3579083‑3579248

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