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% 48.2 / NA 16T2A (ACC→GCC) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/10);  total (6/10)

GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258
                                                                                       |                                                                                     
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                     <  1:293494/90‑1 (MQ=255)
           tttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACt                                                                          <  2:337950/90‑1 (MQ=255)
           tttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACt                                                                          <  2:465285/90‑1 (MQ=255)
            ttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc                                                                         <  2:132015/90‑1 (MQ=255)
            ttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc                                                                         <  2:499941/90‑1 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                    >  1:551399/1‑90 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                    >  2:68022/1‑90 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                  <  1:227535/90‑1 (MQ=255)
                     ttACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCGAc                                                                >  1:273743/1‑90 (MQ=255)
                            tCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATaaaa                                                         >  2:393000/1‑90 (MQ=255)
                                                   aCCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGaa                                  >  2:435374/1‑90 (MQ=255)
                                                        aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                             >  1:184198/1‑90 (MQ=255)
                                                                   tattatTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGc                  <  2:184198/90‑1 (MQ=255)
                                                                        tttCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGAt             <  2:67136/90‑1 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc         <  1:68022/90‑1 (MQ=255)
                                                                                   cATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCa  <  1:576157/90‑1 (MQ=255)
                                                                                       |                                                                                     
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA  >  NZ_CP009273/1259086‑1259258

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