Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,264,970 G→A *419* (TAG→TAA hemA → glutamyl tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,264,9700GA100.0% 131.9 / NA 41*419* (TAG→TAAhemAglutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/20);  total (21/20)

AAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGA  >  NC_000913/1264834‑1265103
                                                                                                                                        |                                                                                                                                     
aaGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACa                                                                                                                                     <  1:1093032/139‑1 (MQ=255)
  gCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGt                                                                                                                                   <  2:1000451/139‑1 (MQ=255)
            aGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTAc                                                                                                                                 >  2:431008/1‑131 (MQ=255)
            aGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTAc                                                                                                                                 <  1:431008/131‑1 (MQ=255)
              gATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCt                                                                                                                       >  2:1561483/1‑139 (MQ=255)
                   ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCtttttt                                                                                                                  >  1:1302777/1‑139 (MQ=255)
                       tGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              <  1:331352/139‑1 (MQ=255)
                                                       aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  2:575555‑M1/1‑112 (MQ=255)
                                                       aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  1:575555‑M1/116‑5 (MQ=255)
                                                        aCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactg                                                                             >  2:1263417‑M1/1‑111 (MQ=255)
                                                                  ttCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaa                                                                   >  1:1097295‑M1/1‑101 (MQ=255)
                                                                    cAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  2:834782‑M1/1‑99 (MQ=255)
                                                                    cAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  1:834782‑M1/103‑5 (MQ=255)
                                                                       cAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGg                                                                                                        >  1:860734‑M1/1‑96 (MQ=255)
                                                                       cAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGg                                                                                                        <  2:860734‑M1/97‑2 (MQ=255)
                                                                          gccgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataag                                                           >  1:1006946‑M1/1‑93 (MQ=255)
                                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                <  1:745887/82‑1 (MQ=255)
                                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                >  2:745887/1‑82 (MQ=255)
                                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagc                                                                 <  2:392506‑M1/129‑41 (MQ=255)
                                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagc                                                                 >  1:392506‑M1/1‑89 (MQ=255)
                                                                                                  cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagccc                                   >  1:88755‑M1/1‑69 (MQ=255)
                                                                                                        aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                             <  2:88755‑M1/139‑77 (MQ=255)
                                                                                                               tGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              >  1:1592093/1‑51 (MQ=255)
                                                                                                               tGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              <  2:1592093/51‑1 (MQ=255)
                                                                                                                 cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgt                            >  2:616445‑M1/1‑54 (MQ=255)
                                                                                                                 cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgt                            <  1:616445‑M1/131‑78 (MQ=255)
                                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  2:1613941‑M1/1‑53 (MQ=255)
                                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  1:1613941‑M1/57‑5 (MQ=255)
                                                                                                                    gACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaat                 <  1:1561483‑M1/139‑89 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaaca                                                                                <  1:831427‑M1/74‑26 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaaca                                                                                >  2:831427‑M1/1‑49 (MQ=255)
                                                                                                                               gCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttccgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatggtgaaatcctgctggagc      <  2:1006946‑M1/139‑100 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaatatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgtt  >  1:730848‑M1/1‑36 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagc                                                          <  2:1262838‑M1/83‑48 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagc                                                          >  1:1262838‑M1/1‑36 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  1:567337‑M1/40‑5 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  2:567337‑M1/1‑36 (MQ=255)
                                                                                                                                    aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgc            >  2:744942‑M1/1‑35 (MQ=255)
                                                                                                                                    aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgc            <  1:744942‑M1/128‑94 (MQ=255)
                                                                                                                                    aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                                         >  1:673002‑M1/1‑35 (MQ=255)
                                                                                                                                    aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                                         <  2:673002‑M1/83‑49 (MQ=255)
                                                                                                                                        |                                                                                                                                     
AAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGA  >  NC_000913/1264834‑1265103

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