Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC100.0% 96.3 / NA 29G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/15);  total (14/15)

CCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785605‑3785772
                                                                                  |                                                                                     
ccATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAc                                                                                <  2:63658/90‑1 (MQ=255)
ccATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAc                                                                                <  2:308152/90‑1 (MQ=255)
 cATCCCCAGGAAGCCCGCCAGTTTTTGTTCAAACTCTTTATGCCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAcc                                                                               <  1:3771/90‑1 (MQ=255)
       cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  1:558133/1‑90 (MQ=255)
       cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  2:306735/1‑90 (MQ=255)
                 gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                               >  2:314675/1‑90 (MQ=255)
                 gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                               >  2:21114/1‑90 (MQ=255)
                           gTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgccg                                                     >  1:507860/1‑90 (MQ=255)
                              cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAACCCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  1:240470/90‑1 (MQ=255)
                               aaGCTCTTTATGGCTGCCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCcgccgc                                                 <  1:392720/90‑1 (MQ=255)
                                        aTGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:203389/90‑1 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      >  1:743968/1‑54 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      <  2:743968/54‑1 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:771764/1‑90 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:511178/1‑90 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:49864/1‑90 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:494054/1‑90 (MQ=255)
                                                   gAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCg                             >  2:439014/1‑90 (MQ=255)
                                                        ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:526072/1‑90 (MQ=255)
                                                          gCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  1:436297/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:114867/1‑90 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:435545/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:518329/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:578896/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:187806/1‑90 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCACCCCGAAACCGTGAGAATCCATCCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  1:334222/90‑1 (MQ=255)
                                                                            aGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGcaca    <  2:85583/90‑1 (MQ=255)
                                                                             gCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  1:685614/90‑1 (MQ=255)
                                                                              ccATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaaa  <  1:657519/90‑1 (MQ=255)
                                                                                  |                                                                                     
CCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785605‑3785772

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