Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	628,654	+C	coding (744/1206 nt)	h16_A0595 ←	Uncharacterized Fe‑S protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	628,653	1	.	C	100.0%	65.3 / NA	18	T249A (ACG→GCG)	h16_A0595	Uncharacterized Fe‑S protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (10/8);  total (10/8)

GCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGT‑CAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAGGA  >  AM260479/628520‑628778                                                                                                                                       |                                                                                                                              gcgAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGT‑c                                                                                                                              <  1:505779/134‑1 (MQ=255)                gggCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGAGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGgcgc                                                                                                                <  1:310304/134‑1 (MQ=255)                       ccAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGt                                                                                                                       >  1:16289/1‑120 (MQ=255)                       ccAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGt                                                                                                                       <  2:16289/120‑1 (MQ=255)                          gCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGgcg                                                                                                      <  2:758802/134‑1 (MQ=255)                               cAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGt                                                                                                 >  1:149117/1‑134 (MQ=255)                                                         gCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTgg                                                                        >  1:493146/1‑133 (MQ=255)                                                            cATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGt                                                                    >  2:966524/1‑134 (MQ=255)                                                                         aaCTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGa                                                                                                           >  1:668353/1‑82 (MQ=255)                                                                         aaCTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGa                                                                                                           <  2:668353/82‑1 (MQ=255)                                                                                                gtgCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcggccggc                                 <  1:898596/133‑1 (MQ=255)                                                                                                gtgCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcggccggc                                 <  2:493146/133‑1 (MQ=255)                                                                                                  gCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGa                                                                                                           <  2:211804/57‑1 (MQ=38)                                                                                                  gCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGa                                                                                                           >  1:211804/1‑57 (MQ=38)                                                                                                  gCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGc                              >  1:250695/1‑134 (MQ=255)                                                                                                            cAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGa                    >  2:993564/1‑134 (MQ=255)                                                                                                               gTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGAc                   >  1:427986/1‑132 (MQ=255)                                                                                                                         gcaGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGc        >  1:582569/3‑134 (MQ=255)                                                                                                                               gcgcgcGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAGga  <  2:250695/134‑1 (MQ=255)                                                                                                                                       |                                                                                                                              GCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGT‑CAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAGGA  >  AM260479/628520‑628778

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