Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,256,301	G→A	S117S (TCC→TCT)	cyuA ←	putative L‑cysteine desulfidase CyuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,256,301	0	G	A	100.0%	38.6 / NA	15	S117S (TCC→TCT)	cyuA	putative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/6);  total (9/6)

GCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAC  >  NC_000913/3256056‑3256515                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                       gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTAAGCGGCGACTCTTTCTCGCCCTCTGCCACCCACGCCTGCTGGGTAAAAAACACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGCTTTACTCGCAAGGTTCCTGGATCTTAACAGAGAc                                                                                                                                                                                                                   >  1:328942/1‑251 (MQ=255)         gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACAATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCccg                                                                                                                                                                                                           >  1:98410/1‑251 (MQ=255)                                       cagcGTCGTTCGGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCCCACACGCCTGCTGGGTAAACACCACCCCATCGGGCGTCTCGATATGCACAATGGTGGTATGCCCGCCGACGATGGTGACACACGCCCACGTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTCCATCGCAAGGTTCCTGGATCTTAACAGGGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTg                                                                                                                                                                             <  1:199577/251‑1 (MQ=255)                                                                            tacgcgtcctcagactcacgcGCCTGCTTGGTAAACCCCACCCCTTCGTTCGTCTCGATAAGCACAATGTTGATATGCCCGCCGACGATGGTGACACTCGCCCACTTCTCACCGTTCCTGACTTTTGCGCGTGAAAAGAGGAGTTCATCTCAATGTCCCTGGATCTTAACAGAGCCTTGCCCCGCCGCCAGCAGTGCTTTGTCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCg                                                                                                                                          <  1:251478/229‑1 (MQ=255)                                                                              cTCGCCCTCTACCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGGTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACGTTGGCGCTTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTcn                                                                                                                                                                                                              <  1:71197/179‑2 (MQ=255)                                                                                                                           gTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAAGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCACGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCAc                                                                                                                          >  1:202218/1‑218 (MQ=255)                                                                                                                                                      aTGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATAGCAAGGTTCCTGGATCTAAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCACCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCACGTGCCGGGAACgggg                                                               >  1:159087/1‑248 (MQ=255)                                                                                                                                                               gacgaTGTTGACACACGCCCACTTCTCACCGTTCCAGTCTTTGGCGCTTGAGAAGAGGATGTCATCGCAAGGGTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATc                                                                                                                                                                                     <  1:310886/123‑1 (MQ=255)                                                                                                                                                               gacgaTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGAc                                                     >  1:299576/1‑251 (MQ=255)                                                                                                                                                                                cccACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCATTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCACTCGGCAGCCCCACCATTCCCGTGCCGGGAACGGGGACGCCCAGACCGTTCTTCATCAGAttt                                    >  1:11308/1‑251 (MQ=255)                                                                                                                                                                                         cACCGTTCCAGACGTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTAAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCGTTCAGCACTTCCATCCCGGCGTTGGCATTTCCACCTACCGCCCCCAGCGCCGCCGCAATCTGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAAc                            <  1:45164/250‑1 (MQ=255)                                                                                                                                                                                                          ncgcgTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTc          >  1:87335/2‑251 (MQ=255)                                                                                                                                                                                                                 agaagaGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAAACTTTCCCCGCCGCCAACAGTGCTTTGGCATCGCCAATTGCATGCGCTGTTGCGTATTTAAGCACTTCCAGCCCGGCGTTGGCATTTACACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCCTTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCCTCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCgg    >  1:177993/1‑250 (MQ=255)                                                                                                                                                                                                                  gaagaGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGan  <  1:284998/251‑2 (MQ=255)                                                                                                                                                                                                                                     ggTTCNNGGNNNTTNNCAGAGACTTTCCCCGCCGCCAGCNGTGCTTTGGCATCGGCAATTGCCTGCGCTNTNGNGTCTNNCAGNACTTCCAGCCNGGCGTTGGCATTTCCACCTAACGCCCCCAgcgc                                                                                                          >  2:255493/1‑128 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                       GCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAC  >  NC_000913/3256056‑3256515

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