Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260480	978,518	A→G	S117G (AGC→GGC)	h16_B0862 →	cytosine deaminase or related metal‑dependent hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260480	978,518	0	A	G	100.0%	95.2 / NA	25	S117G (AGC→GGC)	h16_B0862	cytosine deaminase or related metal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/12);  total (13/12)

ATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAAGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGG  >  AM260480/978401‑978609                                                                                                                      |                                                                                            aTGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGa                                                                             <  1:947012/134‑1 (MQ=255)         gATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGc                                                                                <  2:479925/123‑1 (MQ=255)         gATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGc                                                                                >  1:479925/1‑123 (MQ=255)          aTGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCacctgta                                                                    <  1:77907/134‑1 (MQ=255)                                    ccaccTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCt                                                                                      >  1:1116742/1‑90 (MQ=255)                                    ccaccTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCt                                                                                      <  2:1116742/90‑1 (MQ=255)                                              cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTg                                                                                     >  2:109446/1‑81 (MQ=255)                                              cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTg                                                                                     <  1:109446/81‑1 (MQ=255)                                              cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAg                                                                                  <  1:976340/84‑1 (MQ=255)                                              cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAg                                                                                  >  2:976340/1‑84 (MQ=255)                                              cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                >  1:190190/1‑134 (MQ=255)                                              cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                >  1:146645/1‑134 (MQ=255)                                                   gTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCgcggcg                           >  1:463947/1‑134 (MQ=255)                                                   gTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCgcggcg                           >  2:942934/1‑134 (MQ=255)                                                      cGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCggcggc                         <  2:463947/133‑1 (MQ=255)                                                               gCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                <  1:77087/117‑1 (MQ=255)                                                               gCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                >  2:77087/1‑117 (MQ=255)                                                                         ggCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGc     <  2:146645/134‑1 (MQ=255)                                                                            gAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCgg  >  1:256909/1‑134 (MQ=255)                                                                             aaCCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTcac                    <  1:366786/115‑1 (MQ=255)                                                                             aaCCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTcac                    >  2:366786/1‑115 (MQ=255)                                                                                            tCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCt                        >  1:994889/1‑96 (MQ=255)                                                                                            tCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCt                        >  1:260012/1‑96 (MQ=255)                                                                                            tCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCt                        <  2:260012/96‑1 (MQ=255)                                                                                            tCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCt                        <  2:994889/96‑1 (MQ=255)                                                                                                                      |                                                                                            ATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAAGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGG  >  AM260480/978401‑978609

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