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%	48.4 / NA	16	Y185H (TAT→CAT)	yihS	sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/9);  total (7/9)

TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGG  >  NZ_CP009273/4064479‑4064647                                                                                      |                                                                                    tGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACa                                                                                 <  2:379312/90‑1 (MQ=255)     tattTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCa                                                                             >  1:402741/1‑90 (MQ=255)          cTGGCGACGTCGTGGATAAACACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaaa                                                                        <  1:177991/90‑1 (MQ=255)            ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          >  1:255659/1‑86 (MQ=255)            ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  2:255659/86‑1 (MQ=255)            ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                                      >  1:232061/1‑90 (MQ=255)                               aCGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTg                                                   >  1:359678/1‑90 (MQ=255)                                 ggAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGc                                                 >  1:275932/1‑90 (MQ=255)                                  gAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCa                                                >  1:163542/1‑90 (MQ=255)                                       cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  2:402741/90‑1 (MQ=255)                                                      aTCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTc                            >  1:426677/1‑90 (MQ=255)                                                           gCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGttt                       <  2:232061/90‑1 (MQ=255)                                                                      tCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCt            <  1:47319/90‑1 (MQ=255)                                                                        gtgagtgaCGTCATGAACGATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTc          <  1:347705/90‑1 (MQ=255)                                                                             tgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTccc     <  1:270262/90‑1 (MQ=255)                                                                                cGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAgg  <  2:474108/90‑1 (MQ=255)                                                                                      |                                                                                    TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGG  >  NZ_CP009273/4064479‑4064647

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