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%	166.2 / NA	43	intergenic (‑160/+222)	adk/kdsB	adenylate kinase/3‑Deoxy‑manno‑octulosonate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (21/22);  total (21/22)

GGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTATG  >  AM260479/636978‑637221                                                                                                                                   |                                                                                                                   ggCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gggg                                                                                                                >  2:162794/1‑134 (MQ=255)            aGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgc                                                                                                      >  2:756241/1‑134 (MQ=255)            aGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgc                                                                                                      >  2:22198/1‑134 (MQ=255)                   gggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                               <  2:979605/134‑1 (MQ=255)                   gggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                               <  1:756241/134‑1 (MQ=255)                   gggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                               <  1:162794/134‑1 (MQ=255)                        aCGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCAACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGAc                                                                                          >  2:394257/1‑134 (MQ=255)                         cGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         <  1:1036649/134‑1 (MQ=255)                            cGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGc                                                                                      >  2:1160643/1‑134 (MQ=255)                                                           gATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                               >  1:1045926/1‑94 (MQ=255)                                                           gATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                               <  2:1045926/94‑1 (MQ=255)                                                                     gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTAt                                                                           <  1:686629/104‑1 (MQ=255)                                                                     gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTAt                                                                           >  2:686629/1‑104 (MQ=255)                                                                       cTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCtt                                           >  1:515762/1‑134 (MQ=255)                                                                          cGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTgcg                                        <  1:1160643/134‑1 (MQ=255)                                                                              cAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTccc                                                                      >  1:309916/1‑100 (MQ=255)                                                                              cAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTccc                                                                      <  2:309916/100‑1 (MQ=255)                                                                               aGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgc                                   >  2:925475/1‑134 (MQ=255)                                                                               aGCGCGGGCACTGCTCAGGCGTACCAACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATg                                    <  1:394257/133‑1 (MQ=255)                                                                                gcgcGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         >  2:1107065/1‑79 (MQ=255)                                                                                gcgcGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         <  1:1107065/79‑1 (MQ=255)                                                                                gcgcGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgcc                                  >  1:772423/1‑134 (MQ=255)                                                                                   cgGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                               <  2:304428/134‑1 (MQ=255)                                                                                      gCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGccggccgg                            <  2:515762/134‑1 (MQ=255)                                                                                        aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                          <  2:438943/134‑1 (MQ=255)                                                                                            cTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggccgtatgcg                       <  1:250854/133‑7 (MQ=255)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         <  2:320384/55‑1 (MQ=39)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         >  1:320384/1‑55 (MQ=39)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                                 >  1:900742/1‑79 (MQ=255)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                                 <  2:900742/79‑1 (MQ=255)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCtt                                                           <  2:299073/85‑1 (MQ=255)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCtt                                                           >  1:299073/1‑85 (MQ=255)                                                                                                        cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGc            >  2:868652/1‑132 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCg                                                      <  1:814184/88‑1 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCg                                                      >  2:814184/1‑88 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                               <  1:769456/111‑1 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                               >  2:769456/1‑111 (MQ=255)                                                                                                          cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAg         >  2:634136/1‑133 (MQ=255)                                                                                                              tACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTa    >  2:737534/1‑134 (MQ=255)                                                                                                               aCGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACAcc                                                                                              <  2:697459/43‑1 (MQ=38)                                                                                                               aCGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACAcc                                                                                              >  1:697459/1‑43 (MQ=38)                                                                                                               aCGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         <  2:476534/48‑1 (MQ=38)                                                                                                               aCGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                         >  1:476534/1‑48 (MQ=38)                                                                                                                cggccggcGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTATg  <  2:772423/134‑1 (MQ=255)                                                                                                                                   |                                                                                                                   GGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTATG  >  AM260479/636978‑637221

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