Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,568,113 G→A R406C (CGT→TGT)  glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,568,1130GA100.0% 48.3 / NA 17R406C (CGT→TGT) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/8);  total (9/8)

GGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACGACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATC  >  NC_000913/3567899‑3568312
                                                                                                                                                                                                                      |                                                                                                                                                                                                       
gggAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAACCCTGCATTATCGCTCCTGTTCATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCAt                                                                                                                                                                     >  1:328825/1‑251 (MQ=255)
gggAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTCATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCAt                                                                                                                                                                     >  1:328816/1‑251 (MQ=255)
                    cncgAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAAcaca                                                                                                                                                   >  1:255674/1‑249 (MQ=255)
                             cAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCTGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTGCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTcc                                                                                                                                                                                            <  1:208300/199‑1 (MQ=255)
                                        gCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGTGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAgcgg                                                                                                                             >  1:20901/1‑251 (MQ=255)
                                                       ggTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAATCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGCACGan                                                                                                              <  1:191500/251‑2 (MQ=255)
                                                                              naGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGACTCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCgg                                                                                                                                                            >  1:61992/2‑182 (MQ=255)
                                                                                    gccgGTTTGAAGCAGCGGGAACATCGCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTATCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACATGCATCTTCCTCTGCGTTTTCACCAAGCACCAGGCCTTCCGGAATAACACAAGCACGATTGATGACTCAGCGGCGCAGTCGGCACGAGCGACCTACCCATACTTCCGGTAACAATAc                                                                                 <  1:313116/251‑1 (MQ=255)
                                                                                                          aTCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCNCCANNNNnccnnnnnn                                                                                                                                                                >  2:205241/1‑143 (MQ=255)
                                                                                                          aTCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATAcggcn                                                                             <  1:205241/233‑2 (MQ=255)
                                                                                                                         tGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTCCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTGCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCAcc                                                                                                                                                                       <  1:217034/128‑1 (MQ=255)
                                                                                                                                ccTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGCACGAGCGAcc                                                                                                         <  1:20337/183‑1 (MQ=255)
                                                                                                                                                        cccTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGTACGAGCGACCTAc                                                                                                      >  1:315175/1‑162 (MQ=255)
                                                                                                                                                                             gcgcgTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTCCTCTGCGTTTTCACCACTCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCCCCACCGAACCGGAGATc  >  1:228553/1‑241 (MQ=255)
                                                                                                                                                                                 gTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACAACGTGCATCTTACTCTGGGTTTTCACCAATCACCAt                                                                                                                                                                     >  1:136284/1‑74 (MQ=255)
                                                                                                                                                                                                ccTTCATCTGAACGATAGAAACAACGTTCATCTTCCTCTGCGTTTTCCCCAATCACCATGCCTTCCGGAATAACACAAGCa                                                                                                                                               <  1:147158/81‑1 (MQ=255)
                                                                                                                                                                                                         gNNCGATAGNNACAACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACgg                                                                                                                    <  2:20901/99‑1 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                       
GGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTTATGCCCTAACTTCCGTAGCATTTCGCGCGTTACCAGCACGATGCCTTCTTCTGAACGATAGAAACGACGTGCATCTTCCTCTGCGTTTTCACCAATCACCATGCCTTCCGGAATAACACAAGCACGATCGATGACGCAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATC  >  NC_000913/3567899‑3568312

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