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% 207.0 / NA 63*419* (TAG→TAAhemAglutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/29);  total (34/29)

TCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGA  >  NC_000913/1264850‑1265103
                                                                                                                        |                                                                                                                                     
tCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCttt                                                                                                                     >  1:1137688/1‑139 (MQ=255)
   ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCtttttt                                                                                                                  >  1:458518/1‑139 (MQ=255)
     cATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                <  2:458518/139‑1 (MQ=255)
              tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              >  2:389530/1‑132 (MQ=255)
              tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              <  1:389530/132‑1 (MQ=255)
              tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGga                                                                                                       >  1:1498441‑M1/1‑137 (MQ=255)
                 cTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagga                                                                                                    <  2:1498441‑M1/139‑6 (MQ=255)
                    aCCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggagga                                                                                                 >  2:542333‑M1/1‑131 (MQ=255)
                        cTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataa                                                                                             >  1:114550‑M1/1‑127 (MQ=255)
                          tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  1:209367‑M1/129‑5 (MQ=255)
                          tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  2:209367‑M1/1‑125 (MQ=255)
                              ccATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaat                                                                                       >  2:490097‑M1/1‑121 (MQ=255)
                               cATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaata                                                                                      <  2:114550‑M1/139‑20 (MQ=255)
                                aTGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatat                                                                                     <  1:542333‑M1/139‑21 (MQ=255)
                                  gcgcCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  2:1212916‑M1/121‑5 (MQ=255)
                                  gcgcCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  1:1212916‑M1/1‑117 (MQ=255)
                                   cgcCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaa                                                                                  >  1:1503601‑M1/1‑116 (MQ=255)
                                        aCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                             >  2:210132/1‑107 (MQ=255)
                                        aCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                             <  1:210132/107‑1 (MQ=255)
                                                    cAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagc                                                                 >  1:840933‑M1/1‑99 (MQ=255)
                                                       cAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaat                                                              <  2:840933‑M1/139‑44 (MQ=255)
                                                         ggCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaa                                                                                  <  1:120181‑M1/117‑24 (MQ=255)
                                                         ggCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaa                                                                                  >  2:120181‑M1/1‑94 (MQ=255)
                                                          gccgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataag                                                           <  2:1137688‑M1/139‑47 (MQ=255)
                                                          gccgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataag                                                           >  1:1392428‑M1/1‑93 (MQ=255)
                                                            cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              <  2:193779/86‑1 (MQ=255)
                                                            cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              >  1:193779/1‑86 (MQ=255)
                                                            cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatat                                                                                     >  2:634420‑M1/1‑91 (MQ=255)
                                                            cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatat                                                                                     <  1:634420‑M1/111‑21 (MQ=255)
                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTc                                                                                                                        >  1:559471/1‑74 (MQ=255)
                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTc                                                                                                                        <  2:559471/74‑1 (MQ=255)
                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGt                                                                                                          >  2:509732/1‑88 (MQ=255)
                                                              cccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGt                                                                                                          <  1:509732/88‑1 (MQ=255)
                                                                      ggggATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggc                                               <  1:1462932‑M1/139‑59 (MQ=255)
                                                                         gATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGCAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgag                                            >  1:984237‑M1/1‑78 (MQ=255)
                                                                           tAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcg                                          >  2:919844‑M1/1‑76 (MQ=255)
                                                                                aacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagc                                     >  1:362158‑M1/1‑71 (MQ=255)
                                                                                  cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAggg                                                                                                           <  1:719547/67‑1 (MQ=255)
                                                                                  cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAggg                                                                                                           >  2:719547/1‑67 (MQ=255)
                                                                                  cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                          <  2:1260970‑M1/84‑16 (MQ=255)
                                                                                  cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                          >  1:1260970‑M1/1‑69 (MQ=255)
                                                                                  cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                          >  1:1010937‑M1/1‑69 (MQ=255)
                                                                                  cgCCTGAATACTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggcggaataatgg                                                                                          <  2:1010937‑M1/84‑16 (MQ=255)
                                                                                        aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                             >  2:1156198‑M1/1‑63 (MQ=255)
                                                                                        aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                             <  1:262925‑M1/139‑77 (MQ=255)
                                                                                        aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                             >  2:262925‑M1/1‑63 (MQ=255)
                                                                                        aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                             <  1:1156198‑M1/139‑77 (MQ=255)
                                                                                                 cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              >  2:1200712/1‑49 (MQ=255)
                                                                                                 cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              <  1:1200712/49‑1 (MQ=255)
                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                <  1:1028782/46‑1 (MQ=255)
                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                >  2:1028782/1‑46 (MQ=255)
                                                                                                    gACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaat                 >  1:758062‑M1/1‑51 (MQ=255)
                                                                                                    gACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaat                                                              >  2:607610‑M1/1‑51 (MQ=255)
                                                                                                    gACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaat                                                              <  1:607610‑M1/94‑44 (MQ=255)
                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcc               >  1:526246‑M1/1‑49 (MQ=255)
                                                                                                       aGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcct              <  1:919844‑M1/139‑92 (MQ=255)
                                                                                                         ccTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                >  2:535762/1‑39 (MQ=255)
                                                                                                         ccTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                <  1:535762/39‑1 (MQ=255)
                                                                                                            cGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctgg         <  2:526246‑M1/139‑97 (MQ=255)
                                                                                                              ggCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgc                       <  2:1401825‑M1/123‑83 (MQ=255)
                                                                                                              ggCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgc                       >  1:1401825‑M1/1‑41 (MQ=255)
                                                                                                               gCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagc      <  2:362158‑M1/139‑100 (MQ=255)
                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgt  >  2:322678‑M1/1‑36 (MQ=255)
                                                                                                                        |                                                                                                                                     
TCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGA  >  NC_000913/1264850‑1265103

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