Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 587,628 A→G F336S (TTC→TCC)  nrfB ← cyclic di‑3',5'‑guanylate‑activated glycosyltransferase NrfB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273587,6280AG100.0% 86.0 / NA 24F336S (TTC→TCC) nrfBcyclic di‑3',5'‑guanylate‑activated glycosyltransferase NrfB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/7);  total (17/7)

TTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTT  >  NZ_CP009273/587551‑587714
                                                                             |                                                                                      
ttGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCGGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCgat                                                                           >  2:121826/1‑91 (MQ=255)
  gCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCgatga                                                                         >  2:227127/1‑91 (MQ=255)
      cTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCa                                                                     >  1:129077/1‑91 (MQ=255)
              tCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGAtt                                                             >  1:756102/1‑91 (MQ=255)
                gcgcCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGAtttt                                                           >  1:643414/1‑91 (MQ=255)
                  gcCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTg                                                         >  1:153549/1‑91 (MQ=255)
                            gtagtTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCa                                               <  2:839357/91‑1 (MQ=255)
                                         aGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAACGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTAt                                  >  2:276068/1‑91 (MQ=255)
                                              ggAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCgggg                             >  2:448122/1‑91 (MQ=255)
                                              ggAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCgggg                             >  2:539442/1‑91 (MQ=255)
                                               gAGGTCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGa                            >  2:342899/1‑91 (MQ=255)
                                                   tCCATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAata                        >  2:508430/1‑91 (MQ=255)
                                                    ccATTTATGGGTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAatat                       >  2:855628/1‑91 (MQ=255)
                                                              gTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGcgcacgc             <  1:795936/91‑1 (MQ=255)
                                                              gTTTTAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGcgcacgc             <  2:570932/91‑1 (MQ=255)
                                                                tttAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGaaaa                                       >  1:8640/1‑63 (MQ=255)
                                                                tttAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGaaaa                                       <  2:8640/63‑1 (MQ=255)
                                                                tttAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAg           >  2:804871/1‑91 (MQ=255)
                                                                 ttAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGa          >  1:468706/1‑91 (MQ=255)
                                                                 ttAAAGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGa          <  2:726853/91‑1 (MQ=255)
                                                                     aGCCTTGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCAt      <  2:672176/91‑1 (MQ=255)
                                                                         ttGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTc                                            <  2:181068/49‑1 (MQ=255)
                                                                         ttGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTc                                            >  1:181068/1‑49 (MQ=255)
                                                                         ttGGGAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGttt  >  1:117612/1‑91 (MQ=255)
                                                                             |                                                                                      
TTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTT  >  NZ_CP009273/587551‑587714

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