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

ATGGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTAT  >  AM260479/636976‑637220                                                                                                                                     |                                                                                                                  aTGGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gg                                                                                                                 >  1:925272/1‑134 (MQ=255)            tCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggcggc                                                                                                       >  2:735503/1‑134 (MQ=255)                     gggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                              <  2:190468/134‑1 (MQ=255)                         aaCGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGa                                                                                          >  1:311089/1‑134 (MQ=255)                         aaCGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGa                                                                                          >  1:542557/1‑134 (MQ=255)                          aCGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGAc                                                                                         >  2:554284/1‑134 (MQ=255)                           cGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                        <  2:925272/134‑1 (MQ=255)                           cGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                        <  2:185134/134‑1 (MQ=255)                                    ttGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATc                                                                               >  1:94767/1‑134 (MQ=255)                                                      cGGTCCTGATGCATTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑ggggg                                                                                                              >  2:750318/1‑83 (MQ=255)                                                      cGGTCCTGATGCATTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑ggggg                                                                                                              <  1:750318/84‑2 (MQ=255)                                                       ggTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACggg                                                                                      <  2:827204/108‑1 (MQ=255)                                                       ggTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACggg                                                                                      >  1:827204/1‑108 (MQ=255)                                                       ggTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGcc                                                            <  2:1106726/134‑1 (MQ=255)                                                         tCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCtt                                                          <  2:311089/134‑1 (MQ=255)                                                         tCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCtt                                                          <  2:542557/134‑1 (MQ=255)                                                           cTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTgg                                                        >  2:1098101/1‑134 (MQ=255)                                                           cTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gggggatg                                                                                                           >  1:34132/1‑78 (MQ=255)                                                           cTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑gggggatg                                                                                                           <  2:34132/82‑5 (MQ=255)                                                            tGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgcggc                                                                                                  <  1:814042/91‑1 (MQ=255)                                                            tGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgcggc                                                                                                  >  2:814042/1‑91 (MQ=255)                                                             gATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTa                                                                                 >  1:730715/1‑107 (MQ=255)                                                             gATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTa                                                                                 <  2:730715/107‑1 (MQ=255)                                                                       gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                        >  1:202728/1‑90 (MQ=255)                                                                       gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                        <  2:202728/90‑1 (MQ=255)                                                                       gACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGc                                            >  1:56898/1‑134 (MQ=255)                                                                        aCTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACg                                             <  2:337163/132‑1 (MQ=255)                                                                         cTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCtt                                          >  2:1074853/1‑134 (MQ=255)                                                                                 aGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgc                                  >  2:436441/1‑134 (MQ=255)                                                                                 aGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATg                                   <  1:1098101/133‑1 (MQ=255)                                                                                 aGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTT‑GCCGGAGGCACGCTTGCGCATgcc                                 >  2:352769/1‑134 (MQ=255)                                                                                 aGCGCGGGCACTGCTCAGGCGTA‑CCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgcc                                 <  1:94960/134‑1 (MQ=255)                                                                                     cgGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                              <  1:1074853/134‑1 (MQ=255)                                                                                      gggCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggcgg                                                                                                        <  2:625669/59‑1 (MQ=255)                                                                                      gggCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggcgg                                                                                                        >  1:625669/1‑59 (MQ=255)                                                                                      gggCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTgg                                                        >  1:356338/1‑107 (MQ=255)                                                                                      gggCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTgg                                                        <  2:356338/107‑1 (MQ=255)                                                                                       ggCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggccg                            >  2:917066/1‑134 (MQ=255)                                                                                          aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                         <  2:94767/134‑1 (MQ=255)                                                                                          aCTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCACACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                         <  1:735503/134‑1 (MQ=255)                                                                                                 aGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGctc                    <  2:56898/132‑1 (MQ=255)                                                                                                       aCCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                                <  2:533106/82‑1 (MQ=255)                                                                                                       aCCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCg                                                                >  1:533106/1‑82 (MQ=255)                                                                                                            cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTa                                                                                 <  2:363080/60‑1 (MQ=255)                                                                                                            cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTa                                                                                 >  1:363080/1‑60 (MQ=255)                                                                                                            cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAg        >  1:355123/1‑133 (MQ=255)                                                                                                            cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAg        >  1:400932/1‑133 (MQ=255)                                                                                                              cATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCAt     <  1:554284/134‑1 (MQ=255)                                                                                                                 aCGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTAt  >  2:730732/1‑134 (MQ=255)                                                                                                                          ggcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                              <  2:696187/97‑1 (MQ=255)                                                                                                                          ggcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgccggc                              >  1:696187/1‑97 (MQ=255)                                                                                                                           gcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTAt               >  2:884553/1‑111 (MQ=255)                                                                                                                           gcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTAt               <  1:884553/111‑1 (MQ=255)                                                                                                                            cggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATc              >  2:1062988/1‑111 (MQ=255)                                                                                                                            cggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATc              <  1:1062988/111‑1 (MQ=255)                                                                                                                                     |                                                                                                                  ATGGCGGGAGATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTAT  >  AM260479/636976‑637220

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