Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,937,802	A→G	D112G (GAC→GGC)	fucK →	L‑fuculokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,937,802	0	A	G	100.0%	114.8 / NA	31	D112G (GAC→GGC)	fucK	L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/14);  total (17/14)

CTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCAC  >  NC_000913/2937672‑2937912                                                                                                                                   |                                                                                                               ctgactgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTg                                                                                                        >  2:525132/1‑139 (MQ=255)  tgactgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGa                                                                                                       >  2:116284/1‑139 (MQ=255)   gactgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGaa                                                                                                      <  2:443767/139‑1 (MQ=255)    actgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAAc                                                                                                     >  1:151669/1‑139 (MQ=255)    actgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAAc                                                                                                     >  1:695285/1‑139 (MQ=255)              acaTCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAAtctc                                                                                           >  2:665577/1‑139 (MQ=255)                     cgGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                    >  2:459501/1‑139 (MQ=255)                      ggTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGc                                                                                   <  1:116284/139‑1 (MQ=255)                             ccGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                            >  2:424529/1‑139 (MQ=255)                                                 gATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAg                                                        <  1:30466/139‑1 (MQ=255)                                                 gATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAg                                                        <  2:151669/139‑1 (MQ=255)                                                      cgcTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCttt                                                   >  1:243463/1‑139 (MQ=255)                                                           tGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGttt                                              >  2:104997/1‑139 (MQ=255)                                                             gTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                            <  2:537906/139‑1 (MQ=255)                                                             gTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                            <  2:110574/139‑1 (MQ=255)                                                                gATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                         <  2:316173/139‑1 (MQ=255)                                                                     gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                            <  1:323801/131‑1 (MQ=255)                                                                     gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                            >  2:323801/1‑131 (MQ=255)                                                                          gCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATaa                               >  1:269189/1‑139 (MQ=255)                                                                          gCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAAAAGCAGGGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGAAGGCTATTTCTGGGTGCGGAGCCTTTAGTTTCAATAGGTGATATaa                               >  1:247754/1‑139 (MQ=255)                                                                              tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                            >  2:326444/1‑122 (MQ=255)                                                                              tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                            <  1:326444/122‑1 (MQ=255)                                                                              tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTg                           <  2:243463/139‑1 (MQ=255)                                                                              tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTg                           >  1:631395/1‑139 (MQ=255)                                                                              tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTg                           >  2:777089/1‑139 (MQ=255)                                                                                     tATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGt                    <  2:695285/139‑1 (MQ=255)                                                                                        ccGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGCGGTTga                 >  2:546079/1‑139 (MQ=255)                                                                                            ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGc                                                                                        >  1:710098/1‑64 (MQ=255)                                                                                            ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGc                                                                                        <  2:710098/64‑1 (MQ=255)                                                                                                       aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcac  <  2:110911/139‑1 (MQ=255)                                                                                                       aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcac  <  1:424529/139‑1 (MQ=255)                                                                                                                                   |                                                                                                               CTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCAC  >  NC_000913/2937672‑2937912

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