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,1730AG89.5% 46.0 / ‑0.8 19T2A (ACC→GCC) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (12/5);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGC  >  NZ_CP009273/1259086‑1259261
                                                                                       |                                                                                        
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                        >  2:532499/1‑90 (MQ=255)
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc                                                                                        >  2:265068/1‑90 (MQ=255)
    gcTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCttt                                                                                    >  2:56592/1‑90 (MQ=255)
            ttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTc                                                                            <  1:128688/90‑1 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                       >  2:698602/1‑90 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                       >  2:208891/1‑90 (MQ=255)
                 gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTAt                                                                       >  1:7666/1‑90 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                     <  2:559462/90‑1 (MQ=255)
                   tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa                                                                     <  2:629194/90‑1 (MQ=255)
                              cTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACg                                                          <  1:142216/90‑1 (MQ=255)
                                 tATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCa                                                       <  2:577805/90‑1 (MQ=255)
                                    gatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCt                                                    >  1:203442/1‑90 (MQ=255)
                                              gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGc                                          >  2:28329/1‑90 (MQ=255)
                                                    cccTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAAc                                    <  1:368062/90‑1 (MQ=255)
                                                        aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                                >  1:423769/1‑90 (MQ=255)
                                                        aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                                >  2:449883/1‑90 (MQ=255)
                                                        aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                                >  2:623753/1‑90 (MQ=255)
                                                                            ccGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc            >  2:107081/1‑90 (MQ=255)
                                                                                      gACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGgc  <  1:286976/90‑1 (MQ=255)
                                                                                       |                                                                                        
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGC  >  NZ_CP009273/1259086‑1259261

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