Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,411,884 A→G K78E (AAG→GAG)  yhdJ → DNA adenine methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,411,8840AG100.0% 175.9 / NA 45K78E (AAG→GAG) yhdJDNA adenine methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/31);  total (14/31)

GCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAGGCGAAAAAACACTACGGCTCCA  >  NC_000913/3411743‑3412028
                                                                                                                                             |                                                                                                                                                
gccgAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGgca                                                                                                                                         <  1:395730/151‑1 (MQ=255)
  cgAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGgcagc                                                                                                                                       >  1:603333/1‑151 (MQ=255)
  cgAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGgcagc                                                                                                                                       >  1:539715/1‑151 (MQ=255)
      aGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGt                                                                                                                                   <  2:629630/151‑1 (MQ=255)
                 gatctTTGCCGACCCAACATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGa                                                                                                                          <  1:400878/149‑1 (MQ=255)
                      ttGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACa                                                                                                                       <  1:69270/147‑1 (MQ=255)
                       tGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTAc                                                                                                                   <  1:191590/150‑1 (MQ=255)
                       tGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACg                                                                                                                  >  1:302222/1‑151 (MQ=255)
                        gCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACgg                                                                                                                 <  1:411962/151‑1 (MQ=255)
                          cGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGa                                                                                                                >  1:135285/1‑150 (MQ=255)
                              ccaccaTATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGaa                                                                                                                         >  1:267641/1‑137 (MQ=255)
                              ccaccaTATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGaa                                                                                                                         <  2:267634/137‑1 (MQ=255)
                                     ataACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCtt                                                                                                    >  2:494647/1‑151 (MQ=255)
                                     ataACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCt                                                                                                     <  1:506886/150‑1 (MQ=255)
                                       aaCATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCttt                                                                                                   <  1:196035/150‑1 (MQ=255)
                                       aaCATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTa                                                                                                  <  2:134392/151‑1 (MQ=255)
                                        aCATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTAt                                                                                                 <  2:128779/151‑1 (MQ=255)
                                             ggTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGAtc                                                                                            <  2:306325/151‑1 (MQ=255)
                                              gTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGAtct                                                                                           >  2:215290/1‑151 (MQ=255)
                                                aaaaaTTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTcc                                                                                         <  2:564547/151‑1 (MQ=255)
                                                        tGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGc                                                                                 >  2:582555/1‑151 (MQ=255)
                                                         gATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCa                                                                                <  2:124579/151‑1 (MQ=255)
                                                                 gATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCtttt                                                                         <  1:63074/150‑1 (MQ=255)
                                                                 gATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCttt                                                                          <  1:206976/149‑1 (MQ=255)
                                                                 gATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCttt                                                                          <  1:622705/149‑1 (MQ=255)
                                                                 gATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCtt                                                                           <  2:118716/148‑1 (MQ=255)
                                                                  aTCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCttttt                                                                        >  1:11937/1‑150 (MQ=255)
                                                                   tCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGc                                                                                    <  2:592673/137‑1 (MQ=255)
                                                                   tCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGc                                                                                    >  1:592688/1‑137 (MQ=255)
                                                                      aaGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCAt                                                                   <  1:401752/151‑1 (MQ=255)
                                                                            ggAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCaaa                                                               >  1:618424/1‑149 (MQ=255)
                                                                               aagaagATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCg                                                          <  2:546736/151‑1 (MQ=255)
                                                                                         gTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCg                                                     <  2:541867/146‑1 (MQ=255)
                                                                                               cGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCatat                                          >  1:90564/1‑151 (MQ=255)
                                                                                                gACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATg                                         <  2:586685/151‑1 (MQ=255)
                                                                                                         tttGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTg                                >  2:387837/1‑151 (MQ=255)
                                                                                                          ttGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTg                                <  2:195038/150‑1 (MQ=255)
                                                                                                            gAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGa                              <  1:64384/150‑1 (MQ=255)
                                                                                                              aGTGATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTg                           <  1:266658/151‑1 (MQ=255)
                                                                                                                 gATTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAg                        <  2:298204/151‑1 (MQ=255)
                                                                                                                  aTTGCAGAGTGCCACCGCGTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAgg                       <  1:162020/151‑1 (MQ=255)
                                                                                                                                   cgTTCTGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAGGCGAAAAAACACTACGGc      <  1:186960/151‑1 (MQ=255)
                                                                                                                                        tGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAGGCGAAAAAACACTACGGCTCCa  >  1:215922/1‑150 (MQ=255)
                                                                                                                                        tGAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAGGCGAAAAAACACTACGGCTCCa  <  2:215917/150‑1 (MQ=255)
                                                                                                                                         gAAAGAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAGGCGAAAAAACACTACGGCTCCa  <  2:4637/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATCGATCTCCAGTGCCGCAAGCTTTTTACCATCAAAAGTCGCATCGTCTGGTCATATGACAGTTCTGGAGTGCAGGCGAAAAAACACTACGGCTCCA  >  NC_000913/3411743‑3412028

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