Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,986,686 C→T *399* (TAG→TAA hemY ← putative protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,986,6860CT100.0% 80.7 / NA 26*399* (TAG→TAAhemYputative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/11);  total (15/11)

TGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGC  >  NC_000913/3986553‑3986802
                                                                                                                                     |                                                                                                                    
tGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtg                                                                                                                 >  1:569003/1‑139 (MQ=255)
    cacaACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGc                                                                                                             <  1:757596/139‑1 (MQ=255)
       aaCGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                          <  2:569003/139‑1 (MQ=255)
                tGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                 <  2:15518/139‑1 (MQ=255)
                                  tCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCg                                                                               >  1:732182/1‑139 (MQ=255)
                                     ttcAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCTACGTTAACATCAAACCGTcg                                                                            <  2:452987/137‑1 (MQ=255)
                                          aCCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                       >  1:246297/1‑139 (MQ=255)
                                                         tGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGcttc                                                        >  1:263381/1‑139 (MQ=255)
                                                                          aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTc                                       >  1:597137/1‑139 (MQ=255)
                                                                                            cAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                 <  2:644290/111‑1 (MQ=255)
                                                                                            cAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                 >  1:644290/1‑111 (MQ=255)
                                                                                               gAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGcc                  <  1:703769/139‑1 (MQ=255)
                                                                                                   gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                          <  1:392261/47‑1 (MQ=255)
                                                                                                   gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                          >  2:392261/1‑47 (MQ=255)
                                                                                                   gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCg                                                                                                            <  2:126761/45‑1 (MQ=255)
                                                                                                   gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCg                                                                                                            >  1:126761/1‑45 (MQ=255)
                                                                                                   gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc              >  1:258920/1‑139 (MQ=255)
                                                                                                   gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc              >  2:573900/1‑139 (MQ=255)
                                                                                                       cGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTaa          >  1:262811/1‑139 (MQ=255)
                                                                                                             tGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCgg                                                  >  1:100016/1‑93 (MQ=255)
                                                                                                             tGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCgg                                                  <  2:100016/93‑1 (MQ=255)
                                                                                                               cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCAt                                                                    <  1:385248/73‑1 (MQ=255)
                                                                                                               cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCAt                                                                    >  2:385248/1‑73 (MQ=255)
                                                                                                               cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGc  >  1:362966/1‑139 (MQ=255)
                                                                                                                         ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCGGTCTGTCCAGCGCGTCGGCAAGCCATGCGt            <  1:429351/119‑1 (MQ=255)
                                                                                                                         ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCGGTCTGTCCAGCGCGTCGGCAAGCCATGCGt            >  2:429351/1‑119 (MQ=255)
                                                                                                                                     |                                                                                                                    
TGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGC  >  NC_000913/3986553‑3986802

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