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_0009132,937,8020AG100.0% 80.8 / NA 23D112G (GAC→GGC) fucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/13);  total (10/13)

ACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACT  >  NC_000913/2937685‑2937916
                                                                                                                     |                                                                                                                  
acaTCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCa                                                                                                                  >  2:91549/1‑120 (MQ=255)
acaTCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCa                                                                                                                  <  1:91549/120‑1 (MQ=255)
                                          ctctGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAg                                                     <  2:346646/139‑1 (MQ=255)
                                               gTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                <  2:644301/139‑1 (MQ=255)
                                                  gATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                             <  2:61167/139‑1 (MQ=255)
                                                    tAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                           <  1:618511/139‑1 (MQ=255)
                                                      aGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGtt                                         <  1:615649/139‑1 (MQ=255)
                                                       gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGc                                                                                 <  2:132299/98‑1 (MQ=255)
                                                       gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGc                                                                                 >  1:132299/1‑98 (MQ=255)
                                                       gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTta                                        >  2:661478/1‑139 (MQ=255)
                                                            gCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCgg                                                              >  1:699293/1‑112 (MQ=255)
                                                            gCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCgg                                                              <  2:699293/112‑1 (MQ=255)
                                                                                 ttaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGaaaa              >  1:432597/1‑139 (MQ=255)
                                                                                   aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                  <  1:489811/69‑1 (MQ=255)
                                                                                   aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                  >  2:489811/1‑69 (MQ=255)
                                                                                   aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTc                                                                <  2:263261/87‑1 (MQ=255)
                                                                                   aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTc                                                                >  1:263261/1‑87 (MQ=255)
                                                                                   aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAatc            >  1:512679/1‑139 (MQ=255)
                                                                                             gTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAAct  <  2:512679/139‑1 (MQ=255)
                                                                                                  cgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtatat                                     >  2:452337/1‑99 (MQ=255)
                                                                                                  cgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtatat                                     <  1:452337/99‑1 (MQ=255)
                                                                                                              gTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                           >  1:520453/1‑81 (MQ=255)
                                                                                                              gTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                           <  2:520453/81‑1 (MQ=255)
                                                                                                                     |                                                                                                                  
ACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACT  >  NC_000913/2937685‑2937916

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