Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,090,333 A→G D80G (GAC→GGC)  hisIE → bifunctional phosphoribosyl‑AMP cyclohydrolase/phosphoribosyl‑ATP diphosphatase HisIE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,090,3330AG100.0% 67.7 / NA 22D80G (GAC→GGC) hisIEbifunctional phosphoribosyl‑AMP cyclohydrolase/phosphoribosyl‑ATP diphosphatase HisIE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/8);  total (14/8)

CGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACC  >  NZ_CP009273/2090258‑2090418
                                                                           |                                                                                     
cgcgCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGaca                                                                         <  2:937700/90‑1 (MQ=255)
 gcgcACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGCGTATTGCCCCGGGCTGCGACAACGacac                                                                        >  2:311604/1‑90 (MQ=255)
 gcgcACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGacac                                                                        >  1:563693/1‑90 (MQ=255)
              cGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTg                                                           >  1:712880/1‑90 (MQ=255)
              cGACTGGGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTg                                                           >  1:1040202/1‑90 (MQ=255)
                aCTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGc                                                         >  2:803475/1‑90 (MQ=255)
                aCTGTGGACCAAAGGCGAAACGGCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGc                                                         >  2:464243/1‑90 (MQ=255)
                  tgtgGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGa                                                       >  2:267543/1‑90 (MQ=255)
                   gtgGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGaa                                                      >  1:566924/1‑90 (MQ=255)
                          aaaGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATc                                               >  2:259260/1‑90 (MQ=255)
                               cGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGccc                                          <  1:259260/90‑1 (MQ=255)
                                             aTTTCCTCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCAAACCCCATCGGCCCGCCCTGCCACAAAg                            <  2:1040202/90‑1 (MQ=255)
                                                   tCAACGTAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACca                      <  2:712880/90‑1 (MQ=255)
                                                         tAGTGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGtt                >  2:965590/1‑90 (MQ=255)
                                                            tGAGTATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCt             <  1:192488/90‑1 (MQ=255)
                                                                tATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGAGAATCCCATCGGccc                                          <  2:332775/57‑1 (MQ=255)
                                                                tATTGCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGAGAATCCCATCGGccc                                          >  1:332775/1‑57 (MQ=255)
                                                                    gCCCCGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGac     >  2:101041/1‑90 (MQ=255)
                                                                     ccccGGGCTGCGACAACGCCACGTTACGGGTGCGGGCGAACCCCATCGCCCCGCCCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGaca    <  2:366034/90‑1 (MQ=255)
                                                                       ccGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGac     >  2:690286/1‑87 (MQ=255)
                                                                       ccGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGac     <  1:690286/87‑1 (MQ=255)
                                                                       ccGGGCTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACAcc  >  1:366034/1‑90 (MQ=255)
                                                                           |                                                                                     
CGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACC  >  NZ_CP009273/2090258‑2090418

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