Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	1,464,435	T→C	L40P (CTG→CCG)	h16_A1354 →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	1,464,435	0	T	C	100.0%	108.3 / NA	29	L40P (CTG→CCG)	h16_A1354	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/14);  total (15/14)

TGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGC  >  AM260479/1464318‑1464558                                                                                                                      |                                                                                                                            tGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggc                                                                                                            >  1:522111‑M1/1‑122 (MQ=255)       ccTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggc                                                                                                      <  1:1041529‑M1/134‑19 (MQ=255)          gtgACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcc                                                                                                   >  1:431583‑M1/1‑113 (MQ=255)                  tgtCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgt                                                                                            <  2:431583‑M1/133‑29 (MQ=255)                    tcatcaGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                                                          >  1:500094‑M1/1‑103 (MQ=255)                              gtggGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacgg                                                                                >  1:409012‑M1/1‑93 (MQ=255)                                gggCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcat                                                                             >  2:791344‑M1/1‑91 (MQ=255)                                 ggCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcatt                                                                            <  2:500094‑M1/134‑45 (MQ=255)                                  gCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattg                                                                           >  2:210845‑M1/1‑89 (MQ=255)                                    gAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattggg                                                                         <  1:791344‑M1/134‑48 (MQ=255)                                              agcgtgcCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                              >  1:410620‑M1/4‑78 (MQ=255)                                              agcgtgcCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                              <  2:410620‑M1/101‑27 (MQ=255)                                                       ccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                     <  1:614556‑M1/120‑52 (MQ=255)                                                       ccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                     >  2:614556‑M1/1‑69 (MQ=255)                                                        cacaCGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtgg                                                      <  1:210845‑M1/134‑67 (MQ=255)                                                         acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggc                                                     >  2:822870‑M1/1‑67 (MQ=255)                                                         acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggc                                                     >  1:530119‑M1/1‑67 (MQ=255)                                                                         gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgccc                                     >  1:839508‑M1/1‑51 (MQ=255)                                                                                gCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgac                              <  1:810924‑M1/134‑91 (MQ=255)                                                                                gCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgac                              <  2:522111‑M1/134‑91 (MQ=255)                                                                                          cATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgac                              <  2:1006124‑M1/124‑91 (MQ=255)                                                                                          cATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgac                              >  1:1006124‑M1/1‑34 (MQ=255)                                                                                              aGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatg                 <  1:822870‑M1/133‑104 (MQ=255)                                                                                                           aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcagg    >  2:413890‑M1/1‑17 (MQ=255)                                                                                                            ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcaggca  <  2:839508‑M1/134‑119 (MQ=255)                                                                                                            ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggc                                                     >  1:1064698‑M1/1‑16 (MQ=255)                                                                                                            ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                              <  2:1170912‑M1/42‑27 (MQ=255)                                                                                                            ggTCGACAGCCGCTGGcgcggcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggc                                                     <  2:1064698‑M1/83‑68 (MQ=255)                                                                                                            ggTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGAt                                                                                              >  1:1170912/1‑42 (MQ=14)                                                                                                                      |                                                                                                                            TGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGC  >  AM260479/1464318‑1464558

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