Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,364,361	T→C	intergenic (+101/‑112)	pspE → / → ycjM	thiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,364,361	0	T	C	100.0%	85.7 / NA	25	intergenic (+101/‑112)	pspE/ycjM	thiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/8);  total (17/8)

AGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGT  >  NZ_CP009273/1364275‑1364442                                                                                       |                                                                                  aGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTcccc                                                                                <  1:931279/90‑1 (MQ=255)                         tttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTg                                                        >  2:860782/1‑90 (MQ=255)                            cTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACt                                                     <  1:493919/90‑1 (MQ=255)                                      ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTcccc                                                                                >  2:609471/1‑53 (MQ=255)                                      ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTcccc                                                                                <  1:609471/53‑1 (MQ=255)                                      ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtg                                           >  2:76057/1‑90 (MQ=255)                                      ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtg                                           >  2:162524/1‑90 (MQ=255)                                              aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                   >  1:579397/1‑90 (MQ=255)                                              aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                   >  1:906308/1‑90 (MQ=255)                                              aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                   >  2:1021606/1‑90 (MQ=255)                                              aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                   >  1:127586/1‑90 (MQ=255)                                                ccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                                 >  2:9816/1‑90 (MQ=255)                                                  ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGa                                >  2:900030/1‑89 (MQ=255)                                                  ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGa                                <  1:900030/89‑1 (MQ=255)                                                  ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                               >  2:747896/1‑90 (MQ=255)                                                  ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                               >  2:588972/1‑90 (MQ=255)                                                  ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                               >  2:98773/1‑90 (MQ=255)                                                   cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGTTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                              <  2:1039176/90‑1 (MQ=255)                                                   cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                              <  2:740784/90‑1 (MQ=255)                                                   cATAAATTATTTGCGTCAGCTCACAAATACGCTTTTCCCCTGGT‑AAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAg                             <  1:281252/90‑1 (MQ=255)                                                          tattTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt                       <  1:76057/90‑1 (MQ=255)                                                          tattTGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt                       <  1:162524/90‑1 (MQ=255)                                                           attTGCGTCAGCTCACAAGTACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttt                      >  2:711485/1‑90 (MQ=255)                                                            tttGCGTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt                     >  1:228746/1‑90 (MQ=255)                                                                 gTCAGCTCACAAATACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCg                >  1:577581/1‑90 (MQ=255)                                                                             aTACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtctctc    <  2:906308/90‑1 (MQ=255)                                                                              tACGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCg   >  1:507397/1‑90 (MQ=255)                                                                               aCGCTTTTCCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGt  >  2:630329/1‑90 (MQ=255)                                                                                       |                                                                                  AGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGT  >  NZ_CP009273/1364275‑1364442

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: ‑