Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 858,839 A→G D43D (GAT→GAC ybiY ← glycyl‑radical enzyme activating protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273858,8390AG100.0% 66.6 / NA 19D43D (GAT→GACybiYglycyl‑radical enzyme activating protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/8);  total (11/8)

ACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAGAAATACGACCGTGCG  >  NZ_CP009273/858752‑858919
                                                                                       |                                                                                
aCCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTc                                                                                >  1:367065/1‑90 (MQ=255)
    ttCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCtgc                                                                            >  1:505466/1‑90 (MQ=255)
      cAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTcc                                                                               <  1:428829/85‑1 (MQ=255)
      cAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTcc                                                                               >  2:428829/1‑85 (MQ=255)
        gcgcgcgcTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCtgcgtgc                                                                        >  2:395412/1‑90 (MQ=255)
           cgcgcTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCggg                                                                     >  2:71176/1‑90 (MQ=255)
                      aCTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTcc                                                          <  2:349233/90‑1 (MQ=255)
                                 ggCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGaca                                               >  2:514307/1‑90 (MQ=255)
                                         cgcACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGc                                       <  2:163296/90‑1 (MQ=255)
                                            acaGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGc                                    >  2:45790/1‑90 (MQ=255)
                                              aGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGccc                                  >  2:500722/1‑90 (MQ=255)
                                                  cGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCg                              >  2:19709/1‑90 (MQ=255)
                                                         ttCCAGACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGcc                       >  1:235125/1‑90 (MQ=255)
                                                              gACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTaa                  <  2:235125/90‑1 (MQ=255)
                                                              gACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTaa                  <  2:505466/90‑1 (MQ=255)
                                                              gACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTaa                  <  1:514307/90‑1 (MQ=255)
                                                              gACACAGTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTaa                  <  1:45790/90‑1 (MQ=255)
                                                                    gTCGTGCGTCATACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAGAAATa            >  2:572204/1‑90 (MQ=255)
                                                                              aTACAGCAGGTCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAGAAATACGACCGTGCg  <  1:290952/90‑1 (MQ=255)
                                                                                       |                                                                                
ACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAGAAATACGACCGTGCG  >  NZ_CP009273/858752‑858919

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