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%	128.7 / NA	30	coding (203/1221 nt)	h16_A2693	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/15);  total (15/15)

CCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCG  >  AM260479/2907222‑2907456                                                                                                          |                                                                                                                                  ccGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTCCGAGGCCTACCt                                                                                                              <  2:537598/124‑1 (MQ=255) ccGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTCCGAGGCCTACCt                                                                                                              >  1:537598/1‑124 (MQ=255)   gcgGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTAccc                                                                                                  >  1:1101359/1‑134 (MQ=255)             aaCCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                                                        <  1:213493/134‑1 (MQ=255)                     ctgctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGt                                                                                >  1:118273/1‑134 (MQ=255)                        ctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                                              <  2:769768/133‑1 (MQ=255)                        ctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                                              <  2:118273/133‑1 (MQ=255)                           gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGcc                                                                                                                   <  2:107529/93‑1 (MQ=255)                           gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGcc                                                                                                                   >  1:107529/1‑93 (MQ=255)                           gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGc                                                                          >  2:1150240/1‑134 (MQ=255)                                 atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCCggcg                                                                                                                                >  2:22521/1‑74 (MQ=255)                                 atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCCggcg                                                                                                                                <  1:22521/76‑3 (MQ=255)                                                    aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGaa                                                                                                           <  1:1089421/76‑1 (MQ=255)                                                    aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGaa                                                                                                           >  2:1089421/1‑76 (MQ=255)                                                    aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcg                                                 >  2:1026422/1‑134 (MQ=255)                                                      gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGACCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcggc                                               >  1:512743/1‑134 (MQ=255)                                                        gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGt                                             >  2:1050647/1‑134 (MQ=255)                                                              cGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                                       <  1:531909/134‑1 (MQ=255)                                                                agctgagctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg                                                >  1:427803/1‑123 (MQ=255)                                                                agctgagctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg                                                <  2:427803/123‑1 (MQ=255)                                                                          cTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGcc                                                                                                                   <  1:804116/46‑1 (MQ=12)                                                                          cTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGcc                                                                                                                   >  2:804116/1‑46 (MQ=12)                                                                             ggAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAgg                           <  1:1150240/131‑1 (MQ=255)                                                                                        aagaaCAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGc              <  2:1101359/131‑1 (MQ=255)                                                                                         agaaCAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc                                        >  1:177368/3‑107 (MQ=255)                                                                                         agaaCAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc                                        <  2:177368/105‑1 (MQ=255)                                                                                         agaaCAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCa             >  2:198746/3‑134 (MQ=255)                                                                                              aGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAccgcc        <  1:642558/134‑1 (MQ=255)                                                                                              aGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAccgcc        <  1:111672/134‑1 (MQ=255)                                                                                                  aTTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                                       <  2:122857/99‑1 (MQ=255)                                                                                                  aTTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                                       >  1:122857/1‑99 (MQ=255)                                                                                                    tATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCg  >  2:964483/1‑134 (MQ=255)                                                                                                          |                                                                                                                                  CCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCG  >  AM260479/2907222‑2907456

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