Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,293,106 (G)7→8 coding (1595/1659 nt) nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,293,0991.G100.0% 109.9 / NA 29W530G (TGG→GGG) nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/14);  total (15/14)

CCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCG  >  NC_000913/4292969‑4293223
                                                                                                                                   |                                                                                                                            
ccGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTgggggggg                                                                                                                      <  2:1098670/139‑2 (MQ=255)
  gTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTggggggggag                                                                                                                    >  2:458101/1‑136 (MQ=255)
    caacaaATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTggggggggagga                                                                                                                  <  1:579567/139‑6 (MQ=255)
          atgatgGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGAttgttg                                                                                                             <  1:375584/139‑1 (MQ=255)
                      tCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGgc                                                                                                 >  1:188296/1‑139 (MQ=255)
                      tCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGgc                                                                                                 >  2:395480/1‑139 (MQ=255)
                         aTTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCa                                                                                              <  2:357058/139‑1 (MQ=255)
                                                                         aCTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTtct                                              <  1:395480/139‑1 (MQ=255)
                                                                             ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACg                                          <  2:567174/139‑1 (MQ=255)
                                                                              gCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGc                                         <  1:458101/139‑1 (MQ=255)
                                                                                               ttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTAcc                        >  2:747244/1‑139 (MQ=255)
                                                                                               ttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTgg                                                                                                                            >  1:327792/1‑38 (MQ=255)
                                                                                               ttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTgg                                                                                                                            <  2:327792/38‑1 (MQ=255)
                                                                                                     tgtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       >  2:963281/1‑117 (MQ=255)
                                                                                                     tgtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       <  1:963281/117‑1 (MQ=255)
                                                                                                      gtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTActgc                                   >  1:1200735/1‑121 (MQ=255)
                                                                                                      gtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTActgc                                   <  2:1200735/121‑1 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTtgat                                                                                                           <  1:1001528/39‑1 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTtgat                                                                                                           >  2:1001528/1‑39 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCa                                                                                              >  2:357537/1‑52 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCa                                                                                              <  1:357537/52‑1 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       <  1:391049/107‑1 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       >  2:391049/1‑107 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTAtt                            <  2:626122/118‑1 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTAtt                            >  1:626122/1‑118 (MQ=255)
                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCa       >  2:198650/1‑139 (MQ=255)
                                                                                                                    gtgtGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGtt                    >  1:94090/1‑122 (MQ=255)
                                                                                                                    gtgtGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGtt                    <  2:94090/122‑1 (MQ=255)
                                                                                                                     tgtgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCg  >  1:1074091/1‑139 (MQ=255)
                                                                                                                     tgtgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCg  >  1:250391/1‑139 (MQ=255)
                                                                                                                       tgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCt                                                 <  2:145297/90‑1 (MQ=255)
                                                                                                                       tgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCt                                                 >  1:145297/1‑90 (MQ=255)
                                                                                                                          gcTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGAt                                                                                 <  1:989312/55‑1 (MQ=255)
                                                                                                                          gcTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGAt                                                                                 >  2:989312/1‑55 (MQ=255)
                                                                                                                                   |                                                                                                                            
CCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCG  >  NC_000913/4292969‑4293223

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