Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,291,973 T→C G146G (GGA→GGG ccmH ← heme lyase, CcmH subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,9730TC100.0% 97.4 / NA 28G146G (GGA→GGGccmHheme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/15);  total (13/15)

CGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCAC  >  NC_000913/2291838‑2292107
                                                                                                                                       |                                                                                                                                      
cgcCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCg                                                                                                                                     <  2:949298/139‑1 (MQ=255)
 gcCCTGTCCAGCAACGCCGGAGCGTGTGCCGTGGCCTGCTTCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                                    <  2:1035176/139‑1 (MQ=255)
 gcCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                                    <  1:901238/139‑1 (MQ=255)
                 ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAGTTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    >  2:189893/1‑139 (MQ=255)
                 ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    >  2:909233/1‑139 (MQ=255)
                           tgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                                                  <  1:391611/115‑1 (MQ=255)
                           tgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                                                  >  2:391611/1‑115 (MQ=255)
                                  gCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACg                                                                                                                <  2:123270/126‑1 (MQ=255)
                                  gCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACg                                                                                                                >  1:123270/1‑126 (MQ=255)
                                                            aTAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGa                                                                         >  1:543615/1‑139 (MQ=255)
                                                                     aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                <  2:543615/139‑1 (MQ=255)
                                                                          ggTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGc                                                           <  2:127940/139‑1 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                                    >  2:848260/1‑139 (MQ=255)
                                                                                      cGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCa                                               >  1:183183/1‑139 (MQ=255)
                                                                                            ccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcc                                         >  2:504923/1‑139 (MQ=255)
                                                                                              acaATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaa                                       <  1:848260/139‑1 (MQ=255)
                                                                                                   tAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca                                  <  1:189893/139‑1 (MQ=255)
                                                                                                      cGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                   >  1:758342/1‑55 (MQ=255)
                                                                                                      cGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                   <  2:758342/55‑1 (MQ=255)
                                                                                                           ccacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCac                          >  1:753383/1‑139 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    <  2:705670/41‑1 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    >  1:705670/1‑41 (MQ=255)
                                                                                                                         aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGc                             <  2:25359/122‑1 (MQ=255)
                                                                                                                         aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGc                             >  1:25359/1‑122 (MQ=255)
                                                                                                                           aTAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGGAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcca                                        <  1:89680/109‑1 (MQ=255)
                                                                                                                           aTAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGGAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcca                                        >  2:89680/1‑109 (MQ=255)
                                                                                                                                 aacaTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACagcagc    <  1:504923/139‑1 (MQ=255)
                                                                                                                                   caTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCAc  <  1:1049462/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCAC  >  NC_000913/2291838‑2292107

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