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

ATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTATGCCAT  >  AM260479/636986‑637225                                                                                                                           |                                                                                                                       aTCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggc                                                                                                               <  2:348063/132‑1 (MQ=255) aTCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggc                                                                                                               >  1:348063/1‑132 (MQ=255)    aGAAGGCGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggc                                                                                                               >  1:449287/1‑129 (MQ=255)    aGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATggc                                                                                                               <  2:449287/129‑1 (MQ=255)    aGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgc                                                                                                          >  2:640513/1‑134 (MQ=255)    aGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgc                                                                                                          >  1:488104/1‑134 (MQ=255)       aggaggGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGgcggcgcggc                                                                                                       <  1:258638/134‑1 (MQ=255)            ggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACAcc                                                                                                  <  1:986072/134‑1 (MQ=255)            ggTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACAcc                                                                                                  <  1:986071/134‑1 (MQ=255)                             aCGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCttt                                                                                 <  1:640513/134‑1 (MQ=255)                                   gACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTcc                                                                           <  2:488104/134‑1 (MQ=255)                                                  tGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                             >  2:993482/1‑101 (MQ=255)                                                  tGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                             <  1:993482/101‑1 (MQ=255)                                                               cTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGa                                                                                               >  2:32275/1‑86 (MQ=255)                                                               cTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGa                                                                                               <  1:32275/86‑1 (MQ=255)                                                               cTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCtt                                               >  2:752611/1‑134 (MQ=255)                                                                      cAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATg                                        >  2:856026/1‑134 (MQ=255)                                                                        gcgcGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                                   <  1:266684/73‑1 (MQ=255)                                                                        gcgcGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCacac                                                                                                   >  2:266684/1‑73 (MQ=255)                                                                        gcgcGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATgcc                                      >  1:1045022/1‑134 (MQ=255)                                                                                       aGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATg                                        >  1:334772/1‑117 (MQ=255)                                                                                       aGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATg                                        <  2:334772/117‑1 (MQ=255)                                                                                          cGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTAt                    <  1:856026/134‑1 (MQ=255)                                                                                          cGTACCCACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTAt                    <  1:752611/134‑1 (MQ=255)                                                                                                cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                             <  2:127429/55‑1 (MQ=39)                                                                                                cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACg                                                                                             >  1:127429/1‑55 (MQ=39)                                                                                                cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTAt                                                                               <  1:43950/69‑1 (MQ=255)                                                                                                cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTAt                                                                               >  2:43950/1‑69 (MQ=255)                                                                                                cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTAt                    <  2:160607/128‑1 (MQ=255)                                                                                                cACGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTAt                    >  1:160607/1‑128 (MQ=255)                                                                                                  cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGccggccgg                                <  2:197168/114‑1 (MQ=255)                                                                                                  cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGccggccgg                                >  1:197168/1‑114 (MQ=255)                                                                                                  cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCaa              <  1:63518/132‑1 (MQ=255)                                                                                                  cGCATACGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCaa              >  2:63518/1‑132 (MQ=255)                                                                                                       aCGGCCGGCGGCGGCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTAt       >  1:708338/1‑134 (MQ=255)                                                                                                            cggcggcggCTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTATGCCAt  <  2:1045022/134‑1 (MQ=255)                                                                                                               cggcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                              >  2:307671/1‑103 (MQ=255)                                                                                                               cggcggcTGCCAGGGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGggc                              <  1:307671/103‑1 (MQ=255)                                                                                                                           |                                                                                                                       ATCAGAAGGAGGGTAACGCCGTAACTTGACGCAGGACTTACTGCGGTCCTGATGCGTTCTGACTGCGTTCAGCGCGGGCACTGCTCAGGCGTACCCACGCATACGGCCGGCGGCGGCTGCCA‑GGGGATGGCGGCGCGGCACACCGGACGGGCTCTATCTTTATCTCCCCACCGAGCCTTGGCCGGAGGCACGCTTGCGCATGCCGGCCGGGCGGCTCTGTATCGGCAAGCATTATGCCAT  >  AM260479/636986‑637225

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