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% 65.2 / NA 22*399* (TAG→TAAhemYputative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/11);  total (11/11)

ACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGTTGTTTCAG  >  NC_000913/3986558‑3986820
                                                                                                                                |                                                                                                                                      
acaACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCg                                                                                                                              <  2:258098/139‑1 (MQ=255)
  aaCGTTGAGTGCAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCGTCTCCGTAATAAGAGCAGGCGTTTTTTTATGGAACAGGAAGGCCCCGGAGGCGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                                            <  1:592626/139‑1 (MQ=255)
                          ttgtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCaaa                                                                                                    <  2:471306/139‑1 (MQ=255)
                                tACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGGAGGCGTTGTTTTTAGCATCGGAAAGGTCCCGGATGTGCTTTCCGCCGCGTGGGGGGTAATTACGGTGGCGGGAGATGCTTCAACGCTAACaagataacgtct                                                                                              >  2:599795/1‑128 (MQ=255)
                                tACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                              >  2:131211/1‑139 (MQ=255)
                                                    tGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGcttc                                                                          >  2:620697/1‑139 (MQ=255)
                                                                 gagCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGtt                                                                                                             <  2:116166/91‑1 (MQ=255)
                                                                 gagCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGtt                                                                                                             >  1:116166/1‑91 (MQ=255)
                                                                  agCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCa                                                            <  2:640522/139‑1 (MQ=255)
                                                                   gCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                                   <  2:721918/131‑1 (MQ=255)
                                                                   gCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                                   >  1:721918/1‑131 (MQ=255)
                                                                                              gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAAc                                                                                                                >  2:592271/1‑59 (MQ=255)
                                                                                              gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                                >  2:215019/1‑139 (MQ=255)
                                                                                              gCCCCGGAGGTACTTGCCTCCGAGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAAc                                                                                                                <  1:592271/59‑1 (MQ=255)
                                                                                               ccccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCg                               >  2:775197/1‑139 (MQ=255)
                                                                                                cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcgc                                                                                             >  2:228547/1‑76 (MQ=255)
                                                                                                cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcgc                                                                                             <  1:228547/76‑1 (MQ=255)
                                                                                                cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGt                              >  2:357791/1‑139 (MQ=255)
                                                                                                cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGt                              >  1:465474/1‑139 (MQ=255)
                                                                                                                  cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                              >  2:285004/1‑57 (MQ=255)
                                                                                                                  cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                              <  1:285004/57‑1 (MQ=255)
                                                                                                                    ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACgt          <  2:465474/139‑1 (MQ=255)
                                                                                                                            ggAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGTTGTTTCAg  <  1:620697/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
ACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGTTGTTTCAG  >  NC_000913/3986558‑3986820

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