Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,318,511	Δ1 bp	coding (254/759 nt)	phnK ←	carbon‑phosphorus lyase subunit PhnK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,318,508	0	C	.	100.0%	139.2 / NA	33	coding (257/759 nt)	phnK	carbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (17/16);  total (17/16)

GCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTT  >  NC_000913/4318430‑4318583                                                                               |                                                                            gCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGccc                                                                           <  1:818387/81‑1 (MQ=255) gCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGccc                                                                           <  2:818387/81‑1 (MQ=255) gCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGcc                                                                            <  2:818388/80‑1 (MQ=255)   cATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGcccat                                                                         >  1:818450/1‑79 (MQ=255)   cATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGcccat                                                                         <  2:818389/81‑3 (MQ=255)    aTCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGcccatt                                                                        >  1:818459/1‑78 (MQ=255)       aGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCgg                                                                    >  2:818457/1‑81 (MQ=255)        gCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGt                                                                   <  2:818390/81‑1 (MQ=255)          cGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTAc                                                                 <  1:818392/81‑1 (MQ=255)            cTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACgc                                                               <  1:818393/81‑1 (MQ=255)            cTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACg                                                                >  2:818461/1‑80 (MQ=255)            cTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACg                                                                >  1:818452/1‑80 (MQ=255)               gCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACgcagc                                                            <  2:818393/81‑1 (MQ=255)                ccGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACgcag                                                             >  1:818446/1‑79 (MQ=255)                  gATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAgg                                                         >  1:818442/1‑81 (MQ=255)                             cTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTc                                              <  1:818395/81‑1 (MQ=255)                               gCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCgg                                            >  1:818451/1‑81 (MQ=255)                                ccGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGc                                           >  1:818461/1‑81 (MQ=255)                                 cGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGcc                                          <  2:818395/81‑1 (MQ=255)                                   acacCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTc                                        <  1:818397/81‑1 (MQ=255)                                   acacCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTc                                        <  1:818398/81‑1 (MQ=255)                                    cacCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCg                                       <  1:818399/81‑1 (MQ=255)                                     acCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGc                                      >  1:818445/1‑81 (MQ=255)                                       cTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCt                                     >  1:818443/1‑80 (MQ=255)                                          gcggcgCAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCAt                                  <  1:818400/80‑1 (MQ=255)                                              cgcAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCa                              >  1:818465/1‑80 (MQ=255)                                               gcAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATa                            >  1:818466/1‑81 (MQ=255)                                               gcAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATa                            <  1:818401/81‑1 (MQ=255)                                                cAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATa                            >  1:818456/1‑80 (MQ=255)                                                cAGGCCGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATAc                           <  2:818401/81‑1 (MQ=255)                                                    ccGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGc                        >  2:818454/1‑80 (MQ=255)                                                    ccGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCg                       <  1:818403/81‑1 (MQ=255)                                                     cGTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGa                      <  1:818404/81‑1 (MQ=255)                                                      gTCGAGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGaa                     >  1:818471/1‑81 (MQ=255)                                                          aGTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGt                 >  1:818453/1‑81 (MQ=255)                                                           gTGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGtt                >  2:818453/1‑81 (MQ=255)                                                            tGGATGCTGATGCACCACG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTtc               >  2:818474/1‑81 (MQ=255)                                                                        caccacG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAAtt   <  2:818404/81‑1 (MQ=255)                                                                         accacG‑CCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAAttt  <  1:818406/81‑1 (MQ=255)                                                                               |                                                                            GCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTT  >  NC_000913/4318430‑4318583

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