Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	1,464,433	A→G	Q39Q (CAA→CAG)	h16_A1354 →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	1,464,433	0	A	G	100.0%	98.4 / NA	26	Q39Q (CAA→CAG)	h16_A1354	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/15);  total (11/15)

GGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGC  >  AM260479/1464308‑1464558                                                                                                                              |                                                                                                                              ggAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcg                                                                                                                      >  2:781416‑M1/1‑132 (MQ=255)           tGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggc                                                                                                            <  2:48122‑M1/134‑13 (MQ=255)               cGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcg                                                                                                        <  2:779017‑M1/134‑17 (MQ=255)                              tcatcaGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                                                          >  1:491925‑M1/1‑103 (MQ=255)                                atcaGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGccgct                                                                                                                          <  1:586684/99‑1 (MQ=255)                                atcaGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGccgct                                                                                                                          >  2:586684/1‑99 (MQ=255)                                atcaGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccag                                                                                       <  2:146716‑M1/134‑34 (MQ=255)                                atcaGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                                                          >  2:629468‑M1/1‑101 (MQ=255)                                    gCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagca                                                                                   <  2:808880‑M1/134‑38 (MQ=255)                                     cTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcac                                                                                  >  2:685341‑M1/1‑96 (MQ=255)                                           ggCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcatt                                                                            <  1:629468‑M1/134‑45 (MQ=255)                                                   gTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagc                                                                    >  2:660099‑M1/1‑82 (MQ=255)                                                              ccgccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaata                                                           <  1:660099‑M1/133‑62 (MQ=255)                                                                   acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattggg                                                                         >  2:239733‑M1/1‑67 (MQ=255)                                                                   acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattggg                                                                         <  1:239733‑M1/114‑48 (MQ=255)                                                                   acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggc                                                                                                      <  1:784227‑M1/85‑19 (MQ=255)                                                                   acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggc                                                                                                      >  2:784227‑M1/1‑67 (MQ=255)                                                                       gAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaat                                                 <  1:685341‑M1/134‑72 (MQ=255)                                                                                gCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccg                                          <  2:491925‑M1/132‑79 (MQ=255)                                                                                   gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgccc                                     >  2:440774‑M1/1‑51 (MQ=255)                                                                                     ctctGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccgg                                   <  2:131608‑M1/134‑86 (MQ=255)                                                                                        tGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggt                                  <  1:781416‑M1/132‑87 (MQ=255)                                                                                          gCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgac                              <  1:440774‑M1/134‑91 (MQ=255)                                                                                                  gCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctga                      >  1:561985‑M1/1‑36 (MQ=255)                                                                                                                      ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcaggca  <  2:561985‑M1/134‑119 (MQ=255)                                                                                                                      ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcaggca  >  2:885858‑M1/1‑16 (MQ=255)                                                                                                                              |                                                                                                                              GGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGC  >  AM260479/1464308‑1464558

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