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

GATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGAC  >  NC_000913/1264848‑1265104
                                                                                                                          |                                                                                                                                      
gATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCt                                                                                                                        >  1:1164706/1‑139 (MQ=255)
     ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCtttttt                                                                                                                   >  1:473214/1‑139 (MQ=255)
            aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAggg                                                                                                            >  1:994081/1‑139 (MQ=255)
                    tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggag                                                                                                    <  1:310348‑M1/139‑7 (MQ=255)
                    tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggag                                                                                                    <  2:473214‑M1/139‑7 (MQ=255)
                                   tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                              >  1:402307/1‑114 (MQ=255)
                                   tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                              >  1:402308/1‑114 (MQ=255)
                                   tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                              <  2:402308/114‑1 (MQ=255)
                                   tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                              <  2:402307/114‑1 (MQ=255)
                                   tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatc                                                                                     >  2:554215‑M1/1‑118 (MQ=255)
                                   tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatc                                                                                     >  1:929954‑M1/1‑118 (MQ=255)
                                         aaCGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                           >  1:1266274‑M1/1‑112 (MQ=255)
                                         aaCGAAATCACTTCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                           <  2:1266274‑M1/127‑16 (MQ=255)
                                                     tCAACAGGCCGCCCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaag                                                                   >  1:703317‑M1/1‑100 (MQ=255)
                                                             ccgccCGTGATGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagc                                                           <  2:929954‑M1/139‑48 (MQ=255)
                                                                                    cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagccc                                    >  1:633989‑M1/1‑69 (MQ=255)
                                                                                      ccTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgc                                  >  2:558463‑M1/1‑67 (MQ=255)
                                                                                          aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACTTCATTTTCTTTTTTTACAGGGTGgagaaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagtgaaagcccgcggcg                              <  2:831124‑M1/139‑77 (MQ=255)
                                                                                          aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcg                              <  1:554215‑M1/139‑77 (MQ=255)
                                                                                          aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                           >  1:479124‑M1/1‑63 (MQ=255)
                                                                                          aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatgg                                                                                           <  2:479124‑M1/78‑16 (MQ=255)
                                                                                             atTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtga                           <  1:558463‑M1/139‑80 (MQ=255)
                                                                                                    gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatc                                                                                     >  2:1273267‑M1/1‑53 (MQ=255)
                                                                                                    gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatc                                                                                     <  1:1273267‑M1/74‑22 (MQ=255)
                                                                                                                      aGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgtt  >  2:822126‑M1/1‑35 (MQ=255)
                                                                                                                          |                                                                                                                                      
GATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGAC  >  NC_000913/1264848‑1265104

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