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% 132.6 / NA 41*407* (TAG→TAA) wcaLputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/16);  total (25/16)

TGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGT  >  NC_000913/2115771‑2116039
                                                                                                                                        |                                                                                                                                    
tgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGTGTTCTCCTTta                                                                                                                                    >  2:1219433/1‑139 (MQ=255)
tgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTta                                                                                                                                    >  2:1041202/1‑139 (MQ=255)
     tGTAGTCGTTAATGTTAACGCTTTGTCGGCCGGGAAGCGCGCGGGCGAAAGGTGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAg                                                                                                                               >  1:1086851/1‑139 (MQ=255)
       tAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGcc                                                                                                                             >  2:220088/1‑139 (MQ=255)
              tAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcagca                                                                                                                      >  2:359422/1‑139 (MQ=255)
                           ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcagc                                                                                                                       <  1:1207054/125‑1 (MQ=255)
                           ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTATCCTTTATAAAGCCTgcagc                                                                                                                       >  2:1207054/1‑125 (MQ=255)
                             gTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGa                                                                                                             <  2:1657651/133‑1 (MQ=255)
                             gTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGa                                                                                                             >  1:1657651/1‑133 (MQ=255)
                             gTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                       >  2:1164126/1‑139 (MQ=255)
                                          cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                          <  1:220088/139‑1 (MQ=255)
                                                         ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                                                           <  1:359422/139‑1 (MQ=255)
                                                           tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGtt                                                                         >  2:419353/1‑139 (MQ=255)
                                                           tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGtt                                                                         >  1:1824302/1‑139 (MQ=255)
                                                                  gaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                                                                  >  1:1396015/1‑139 (MQ=255)
                                                                    aGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttttt                                                                >  2:426892/1‑139 (MQ=255)
                                                                     ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                                                                <  2:1651071/106‑1 (MQ=255)
                                                                     ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                                                                >  1:1651071/1‑106 (MQ=255)
                                                                       cAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgc                                                             >  2:1558398/1‑139 (MQ=255)
                                                                          gcgcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtata                                                                                                                                  >  2:274752/1‑67 (MQ=255)
                                                                          gcgcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtata                                                                                                                                  <  1:274752/67‑1 (MQ=255)
                                                                            gcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGtt                                                                                                          <  1:212032/89‑1 (MQ=255)
                                                                            gcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGtt                                                                                                          >  2:212032/1‑89 (MQ=255)
                                                                                     gCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtat                                                                                                                                   <  1:261379/55‑1 (MQ=255)
                                                                                     gCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtat                                                                                                                                   >  2:261379/1‑55 (MQ=255)
                                                                                     gCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTgacg                                               >  1:275409/1‑139 (MQ=255)
                                                                                                cGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCaa                                    >  1:1209617/1‑139 (MQ=255)
                                                                                                  tGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                       >  1:195362/1‑70 (MQ=255)
                                                                                                  tGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                       <  2:195362/70‑1 (MQ=255)
                                                                                                    cGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                          <  2:149059/81‑1 (MQ=255)
                                                                                                    cGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                          >  1:149059/1‑81 (MQ=255)
                                                                                                     gTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                       <  1:713930/67‑1 (MQ=255)
                                                                                                     gTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                       >  2:713930/1‑67 (MQ=255)
                                                                                                     gTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcgt                               >  2:389470/1‑139 (MQ=255)
                                                                                                     gTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcgt                               <  1:1164126/139‑1 (MQ=255)
                                                                                                                tttttAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGt                    <  1:426892/139‑1 (MQ=255)
                                                                                                                 ttttAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGtt                   >  1:943554/1‑139 (MQ=255)
                                                                                                                       aTGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgc                                                                                         <  1:776854/63‑1 (MQ=255)
                                                                                                                       aTGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgc                                                                                         >  2:776854/1‑63 (MQ=255)
                                                                                                                                gTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATCGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCa    <  1:1219433/139‑1 (MQ=255)
                                                                                                                                  tctcCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGt  <  2:1824302/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGT  >  NC_000913/2115771‑2116039

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: ‑