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% 110.8 / NA 32G146G (GGA→GGGccmHheme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/15);  total (17/15)

GCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACC  >  NC_000913/2291839‑2292096
                                                                                                                                      |                                                                                                                           
gcCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                         <  1:264307/139‑1 (MQ=255)
                    aGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACg                                                                                                     <  2:805453/139‑1 (MQ=255)
                                                 tCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATAccc                                                                                                                           <  1:1255798/88‑1 (MQ=255)
                                                       gctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCg                                                                                                                          <  2:755893/83‑1 (MQ=255)
                                                       gctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCg                                                                                                                          >  1:755893/1‑83 (MQ=255)
                                                        ctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAcc                                                                 >  1:1695/1‑139 (MQ=255)
                                                                        tGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                 <  2:1435067/139‑1 (MQ=255)
                                                                          gtagtaGCTGACGCCAGCCACAATTAGCGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                        <  2:305371/82‑1 (MQ=255)
                                                                          gtagtaGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                        >  1:305371/1‑82 (MQ=255)
                                                                          gtagtaGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCg                                               >  1:1253175/1‑139 (MQ=255)
                                                                                gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                     >  2:560791/1‑127 (MQ=255)
                                                                                gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                     <  1:560791/127‑1 (MQ=255)
                                                                                gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                         >  1:428781/1‑139 (MQ=255)
                                                                                                    gcGCCAACACAGTACCCGGCAGATAACCAACATACCCGGCCGCCTCACCTTCGGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGCGCGGGCGTATATGACCCAACCgccaatgccaat                     <  1:713525/138‑1 (MQ=255)
                                                                                                         accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaaa                                                                                       >  1:1275704/1‑68 (MQ=255)
                                                                                                         accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaaa                                                                                       <  2:1275704/68‑1 (MQ=255)
                                                                                                         accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                           <  2:1729903/80‑1 (MQ=255)
                                                                                                         accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                           >  1:1729903/1‑80 (MQ=255)
                                                                                                         accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCa                >  2:872656/1‑139 (MQ=255)
                                                                                                                tACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGc                                                    >  1:1245686/1‑96 (MQ=255)
                                                                                                                tACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGc                                                    <  2:1245686/96‑1 (MQ=255)
                                                                                                                tACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTgg         >  2:284525/1‑139 (MQ=255)
                                                                                                                 aCCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCtt                                                                                                  <  1:537349/49‑1 (MQ=255)
                                                                                                                 aCCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCtt                                                                                                  >  2:537349/1‑49 (MQ=255)
                                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                   >  2:1608164/1‑44 (MQ=255)
                                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                   <  1:1608164/44‑1 (MQ=255)
                                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                     >  1:1784895/1‑90 (MQ=255)
                                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                     <  2:1784895/90‑1 (MQ=255)
                                                                                                                       cAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCAcc  >  2:556164/1‑139 (MQ=255)
                                                                                                                            aaaCAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcc                              <  2:1655794/106‑1 (MQ=255)
                                                                                                                            aaaCAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcc                              >  1:1655794/1‑106 (MQ=255)
                                                                                                                                  caTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                     <  2:1824493/77‑1 (MQ=255)
                                                                                                                                  caTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                     >  1:1824493/1‑77 (MQ=255)
                                                                                                                                      |                                                                                                                           
GCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACC  >  NC_000913/2291839‑2292096

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