Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	637,111	+G	intergenic (‑164/+218)	adk ← / ← kdsB	adenylate kinase/3‑Deoxy‑manno‑octulosonate cytidylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	637,107	1	.	G	100.0%	137.7 / NA	36	intergenic (‑160/+222)	adk/kdsB	adenylate kinase/3‑Deoxy‑manno‑octulosonate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/21);  total (15/21)

GGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAG  >  AM260479/636978‑637214                                                                                                                                   |                                                                                                            ggCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gggg                                                                                                         >  2:981993/1‑134 (MQ=255) ggCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gggg                                                                                                         >  1:373901/1‑134 (MQ=255)  gCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gggg                                                                                                         >  2:503711/1‑133 (MQ=255)          tCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggcggc                                                                                                 >  2:1047006/1‑134 (MQ=255)                   gggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        <  1:655371/134‑1 (MQ=255)                   gggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        <  1:503711/134‑1 (MQ=255)                         cGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                  <  1:919861/134‑1 (MQ=255)                                  ttGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATc                                                                         >  1:911640/1‑134 (MQ=255)                                    gACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCt                                                                        <  2:225011/133‑1 (MQ=255)                                            aCTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTccc                                                               <  1:611506/134‑1 (MQ=255)                                                 cTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCTCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                          <  1:981993/134‑1 (MQ=255)                                                 cTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                          <  2:373901/134‑1 (MQ=255)                                                    cGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATc                                                                         <  1:355676/116‑1 (MQ=255)                                                    cGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATc                                                                         <  1:1016964/116‑1 (MQ=255)                                                    cGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATc                                                                         >  2:355676/1‑116 (MQ=255)                                                    cGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATc                                                                         >  2:1016964/1‑116 (MQ=255)                                                          tGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgcg                                                                                              >  1:374771/1‑89 (MQ=255)                                                          tGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgcg                                                                                              <  2:374771/89‑1 (MQ=255)                                                             tGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                  >  2:845960/1‑98 (MQ=255)                                                             tGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                  <  1:845960/98‑1 (MQ=255)                                                                     gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        <  2:774185/84‑1 (MQ=255)                                                                     gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        >  1:774185/1‑84 (MQ=255)                                                                     gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGc                                      >  1:735611/1‑134 (MQ=255)                                                                              cAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                          <  1:1144152/105‑1 (MQ=255)                                                                              cAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                          >  2:1144152/1‑105 (MQ=255)                                                                                   cgGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                        <  2:735611/134‑1 (MQ=255)                                                                                   cgGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                        <  1:1047006/134‑1 (MQ=255)                                                                                        aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        >  2:1084449/1‑65 (MQ=255)                                                                                        aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        <  1:1084449/65‑1 (MQ=255)                                                                                        aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTa                                                                           <  1:679149/78‑1 (MQ=255)                                                                                        aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTa                                                                           >  2:679149/1‑78 (MQ=255)                                                                                            cTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggcggc                <  1:12767/133‑1 (MQ=255)                                                                                               aGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        <  2:1180179/58‑1 (MQ=39)                                                                                               aGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                        >  1:1180179/1‑58 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTTGCCGGAGGCACGCTTGCGCATgccggc                        >  2:119141/1‑111 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                        <  1:119141/111‑1 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAg  >  1:217089/1‑133 (MQ=255)                                                                                                                         gcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                   >  2:1160956/1‑101 (MQ=255)                                                                                                                         gcggcTGCCAGGGGGATGGCGGCGCAGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                   <  1:1160956/101‑1 (MQ=255)                                                                                                                                   |                                                                                                            GGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAG  >  AM260479/636978‑637214

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