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% 179.0 / NA 56*478* (TAG→TAAglgAglycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/27);  total (29/27)

ATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATA  >  NC_000913/3566465‑3566732
                                                                                                                                       |                                                                                                                                    
aTGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTAt                                                                                                                                   >  1:1193321/1‑139 (MQ=255)
     tATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCaa                                                                                                                              >  1:1716958/1‑139 (MQ=255)
                                        aTCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGa                                                                                           <  1:1881974/139‑1 (MQ=255)
                                          cAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAcc                                                                                         >  1:530797/1‑139 (MQ=255)
                                             cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                      >  1:1775944/1‑139 (MQ=255)
                                                        tAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATc                                                                           <  1:1795040/139‑1 (MQ=255)
                                                                  aaGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATa                                                                 >  1:1257981/1‑139 (MQ=255)
                                                                       cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcg                                                                                                 >  2:628446/1‑102 (MQ=255)
                                                                       cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcg                                                                                                 <  1:628446/102‑1 (MQ=255)
                                                                       cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTg                                                            >  2:32347/1‑139 (MQ=255)
                                                                         tCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGAATTCTAATAGACGTTCCCTGAAAATTATTTGAAGCGATATTAAGCCTCACGTTACACCTTCGCCTTAACCTCCCAGCAAAATTCCATTGCAATAGCCTgtc                                                          >  1:391845/1‑137 (MQ=255)
                                                                                  aGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                                                          <  1:382141/98‑1 (MQ=255)
                                                                                  aGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                                                          >  2:382141/1‑98 (MQ=255)
                                                                                  aGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACaa                                                 >  1:1731828/1‑139 (MQ=255)
                                                                                    cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                               <  2:1731828/139‑1 (MQ=255)
                                                                                       gggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                            <  1:1744322/139‑1 (MQ=255)
                                                                                       gggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                            >  1:211635/1‑139 (MQ=255)
                                                                                         gCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTaaaa                                                                             >  2:1263931/1‑104 (MQ=255)
                                                                                         gCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTaaaa                                                                             <  1:1263931/104‑1 (MQ=255)
                                                                                         gCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                         >  2:1199819/1‑108 (MQ=255)
                                                                                         gCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                         <  1:1199819/108‑1 (MQ=255)
                                                                                            aTGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAg                                       <  1:965322/139‑1 (MQ=255)
                                                                                                 taagTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAg                                  <  2:391845/136‑1 (MQ=255)
                                                                                                   tGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATATCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATGCATTGACATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGa                                <  2:1257981/139‑1 (MQ=255)
                                                                                                       aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     <  1:195765/82‑1 (MQ=255)
                                                                                                       aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     >  2:195765/1‑82 (MQ=255)
                                                                                                       aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     >  1:475456/1‑82 (MQ=255)
                                                                                                       aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     <  2:475456/82‑1 (MQ=255)
                                                                                                       aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTa                                                                                <  2:87020/87‑1 (MQ=255)
                                                                                                       aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTa                                                                                >  1:87020/1‑87 (MQ=255)
                                                                                                         cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCt                                                                                                                >  1:1651764/1‑53 (MQ=255)
                                                                                                         cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCt                                                                                                                <  2:1651764/53‑1 (MQ=255)
                                                                                                         cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  2:1775944/139‑1 (MQ=255)
                                                                                                         cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  1:32347/139‑1 (MQ=255)
                                                                                                            agaaTTAATATAGGCGTTTCCTGAAAATTATTTCAAGAGATAGAAAAGCTCCCGGTACGACTTGGCCGCGACTTGCCAGCTAAAATCCACTGCCATAGTCTGCCTTTGCATAAATGAATACAGTGAAGGACGGGACcac                       >  1:157119/4‑139 (MQ=255)
                                                                                                            aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                            <  2:166529/118‑1 (MQ=255)
                                                                                                            aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                            >  1:166529/1‑118 (MQ=255)
                                                                                                            aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                       >  2:556119/1‑139 (MQ=255)
                                                                                                             gCATTCATATAGGCGTTTCCTGAAAATTATTTCAAg                                                                                                                             <  2:1579804/36‑1 (MQ=255)
                                                                                                             gCATTCATATAGGCGTTTCCTGAAAATTATTTCAAg                                                                                                                             >  1:1579804/1‑36 (MQ=255)
                                                                                                             gCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTg                                                                                       <  2:1374943/74‑1 (MQ=255)
                                                                                                             gCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTg                                                                                       >  1:1374943/1‑74 (MQ=255)
                                                                                                              cATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                         >  2:362144/1‑87 (MQ=255)
                                                                                                              cATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                         <  1:362144/87‑1 (MQ=255)
                                                                                                                      taGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAg             >  2:792240/1‑139 (MQ=255)
                                                                                                                        ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCa           >  1:880713/1‑139 (MQ=255)
                                                                                                                        ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCa           >  2:218232/1‑139 (MQ=255)
                                                                                                                           gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                            <  2:220824/55‑1 (MQ=255)
                                                                                                                           gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                            >  1:220824/1‑55 (MQ=255)
                                                                                                                            tttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGca                                                    <  1:1749422/94‑1 (MQ=255)
                                                                                                                            tttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGca                                                    >  2:1749422/1‑94 (MQ=255)
                                                                                                                              tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAg                                  >  1:1771325/1‑110 (MQ=255)
                                                                                                                              tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAg                                  <  2:1771325/110‑1 (MQ=255)
                                                                                                                               ccTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa    <  2:211635/139‑1 (MQ=255)
                                                                                                                               ccTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa    <  2:1193321/139‑1 (MQ=255)
                                                                                                                                 tGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATa  <  1:792240/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
ATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATA  >  NC_000913/3566465‑3566732

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