Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,259,173 A→G T2A (ACC→GCC)  hemA → glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,259,1730AG100.0% 78.7 / NA 23T2A (ACC→GCC) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/8);  total (15/8)

GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258
                                                                                       |                                                                                     
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     >  2:740981/1‑90 (MQ=255)
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     >  1:878030/1‑90 (MQ=255)
           tttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACt                                                                          <  1:643170/90‑1 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                    >  2:3038/1‑90 (MQ=255)
                                  atgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCAc                                                   >  2:579705/1‑90 (MQ=255)
                                     atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg                                                >  2:249043/1‑90 (MQ=255)
                                        cAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTAt                                             >  1:869306/1‑90 (MQ=255)
                                           gCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGc                                          <  1:813536/90‑1 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc                                                                         <  1:665837/56‑1 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc                                                                         >  2:665837/1‑56 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGc                                       >  2:857913/1‑90 (MQ=255)
                                                              gTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGtt                       >  1:711107/1‑90 (MQ=255)
                                                                   tattatTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGc                  <  1:348200/90‑1 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACgg                                                      >  1:304478/1‑51 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACgg                                                      <  2:304478/51‑1 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga                                   >  2:908256/1‑70 (MQ=255)
                                                                      tattTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga                                   <  1:908256/70‑1 (MQ=255)
                                                                       attTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGa              >  2:280829/1‑90 (MQ=255)
                                                                        tttCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGAt             <  1:74484/90‑1 (MQ=255)
                                                                           cccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAg          >  1:691812/1‑90 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         >  2:471907/1‑90 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         >  1:52541/1‑90 (MQ=255)
                                                                                   cATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCa  <  1:92728/90‑1 (MQ=255)
                                                                                       |                                                                                     
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258

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