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

AAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACC  >  NC_000913/1264834‑1265105
                                                                                                                                        |                                                                                                                                       
aaGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACa                                                                                                                                       <  2:1518277/139‑1 (MQ=255)
                   ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCtttttt                                                                                                                    <  2:812132/139‑1 (MQ=255)
                   ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATCTTCtttttt                                                                                                                    >  1:1109945/1‑139 (MQ=255)
                          aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                  <  2:1015535/134‑1 (MQ=255)
                          aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                  >  1:1015535/1‑134 (MQ=255)
                          aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  2:981616/136‑1 (MQ=255)
                          aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  1:981616/1‑136 (MQ=255)
                          aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGCGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAggg                                                                                                             >  1:1362236/1‑139 (MQ=255)
                          aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGCGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAggg                                                                                                             <  2:1362236/139‑1 (MQ=255)
                             cTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGg                                                                                                          <  1:734170‑M1/139‑2 (MQ=255)
                              tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  2:1733652/1‑132 (MQ=255)
                              tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  1:1733652/132‑1 (MQ=255)
                              tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGga                                                                                                         <  2:1288133‑M1/139‑3 (MQ=255)
                                aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                       <  2:1695872‑M1/139‑5 (MQ=255)
                                aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                       <  2:1269484‑M1/139‑5 (MQ=255)
                                    accacttgaTCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggagga                                                                                                   >  2:949398‑M1/5‑131 (MQ=255)
                                        cTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACAGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataa                                                                                               >  2:1035733‑M1/1‑127 (MQ=255)
                                           gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTCCATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                            <  1:1056834‑M1/139‑16 (MQ=255)
                                           gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                            <  2:11555‑M1/139‑16 (MQ=255)
                                                aTGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                  >  2:553187/1‑112 (MQ=255)
                                                aTGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                                  <  1:553187/112‑1 (MQ=255)
                                                       aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacact                                                                                >  2:606097‑M1/1‑112 (MQ=255)
                                                                aCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  1:888294/98‑1 (MQ=255)
                                                                aCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  2:888294/1‑98 (MQ=255)
                                                                    cAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaac                                                                                   <  2:1520737‑M1/123‑25 (MQ=255)
                                                                    cAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaac                                                                                   >  1:1520737‑M1/1‑99 (MQ=255)
                                                                          gccgccCGTGACAGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataag                                                             <  1:1035733‑M1/139‑47 (MQ=255)
                                                                           ccgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagc                                                            >  1:1733097‑M1/1‑92 (MQ=255)
                                                                            cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                                           >  2:1348977‑M1/1‑91 (MQ=255)
                                                                               ccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGg                                                                                                          <  1:327662‑M1/89‑2 (MQ=255)
                                                                               ccGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGg                                                                                                          >  2:327662‑M1/1‑88 (MQ=255)
                                                                                    aCGGGGATAACGAACGCCTGAATATTCTGCGTGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  2:765388/78‑1 (MQ=255)
                                                                                    aCGGGGATAACGAACGCCTGAATATTCTGCGTGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  1:765388/1‑78 (MQ=255)
                                                                                              cgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaaca                                                                                  <  1:256896‑M1/98‑26 (MQ=255)
                                                                                              cgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaaca                                                                                  >  2:256896‑M1/1‑73 (MQ=255)
                                                                                                   gCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaata                                                                                        >  1:215941‑M1/1‑68 (MQ=255)
                                                                                                   gCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaata                                                                                        <  2:215941‑M1/87‑20 (MQ=255)
                                                                                                         atatTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgt                              >  2:1280566‑M1/1‑62 (MQ=255)
                                                                                                           atTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtga                            <  2:1109945‑M1/139‑80 (MQ=255)
                                                                                                               tGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggactatcaacactggttacgtgaagcaataagccaacttgaggcgagcgaaagcccgcggcg                               >  2:762823‑M1/1‑56 (MQ=255)
                                                                                                               tGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaactctggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                               <  1:762823‑M1/132‑77 (MQ=255)
                                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  1:26942/48‑1 (MQ=255)
                                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  2:26942/1‑48 (MQ=255)
                                                                                                                  gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaa                     >  2:1441453‑M1/1‑53 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  1:1034895/1‑44 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                >  1:1352958/1‑44 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  2:1352958/44‑1 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                                <  2:1034895/44‑1 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtg                             >  1:1313653‑M1/1‑49 (MQ=255)
                                                                                                                      cAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtg                             <  2:1313653‑M1/127‑79 (MQ=255)
                                                                                                                         ccTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgc              <  1:1280566‑M1/139‑94 (MQ=255)
                                                                                                                           tCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctg            >  1:1108064‑M1/1‑44 (MQ=255)
                                                                                                                             gggCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctgga          <  2:8251‑M1/139‑98 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgt    >  1:1378621‑M1/1‑36 (MQ=255)
                                                                                                                                   gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgt    >  2:924548‑M1/1‑36 (MQ=255)
                                                                                                                                    aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaaca                                                                                  <  2:137978‑M1/60‑26 (MQ=255)
                                                                                                                                    aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaaca                                                                                  >  1:137978‑M1/1‑35 (MQ=255)
                                                                                                                                     gTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgtta  <  2:1378621‑M1/139‑106 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACC  >  NC_000913/1264834‑1265105

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