Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779_cfb	3,542,162	A→G	intergenic (+204/‑247)	rhaA → / → rhaD	L‑rhamnose isomerase/rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779_cfb	3,542,162	0	A	G	100.0%	22.9 / NA	11	intergenic (+204/‑247)	rhaA/rhaD	L‑rhamnose isomerase/rhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (4/7);  total (4/7)

GAGAAAGAAATTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCAGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGACAACACCGAAT  >  NC_007779_cfb/3541929‑3542391                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      gagaAAGAAATTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                                                                                            <  1:24261/245‑1 (MQ=255) gagaAAGAAATTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                                                                                            >  2:24261/1‑245 (MQ=255)  agaaagaaATTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                                                                                                                                                                                                                      <  1:764663/250‑1 (MQ=255)   gaaagaaaTTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTGGGCCTGATAAGACGCGACAGCGTCGCACCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCtg                                                                                                                                                                                                                     <  1:896480/250‑2 (MQ=255)   gaaagaaaTCTTGAGTGGCCGCGGGTAACCACTGCCTGAGGCGGCGCGAGCGCTTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                                                                                                                     <  2:531764/250‑1 (MQ=14)    aaagaaaTTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                                                                                                                                                                                                                    <  1:815828/250‑1 (MQ=255)           ttttGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCgt                                                                                                                                                                                                             >  1:1002808/1‑249 (MQ=255)                 gTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGGGCGCCTTATCCGGCCTACGGGGCgtgcatgt                                                                                                                                                                                                       >  2:1396151/1‑243 (MQ=11)                 gTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCgtgcatcc                                                                                                                                                                                                       <  1:869989/250‑8 (MQ=255)                                                                                                                                                                                                                tctgtaaacactgCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCTGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGacaaca        <  1:802032/239‑1 (MQ=11)                                                                                                                                                                                                                   gtaaacactgCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATCGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGACAACACCg     >  1:73322/9‑250 (MQ=11)                                                                                                                                                                                                                      aacactgCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGACAACACCGAAt  >  2:279830/6‑250 (MQ=255)                                                                                                                                                                                                                                    cGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGACAACAcc      >  1:269804/1‑232 (MQ=255)                                                                                                                                                                                                                                    cGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGACAACAcc      <  2:269804/232‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      GAGAAAGAAATTTTGAGTCGCCGCGGGTAAACACTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCCGGCCTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCACAAGTGCCTTATCAGGCCTACAGGTCGGCAATAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTCGTGCATCCGACAACACCGAAT  >  NC_007779_cfb/3541929‑3542391

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