Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,295,852	G→A	R157C (CGT→TGT)	ccmC ←	heme trafficking system membrane protein CcmC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,295,852	0	G	A	100.0%	163.8 / NA	48	R157C (CGT→TGT)	ccmC	heme trafficking system membrane protein CcmC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/22);  total (26/22)

CGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAG  >  NC_000913/2295718‑2295988                                                                                                                                       |                                                                                                                                         cggcgAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCg                                                                                                                                     <  2:155454/140‑1 (MQ=255)      aaCGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGTCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGa                                                                                                                                >  2:128946/1‑140 (MQ=255)      aaCGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGa                                                                                                                                >  2:541578/1‑140 (MQ=255)         gCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGAcgg                                                                                                                             <  1:689251/140‑1 (MQ=255)                        tACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCg                                                                                                              <  2:378343/140‑1 (MQ=255)                        tACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCg                                                                                                              <  2:249661/140‑1 (MQ=255)                             tgctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCca                                                                                                         >  2:129119/1‑140 (MQ=255)                                  cATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCTGCCAGACGGCGGTCGTCGAAGGCGTGCCACAggg                                                                                                    >  2:657216/1‑140 (MQ=255)                                  cATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAggg                                                                                                    >  2:577302/1‑140 (MQ=255)                                   aTCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGc                                                                                                   <  1:132154/140‑1 (MQ=255)                                     ccGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCAGCCAGACGGCGGTCGTCGAAGACGTGCCACAGGGCaa                                                                                                 >  1:225725/1‑140 (MQ=255)                                            gATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACAcc                                                                                          <  2:567291/140‑1 (MQ=255)                                             aTCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACAccc                                                                                         >  2:684483/1‑140 (MQ=255)                                                    aTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACa                                                                                  >  2:156061/1‑140 (MQ=255)                                                          ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAAGCGTGCCTCAGGGCTATCACATCCACATACAAAAAca                                                                            >  1:208404/1‑140 (MQ=255)                                                                 caccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAc                                                                     <  2:405852/140‑1 (MQ=255)                                                                    cacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAg                                                                  >  2:518777/1‑140 (MQ=255)                                                                      cTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCATACCTACGTACAAATACAGCAGCACCAGtt                                                                >  1:345380/1‑140 (MQ=255)                                                                      cTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt                                                                >  2:537786/1‑140 (MQ=255)                                                                      cTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt                                                                >  1:67382/1‑140 (MQ=255)                                                                      cTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt                                                                >  2:645055/1‑140 (MQ=255)                                                                         cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAg                                                             >  2:674558/1‑140 (MQ=255)                                                                               gTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCa                                                       <  1:156061/140‑1 (MQ=255)                                                                                 aaTGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACGCCCCCATACAAAAACATCAGCGCCAGTTCAGAAGTCAGa                                                     >  2:222929/1‑140 (MQ=255)                                                                                     aataatCGGCAGATTCACCACGCCAATCAGCACCAGGTTACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTg                                                 <  1:290115/140‑1 (MQ=255)                                                                                     aataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTg                                                 <  2:518884/140‑1 (MQ=255)                                                                                      ataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGc                                                <  1:518777/140‑1 (MQ=255)                                                                                      ataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGc                                                <  1:577302/140‑1 (MQ=255)                                                                                          tCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACACAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                            >  1:259808/1‑140 (MQ=255)                                                                                          tCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                            >  2:469159/1‑140 (MQ=255)                                                                                          tCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                            >  2:15028/1‑140 (MQ=255)                                                                                            ggCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCa                                          <  1:129119/140‑1 (MQ=255)                                                                                            ggCAGATTCACCACGCCAATCAGCAACAGAATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCa                                          <  1:317688/140‑1 (MQ=255)                                                                                               aGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAc                                       >  1:562340/1‑140 (MQ=255)                                                                                                          gCCAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGc                            <  2:562340/140‑1 (MQ=255)                                                                                                           ccAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAGGTCAGACGTGCATCCCATACCCACCAGGTGcc                           >  2:462992/1‑140 (MQ=255)                                                                                                           ccAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGcc                           >  2:386589/1‑140 (MQ=255)                                                                                                            cAATCAGCACCAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGccc                          >  1:35880/1‑140 (MQ=255)                                                                                                                     ccAGGATACCTGCCGCACAGCCCGCTAGATGGCGGTCCTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGt                 <  1:222929/140‑1 (MQ=255)                                                                                                                     ccAGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGt                 <  1:469159/140‑1 (MQ=255)                                                                                                                       gggATACCTGCCGCATAGCCCGGCATACGTCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGttt               <  1:242220/139‑1 (MQ=255)                                                                                                                       aGGATACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGttt               <  1:541578/140‑1 (MQ=255)                                                                                                                          aTACCTGCCGCACAGCCTGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTcc            <  1:657216/140‑1 (MQ=255)                                                                                                                          aTACCTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTcc            <  2:35880/140‑1 (MQ=255)                                                                                                                             ccTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCa         >  2:26153/1‑140 (MQ=255)                                                                                                                             ccTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCa         >  2:173740/1‑140 (MQ=255)                                                                                                                              cTGCCGCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCAt        >  2:467406/1‑140 (MQ=255)                                                                                                                                  cGCACAGCCCGCCAGACCGCGGTCGTCGAAGTCGTTCCACAGGGCAATTACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag    <  2:225725/140‑1 (MQ=255)                                                                                                                                   gCACAGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACATGCATCCCATACGCACGAGGTGCCCGACATGGGTTTTCCCCATGCagc   >  2:242220/1‑139 (MQ=255)                                                                                                                                    cacaGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCAAACCGACCAGGTGCCCCGCATCCGTTTTCCCCATGCagag  >  2:221151/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         CGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAG  >  NC_000913/2295718‑2295988

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