Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,318,564	A→G	R67R (CGT→CGC)	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,564	0	A	G	100.0%	189.3 / NA	54	R67R (CGT→CGC)	phnK	carbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/29);  total (25/29)

CTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGC  >  NC_000913/4318441‑4318700                                                                                                                            |                                                                                                                                         cTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGaa                                                                                                                          >  2:58978/1‑140 (MQ=255)  tCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAAt                                                                                                                         >  2:367996/1‑140 (MQ=255)  tCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAAt                                                                                                                         >  1:923234/1‑140 (MQ=255)      cGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTcc                                                                                                                     >  1:1035559/1‑140 (MQ=255)                        acacCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAgg                                                                                                   <  2:578270/140‑1 (MQ=255)                            cTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGgcgcg                                                                                               >  2:1058801/1‑140 (MQ=255)                                  gcgcAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGgaga                                                                                         <  2:923234/140‑1 (MQ=255)                                    gcAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGTGCGGAGATc                                                                                       >  1:311755/1‑140 (MQ=255)                                    gcAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATc                                                                                       >  2:670015/1‑140 (MQ=255)                                       ggCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGAc                                                                                    <  2:926865/140‑1 (MQ=255)                                           gTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTca                                                                                <  1:536531/140‑1 (MQ=255)                                                     tGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGt                                                                      <  2:428044/140‑1 (MQ=255)                                                          aTGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                 <  1:367996/140‑1 (MQ=255)                                                          aTGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                 >  2:1391911/1‑140 (MQ=255)                                                           tGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                >  2:1113464/1‑140 (MQ=255)                                                             caccacGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                              >  1:73037/1‑140 (MQ=255)                                                                  cGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCg                                                         >  1:1347663/1‑140 (MQ=255)                                                                  cGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCg                                                         <  1:257974/140‑1 (MQ=255)                                                                  cGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCg                                                         >  1:716910/1‑140 (MQ=255)                                                                    ccccATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                       <  1:895875/140‑1 (MQ=255)                                                                    ccccATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                       <  2:605563/140‑1 (MQ=255)                                                                     cccATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAt                                                      >  1:1363249/1‑140 (MQ=255)                                                                      ccATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAtt                                                     <  1:1344505/140‑1 (MQ=255)                                                                       cATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTc                                                    <  1:670015/140‑1 (MQ=255)                                                                       cATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTc                                                    <  1:1391911/140‑1 (MQ=255)                                                                         ttCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTccc                                                  >  1:166580/1‑140 (MQ=255)                                                                             gTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAc                                              >  1:1422257/1‑140 (MQ=255)                                                                              tACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACa                                             <  1:1313923/140‑1 (MQ=255)                                                                                 gcagcaGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATg                                          >  1:1230489/1‑140 (MQ=255)                                                                                  cagcagGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTCCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGc                                         <  1:1369479/140‑1 (MQ=255)                                                                                   agcagGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGcc                                        <  2:1347663/140‑1 (MQ=255)                                                                                   agcagGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGcc                                        >  2:774415/1‑140 (MQ=255)                                                                                    gcagGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                       <  1:58978/140‑1 (MQ=255)                                                                                       ggCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGc                                    >  1:493873/1‑140 (MQ=255)                                                                                       ggCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGc                                    >  2:863508/1‑140 (MQ=255)                                                                                        gCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                   <  1:863508/140‑1 (MQ=255)                                                                                        gCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                   >  1:1126565/1‑140 (MQ=255)                                                                                          gACGGCGGTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                 <  2:778691/140‑1 (MQ=255)                                                                                               gggTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcccc                            <  2:423719/139‑1 (MQ=255)                                                                                                ggTCGGCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccccc                           <  2:1035559/140‑1 (MQ=255)                                                                                                     gCCTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                      <  2:166580/140‑1 (MQ=255)                                                                                                       cTCGCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                    <  2:493873/140‑1 (MQ=255)                                                                                                          gCTCATTGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAAt                 <  2:716910/140‑1 (MQ=255)                                                                                                                tGCATACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag           >  2:455621/1‑140 (MQ=255)                                                                                                                    tACAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGAc       <  1:1058801/140‑1 (MQ=255)                                                                                                                     aCAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa      <  2:73037/140‑1 (MQ=255)                                                                                                                     aCAGCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa      <  2:1363249/140‑1 (MQ=255)                                                                                                                        gCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa      >  2:743715/1‑137 (MQ=255)                                                                                                                        gCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa      <  1:743715/137‑1 (MQ=255)                                                                                                                        gCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg   >  2:83083/1‑140 (MQ=255)                                                                                                                        gCGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATCCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg   >  2:431103/1‑140 (MQ=255)                                                                                                                         cGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                    <  1:155615/122‑1 (MQ=255)                                                                                                                         cGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa      <  1:1045690/136‑1 (MQ=255)                                                                                                                         cGAGCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc  >  2:820486/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                         CTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCCCATTCGGTACGCAGCAGGCGACGGCGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGC  >  NC_000913/4318441‑4318700

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