Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,646,214	(G)6→7	coding (362/846 nt)	sseA →	3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,646,208	1	.	G	100.0%	99.4 / NA	27	L119R (CTG→CGG)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (13/14);  total (13/14)

TTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGC  >  NZ_CP009273/2646075‑2646340                                                                                                                                       |                                                                                                                                     ttCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggggg                                                                                                                                 >  1:667378/1‑139 (MQ=255)  tCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggg                                                                                                                                >  1:18652/1‑139 (MQ=255)     ccGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtg                                                                                                                             <  1:758526/139‑4 (MQ=255)          gCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                                                          <  1:246531/138‑1 (MQ=255)           cGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTg                                                                                                                        <  2:365954/139‑1 (MQ=255)                                  cAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                                                                 <  2:118691/139‑1 (MQ=255)                                  cAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                                                                 <  1:449757/139‑1 (MQ=255)                                              cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                                                       >  2:807478/1‑137 (MQ=255)                                                 aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                                                                    <  1:280897/137‑1 (MQ=255)                                                     tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                                                              >  1:843277/1‑139 (MQ=255)                                                               aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                                                                    >  1:704280/1‑139 (MQ=255)                                                                         ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                                                             <  2:357754/136‑1 (MQ=255)                                                                                                 aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTg                                    >  1:10000/1‑137 (MQ=255)                                                                                                    ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGcc                               <  2:18652/139‑1 (MQ=255)                                                                                                      gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgt                             >  2:439726/1‑139 (MQ=255)                                                                                                       cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtg                            >  1:799026/1‑139 (MQ=255)                                                                                                       cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtg                            >  1:253420/1‑139 (MQ=255)                                                                                                       cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtg                            >  2:473227/1‑139 (MQ=255)                                                                                                                  tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAAc                 <  2:833758/139‑1 (MQ=255)                                                                                                                    tAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCg               <  1:749774/139‑1 (MQ=255)                                                                                                                      gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                                                                    <  2:241686/68‑1 (MQ=255)                                                                                                                      gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                                                                    >  1:241686/1‑68 (MQ=255)                                                                                                                           aGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtg                            >  2:252256/1‑119 (MQ=255)                                                                                                                           aGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtg                            <  1:252256/119‑1 (MQ=255)                                                                                                                           aGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtat         >  1:487520/1‑138 (MQ=255)                                                                                                                            gtgtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtatta       >  2:274908/1‑139 (MQ=255)                                                                                                                              gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                                                      <  1:741776/58‑1 (MQ=255)                                                                                                                              gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                                                      >  2:741776/1‑58 (MQ=255)                                                                                                                                 gATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGc  <  2:843277/139‑1 (MQ=255)                                                                                                                                 gATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGc  <  1:136215/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     TTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGC  >  NZ_CP009273/2646075‑2646340

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