Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,072,783	C→T	R263R (CGG→CGA)	yihU ←	3‑sulfolactaldehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,072,783	0	C	T	100.0%	57.0 / NA	21	R263R (CGG→CGA)	yihU	3‑sulfolactaldehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/5);  total (16/5)

GGTGCCCCAGCCGCGGCAAACATCATGCGCATGGCTTCGCGCTGATCGACCGCCAGCATGGCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACAAAACGGCGGCTTCTGCCG  >  NC_000913/4072542‑4073018                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            ggTGCCCCAGCCGCGGCAAATATCATGCGCATGGCTTCGCGCTGATCGACCGCCAGCATGGCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAggcg                                                                                                                                                                                                                                    <  1:212507/251‑1 (MQ=255)          gCCGCGGCAAATATCATGCGCATGGCTTCGCGCTGATCGACCGCCAGCATGGCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGc                                                                                                                                                                                                                           >  1:247940/1‑251 (MQ=255)                gCAAATATCATGCGCATGGCTTCGCGCTGATCGACAGCCAGCATGGCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTCCTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGCCCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCCGCGGCAtt                                                                                                                                                                                                                     >  1:91504/1‑251 (MQ=255)                            cgcATGGCTTCGCGCTGATCGACCGCCAGCATGGCAAAAACGCCCGATGCGCGCGTAATGTCGAAGATGGTGTCCTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCCGCACTCACACAGACTTGTTCCAGAATGGCGGCACAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCCGCGGCATTGGCACATGCAGc                                                                                                                                                                                                         >  1:317123/1‑251 (MQ=255)                                               cGACCGCCAGCATGGCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATACCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGAACGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAggg                                                                                                                                                                                      >  1:230550/1‑251 (MQ=255)                                                             gCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATcc                                                                                                                                                                        >  1:2182/1‑251 (MQ=255)                                                                                                                        ccTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGAACTGTTCCCGAATGGCGGACCAGTACTGGCGAACGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGGGAAAGATCTCCGCTGAGGACTTTGTTTgg                                                                                                             >  1:219600/1‑251 (MQ=255)                                                                                                                                                             tCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGTCGGCCCCCAGCGGCATTGGc                                                                                                                                                                                                                  <  1:250027/113‑1 (MQ=255)                                                                                                                                                                             gACCTGTTCCAGAATGGCGGACCAGTNCCGGCGNNCGNGACCCGCTGCGCGCGCCTGNCTNNNANNNTNTAGTGNGGNNGCCCCCAGCNGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTANGGGCAAGATNNATCATGAAGGCGGGGGAAAGANATCNNCNNANGCCnttn                                                                                                                  >  2:75456/1‑192 (MQ=255)                                                                                                                                                                                              cGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACtttn                                                         <  1:85058/233‑2 (MQ=255)                                                                                                                                                                                                       ccTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAgg                              >  1:90421/1‑251 (MQ=255)                                                                                                                                                                                                       ccTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACAGGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTGGGCCATGAAGTGGTGAAGTGAcn                                                                                        <  1:261207/193‑2 (MQ=255)                                                                                                                                                                                                              nccgcgaacCGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGCGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACAt                                               >  1:78679/2‑227 (MQ=255)                                                                                                                                                                                                                     nccGCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACAAAAc                >  1:199164/2‑251 (MQ=255)                                                                                                                                                                                                                       ngCTGCGCGCGCCTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACAAAAcgg              >  1:96200/2‑251 (MQ=255)                                                                                                                                                                                                                         cgaagagcgcATGGCTATAAACCTCTAGTGACGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCCGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCCCATCACTTTGACGGCAACATAGAAg                                          >  1:115793/7‑221 (MQ=255)                                                                                                                                                                                                                                  ccTGGCTATAAACCTCTCGTGAGGCGGCCCCCCGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGcc                                                                                                           >  1:276357/1‑147 (MQ=255)                                                                                                                                                                                                                                  ccTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCATTATGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                                                                                                                     >  1:283834/1‑121 (MQ=255)                                                                                                                                                                                                                                   cTGGCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGCTCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACAAAACGGCGGCTTCTGCCg  >  1:164083/1‑251 (MQ=255)                                                                                                                                                                                                                                      gCTATAACCCTCTCGTGGGGCGGCCCCCAGCGGCATTGGCACATGCATCTGGGTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGt                                                                         <  1:169352/177‑1 (MQ=255)                                                                                                                                                                                                                                      gCTATAAACCTCTCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCTGCTGAGGACTNTGTNNNNNCann                                                                                                        >  2:200639/1‑144 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            GGTGCCCCAGCCGCGGCAAACATCATGCGCATGGCTTCGCGCTGATCGACCGCCAGCATGGCAAAACCGCCCGATGCGCGCGTAATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACAAAACGGCGGCTTCTGCCG  >  NC_000913/4072542‑4073018

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