Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,115,907 C→T *407* (TAG→TAA wcaL ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,115,9070CT100.0% 33.9 / NA 12*407* (TAG→TAAwcaLputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/6);  total (6/6)

GGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAG  >  NC_000913/2115772‑2116020
                                                                                                                                       |                                                                                                                 
gggTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtat                                                                                                                <  1:841051/139‑1 (MQ=255)
                                           cgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                                                     >  1:166881/1‑139 (MQ=255)
                                                                  aaGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                                              >  2:459941/1‑139 (MQ=255)
                                                                      cAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTg                                                                               <  1:823899/102‑1 (MQ=255)
                                                                      cAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTg                                                                               >  2:823899/1‑102 (MQ=255)
                                                                         gcgcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcg                                       >  2:571285/1‑139 (MQ=255)
                                                                                  ttGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTga                              <  1:571285/139‑1 (MQ=255)
                                                                                                   cGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcg             <  2:166881/139‑1 (MQ=255)
                                                                                                    gTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAg                                                                                         >  1:244138/1‑62 (MQ=255)
                                                                                                    gTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAg                                                                                         <  2:244138/62‑1 (MQ=255)
                                                                                                            agTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTcc    >  1:540027/1‑139 (MQ=255)
                                                                                                              ttttttAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAg  <  2:540027/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                 
GGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAG  >  NC_000913/2115772‑2116020

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: