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

CCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCAC  >  NC_000913/2291855‑2292107
                                                                                                                      |                                                                                                                                      
ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    >  1:329191/1‑139 (MQ=255)
          tgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGt                                                                                                          <  1:345342/139‑1 (MQ=255)
                    tgctgcCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATAccc                                                                                                                                      >  2:640127/1‑101 (MQ=255)
                                                    aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                <  2:447437/139‑1 (MQ=255)
                                                      tCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCg                                                              <  1:703138/139‑1 (MQ=255)
                                                       cTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    <  2:651434/84‑1 (MQ=255)
                                                       cTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                                    >  1:651434/1‑84 (MQ=255)
                                                          gtagtaGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCg                                                          >  1:292255/1‑139 (MQ=255)
                                                                 cTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGa                                                                         <  1:62558/117‑1 (MQ=255)
                                                                 cTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGa                                                                         >  2:62558/1‑117 (MQ=255)
                                                                    aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                              <  2:155842/77‑1 (MQ=255)
                                                                    aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                              >  1:155842/1‑77 (MQ=255)
                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                   <  1:355195/39‑1 (MQ=255)
                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                   >  2:355195/1‑39 (MQ=255)
                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                                      <  2:27285/68‑1 (MQ=255)
                                                                                                     ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                                      >  1:27285/1‑68 (MQ=255)
                                                                                                                  caTACCCGGCACGCTGACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCAc  <  2:208584/139‑1 (MQ=255)
                                                                                                                      |                                                                                                                                      
CCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCAC  >  NC_000913/2291855‑2292107

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