Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	2,907,329	Δ1 bp	coding (205/1221 nt)	h16_A2693 →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	2,907,327	0	G	.	100.0%	86.2 / NA	20	coding (203/1221 nt)	h16_A2693	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/11);  total (9/11)

CGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCA  >  AM260479/2907195‑2907435                                                                                                                                     |                                                                                                             cgcAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCCggc                                                                                                            <  2:847218/134‑2 (MQ=255)                   ccgccgGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCt                                                                                         >  2:322436/1‑133 (MQ=255)                                        aaCCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                                   <  2:266740/134‑1 (MQ=255)                                        aaCCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                                   <  1:85838/134‑1 (MQ=255)                                                      gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGc                                                     >  1:613751/1‑134 (MQ=255)                                                       ccGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGcc                                                    <  1:494611/134‑1 (MQ=255)                                                                              gaggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggc                             <  1:322436/134‑1 (MQ=255)                                                                               aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcg                            >  2:282395/1‑134 (MQ=255)                                                                                 gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCt                                         <  1:202877/119‑1 (MQ=255)                                                                                 gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCt                                         >  2:202877/1‑119 (MQ=255)                                                                                 gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcggc                          >  1:611093/1‑134 (MQ=255)                                                                                   gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                                   <  2:275894/91‑1 (MQ=255)                                                                                   gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                                   >  1:275894/1‑91 (MQ=255)                                                                                         cGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                  <  1:783785/134‑1 (MQ=255)                                                                                         cGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                  >  1:431833/1‑134 (MQ=255)                                                                                          gagctgagctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTc                 <  2:431833/134‑1 (MQ=255)                                                                                             ctgagctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCg               >  1:63329/1‑133 (MQ=255)                                                                                             ctgagctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGc              <  2:611093/134‑1 (MQ=255)                                                                                                         gAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGgcgcg                                            >  2:1019096/1‑92 (MQ=255)                                                                                                         gAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGgcgcg                                            <  1:1019096/92‑1 (MQ=255)                                                                                                         gAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCa  <  1:282395/134‑1 (MQ=255)                                                                                                                                     |                                                                                                             CGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCA  >  AM260479/2907195‑2907435

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