Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,064,564	A→G	Y185H (TAT→CAT)	yihS ←	sulfoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,064,564	0	A	G	100.0%	65.3 / NA	21	Y185H (TAT→CAT)	yihS	sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/12);  total (9/12)

ATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGA  >  NZ_CP009273/4064478‑4064648                                                                                       |                                                                                     aTGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAc                                                                                   <  2:570749/90‑1 (MQ=255)         ttCTGGCGGCGTCGTGGATATTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                           <  2:493011/90‑1 (MQ=255)         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                           <  2:445795/90‑1 (MQ=255)                  cgtcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCa                                                                  >  2:266917/1‑90 (MQ=255)                        ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                            <  2:443994/90‑1 (MQ=255)                               cACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATAtt                                                     >  2:601441/1‑90 (MQ=255)                                aCGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTg                                                    >  1:211990/1‑90 (MQ=255)                                         aCGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgccg                                           >  2:186622/1‑90 (MQ=255)                                          cgcgAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgccgc                                          <  1:601397/90‑1 (MQ=255)                                             gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCcgccgc                                       <  2:20183/90‑1 (MQ=255)                                                tCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGt                                    <  2:201891/90‑1 (MQ=255)                                                tCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGt                                    <  2:502002/90‑1 (MQ=255)                                                  gcgcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAc                                  <  2:117995/90‑1 (MQ=255)                                                       aTCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTc                             >  2:435837/1‑90 (MQ=255)                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg                   >  2:410991/1‑90 (MQ=255)                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg                   >  1:89378/1‑90 (MQ=255)                                                                  ttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGa                  <  1:601441/90‑1 (MQ=255)                                                                      gTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGc              >  2:202779/1‑90 (MQ=255)                                                                          tgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCg          >  1:497732/1‑90 (MQ=255)                                                                              tgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTccc      <  2:193203/90‑1 (MQ=255)                                                                                  gTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGa  <  1:410991/90‑1 (MQ=255)                                                                                       |                                                                                     ATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGA  >  NZ_CP009273/4064478‑4064648

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