Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,566,600 C→T *478* (TAG→TAA glgA ← glycogen synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,566,6000CT100.0% 101.6 / NA 32*478* (TAG→TAAglgAglycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/14);  total (18/14)

GACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATA  >  NC_000913/3566477‑3566732
                                                                                                                           |                                                                                                                                    
gacgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAg                                                                                                                       >  2:29590/1‑139 (MQ=255)
 acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGt                                                                                                                      >  1:548714/1‑139 (MQ=255)
         gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                              >  2:662879/1‑139 (MQ=255)
                       tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgc                                                                                                <  2:152809/139‑1 (MQ=255)
                                           aTAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAAt                                                                            >  1:105310/1‑139 (MQ=255)
                                              gagTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                         <  2:549327/139‑1 (MQ=255)
                                                   cTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGcc                                                                    <  2:172407/139‑1 (MQ=255)
                                                           cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTg                                                            <  2:105310/139‑1 (MQ=255)
                                                              cTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACg                                                         >  2:221812/1‑139 (MQ=255)
                                                                        cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                                                          <  2:714773/96‑1 (MQ=255)
                                                                        cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                                                          >  1:714773/1‑96 (MQ=255)
                                                                        cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                               >  2:87983/1‑139 (MQ=255)
                                                                                       tGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                              >  1:329686/1‑61 (MQ=255)
                                                                                       tGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                              <  2:329686/61‑1 (MQ=255)
                                                                                          aaaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACgg                             <  2:548714/139‑1 (MQ=255)
                                                                                             cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  1:662879/139‑1 (MQ=255)
                                                                                                  cATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAg                     >  2:651766/1‑139 (MQ=255)
                                                                                                            ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  1:386455/124‑1 (MQ=255)
                                                                                                            ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          >  2:386455/1‑124 (MQ=255)
                                                                                                            ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCa           >  2:727129/1‑139 (MQ=255)
                                                                                                             gCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          >  1:285500/1‑123 (MQ=255)
                                                                                                             gCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  2:285500/123‑1 (MQ=255)
                                                                                                               gTTTCCTGAAAATTATTTCAAGCGATATTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGt        >  1:596646/1‑139 (MQ=255)
                                                                                                               gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGt        >  1:236520/1‑139 (MQ=255)
                                                                                                                tttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTc       >  2:98862/1‑139 (MQ=255)
                                                                                                                 ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     >  1:768893/1‑60 (MQ=255)
                                                                                                                 ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     <  2:768893/60‑1 (MQ=255)
                                                                                                                 ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACg                                                         <  2:208560/88‑1 (MQ=255)
                                                                                                                 ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACg                                                         >  1:208560/1‑88 (MQ=255)
                                                                                                                  tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                                  <  2:303976/94‑1 (MQ=255)
                                                                                                                  tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                                  >  1:303976/1‑94 (MQ=255)
                                                                                                                     tGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATa  <  1:221812/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                    
GACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATA  >  NC_000913/3566477‑3566732

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