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% 84.8 / NA 27*407* (TAG→TAAwcaLputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/13);  total (14/13)

TAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCG  >  NC_000913/2115778‑2116041
                                                                                                                                 |                                                                                                                                      
tAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGcc                                                                                                                               >  1:1122288/1‑139 (MQ=255)
         aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAg                                                                                                                      <  2:1122288/139‑1 (MQ=255)
              aaCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTGAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 <  2:512224/139‑1 (MQ=255)
              aaCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 <  1:881925/139‑1 (MQ=255)
              aaCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAACAAGCTGGc                                                                                                                 <  2:1269167/139‑1 (MQ=255)
                                   cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  1:1096428/139‑1 (MQ=255)
                                     cgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 >  2:199579/1‑116 (MQ=255)
                                     cgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 <  1:199579/116‑1 (MQ=255)
                                         ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGtt                                                                                      >  2:907601/1‑139 (MQ=255)
                                                            aaGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                                                                   >  1:648775/1‑139 (MQ=255)
                                                             aGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgca                                                                                                                           >  1:155827/1‑82 (MQ=255)
                                                             aGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgca                                                                                                                           <  2:155827/82‑1 (MQ=255)
                                                              ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  1:32605/112‑1 (MQ=255)
                                                              ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            >  2:32605/1‑112 (MQ=255)
                                                                      ccGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcgc                                                         >  1:158401/1‑139 (MQ=255)
                                                                         aGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcagca                                                                                                                        <  1:1262890/73‑1 (MQ=255)
                                                                         aGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcagca                                                                                                                        >  2:1262890/1‑73 (MQ=255)
                                                                                    aGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  2:1259126/90‑1 (MQ=255)
                                                                                    aGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            >  1:1259126/1‑90 (MQ=255)
                                                                                    aGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGa                                           >  2:1197142/1‑139 (MQ=255)
                                                                                           tGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAAtt                                    <  1:1197142/139‑1 (MQ=255)
                                                                                                     gagTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTc                          >  2:798688/1‑139 (MQ=255)
                                                                                                      agTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            >  1:274630/1‑72 (MQ=255)
                                                                                                      agTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  2:274630/72‑1 (MQ=255)
                                                                                                          ttttAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGtt                     >  1:908030/1‑139 (MQ=255)
                                                                                                                aTGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTaa               >  1:463581/1‑139 (MQ=255)
                                                                                                                 tGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTaaa              <  2:865642/139‑1 (MQ=255)
                                                                                                                            ctccttTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcgtcg                               >  2:1153336/1‑111 (MQ=255)
                                                                                                                            ctccttTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcgtcg                               <  1:1153336/111‑1 (MQ=255)
                                                                                                                             tccttTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCg  >  2:445287/1‑139 (MQ=255)
                                                                                                                             tccttTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCg  >  1:492874/1‑139 (MQ=255)
                                                                                                                              ccttTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                         >  2:1237118/1‑35 (MQ=255)
                                                                                                                              ccttTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                         <  1:1237118/35‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCG  >  NC_000913/2115778‑2116041

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