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%	77.0 / NA	21	T249A (ACG→GCG)	h16_A0595	Uncharacterized Fe‑S protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (11/10);  total (11/10)

GCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGT‑CAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAG  >  AM260479/628541‑628776                                                                                                                  |                                                                                                                            gCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGAt                                                                                                        <  2:528908/134‑1 (MQ=255)     gCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGgcg                                                                                                    <  1:554579/134‑1 (MQ=255)      cTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGgcgc                                                                                                   >  1:844677/1‑134 (MQ=255)         gCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGgcgcg                                                                                                  <  2:799461/132‑1 (MQ=255)                                ggTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCgg                                                                         <  2:844677/134‑1 (MQ=255)                                            gcgcgcGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTAc                                                                                                                   <  2:761250/80‑1 (MQ=255)                                            gcgcgcGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTAc                                                                                                                   >  1:761250/1‑80 (MQ=255)                                                   gAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGc                                                      >  1:500463/1‑134 (MQ=255)                                                                    ggCCCTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcg                                     >  1:708077/1‑134 (MQ=255)                                                                       ccTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGT‑c                                                                                                                            >  1:804656/1‑43 (MQ=38)                                                                       ccTTGTGCTCGATGGTCAGGTAGGAGATGC‑AGCGGCGCGCGT‑c                                                                                                                            <  2:804656/43‑1 (MQ=38)                                                                                       cAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAg                                                                            <  2:650020/76‑1 (MQ=255)                                                                                       cAGGTAGGAGATGC‑AGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAg                                                                            >  1:650020/1‑76 (MQ=255)                                                                                                         ggcgCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGt                                                         >  1:680432/1‑78 (MQ=255)                                                                                                         ggcgCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGt                                                         <  2:680432/78‑1 (MQ=255)                                                                                                         ggcgCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAg  >  1:646511/1‑133 (MQ=255)                                                                                                         ggcgCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAg  >  1:215373/1‑133 (MQ=255)                                                                                                            gcgcgTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCa                                                           <  1:778467/73‑1 (MQ=255)                                                                                                            gcgcgTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCa                                                           <  1:933797/73‑1 (MQ=255)                                                                                                            gcgcgTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCa                                                           >  2:778467/1‑73 (MQ=255)                                                                                                            gcgcgTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCa                                                           >  2:933797/1‑73 (MQ=255)                                                                                                            gcgcgTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAg  >  1:705698/1‑130 (MQ=255)                                                                                                            gcgcgTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAg  <  2:705698/130‑1 (MQ=255)                                                                                                                  |                                                                                                                            GCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGT‑CAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAG  >  AM260479/628541‑628776

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