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

CCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCG  >  NC_000913/1264837‑1265089
                                                                                                                                     |                                                                                                                       
ccATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTa                                                                                                                    <  2:1012852/139‑1 (MQ=255)
           gATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCt                                                                                                         >  1:1121021/1‑139 (MQ=255)
                  cATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTCACAGTACATCATTTTCTTTTTTTa                                                                                                  <  2:1121021/139‑1 (MQ=255)
                    tGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                <  1:1460572/139‑1 (MQ=255)
                    tGGAAACGGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACTAACGCCTGCATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                <  1:1042387/139‑1 (MQ=255)
                             aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       <  2:1067743‑M1/139‑5 (MQ=255)
                             aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       <  1:766608‑M1/139‑5 (MQ=255)
                             aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       <  2:121321‑M1/139‑5 (MQ=255)
                                 aCCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggagga                                                                                   >  1:975785‑M1/1‑131 (MQ=255)
                                  ccGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaa                                                                                  >  1:899856‑M1/1‑130 (MQ=255)
                                      ttGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataat                                                                              <  1:713544‑M1/139‑14 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                               <  2:428891/121‑1 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                               >  1:428891/1‑121 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaata                                                                                <  2:487306‑M1/136‑12 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaata                                                                                >  1:487306‑M1/1‑125 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       >  2:1746432‑M1/1‑125 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       <  1:1746432‑M1/129‑5 (MQ=255)
                                       tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTCgaggaggaataatg                                                                             >  2:1129965‑M1/1‑125 (MQ=255)
                                        gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                            <  1:1129965‑M1/139‑16 (MQ=255)
                                        gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                            <  2:848934‑M1/139‑16 (MQ=255)
                                        gATCCACGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                            <  1:1707206‑M1/139‑16 (MQ=255)
                                                                 cAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagga                                                                                      <  1:1386682‑M1/104‑6 (MQ=255)
                                                                 cAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagga                                                                                      >  2:1386682‑M1/1‑99 (MQ=255)
                                                                       gccgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataag                                             >  1:1385201‑M1/1‑93 (MQ=255)
                                                                         cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                           >  1:1352734‑M1/1‑91 (MQ=255)
                                                                         cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                           >  1:794716‑M1/1‑91 (MQ=255)
                                                                         cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       <  2:1749068‑M1/95‑5 (MQ=255)
                                                                         cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       >  1:231137‑M1/1‑91 (MQ=255)
                                                                         cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       <  2:231137‑M1/95‑5 (MQ=255)
                                                                         cgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                       >  1:1749068‑M1/1‑91 (MQ=255)
                                                                                 atggggatAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACACTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcag                                   <  2:975785‑M1/137‑57 (MQ=255)
                                                                                         aacgaacgCCTGGATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaat                                                <  1:73459‑M1/118‑44 (MQ=255)
                                                                                         aacgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggcggaataatggaatatcaacactggttacgtgaagcaat                                                >  2:73459‑M1/1‑75 (MQ=255)
                                                                                          acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                <  1:1052559/69‑1 (MQ=255)
                                                                                          acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                >  2:1052559/1‑69 (MQ=255)
                                                                                          acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaa                                                                    <  2:358464‑M1/97‑24 (MQ=255)
                                                                                          acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaa                                                                    >  1:358464‑M1/1‑74 (MQ=255)
                                                                                          acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                            <  1:1665168‑M1/89‑16 (MQ=255)
                                                                                          acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                            >  2:1665168‑M1/1‑74 (MQ=255)
                                                                                              acgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                  >  1:725519/1‑63 (MQ=255)
                                                                                              acgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                  <  2:725519/63‑1 (MQ=255)
                                                                                                 ccTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                >  2:1762522/1‑62 (MQ=255)
                                                                                                 ccTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                <  1:1762522/62‑1 (MQ=255)
                                                                                                     aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                           <  1:793708‑M1/111‑49 (MQ=255)
                                                                                                     aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagcc                                           >  2:793708‑M1/1‑63 (MQ=255)
                                                                                                            tGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggagg                                                                                    >  2:1725941‑M1/1‑56 (MQ=255)
                                                                                                            tGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggagg                                                                                    <  1:1725941‑M1/63‑8 (MQ=255)
                                                                                                              cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                  <  1:724782/47‑1 (MQ=255)
                                                                                                              cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTa                                                                                                  >  2:724782/1‑47 (MQ=255)
                                                                                                              cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGt                                                                                            <  1:110259/53‑1 (MQ=255)
                                                                                                              cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGt                                                                                            >  2:110259/1‑53 (MQ=255)
                                                                                                               gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaa     >  1:499223‑M1/1‑53 (MQ=255)
                                                                                                                  aCAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatc  <  2:1385201‑M1/139‑90 (MQ=255)
                                                                                                                                     |                                                                                                                       
CCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCG  >  NC_000913/1264837‑1265089

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