Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,321,059	C→T	W63* (TGG→TGA)	phnH ←	carbon‑phosphorus lyase core complex subunit PhnH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,321,059	0	C	T	100.0%	138.8 / NA	40	W63* (TGG→TGA)	phnH	carbon‑phosphorus lyase core complex subunit PhnH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

GGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGA  >  NC_000913/4320924‑4321191                                                                                                                                        |                                                                                                                                     ggCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTcaca                                                                                                                                  <  1:1004128/140‑1 (MQ=255)   cGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGCATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACAcc                                                                                                                                >  2:11235/1‑140 (MQ=255)      tGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGc                                                                                                                             >  1:547946/1‑140 (MQ=255)             ttCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGt                                                                                                                      >  1:896757/1‑140 (MQ=255)             ttCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGt                                                                                                                      >  1:7334/1‑140 (MQ=255)                   ggAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCgg                                                                                                                <  2:681796/140‑1 (MQ=255)                     aaaTCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGcc                                                                                                              <  2:1136145/140‑1 (MQ=255)                            cTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCAGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTca                                                                                                       >  2:717706/1‑140 (MQ=255)                             tCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTcag                                                                                                      <  1:337456/140‑1 (MQ=255)                                      aCCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCt                                                                                             <  2:728750/140‑1 (MQ=255)                                           gAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtgg                                                                                        <  1:717706/140‑1 (MQ=255)                                                     tGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATAtt                                                                              <  2:193126/140‑1 (MQ=255)                                                                 ctgACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCa                                                                  >  2:130773/1‑140 (MQ=255)                                                                     ccAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCcg                                                              <  1:130773/140‑1 (MQ=255)                                                                         cggcgCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAACAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGtt                                                          <  2:582623/140‑1 (MQ=255)                                                                                    tATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGc                                               <  1:11235/140‑1 (MQ=255)                                                                                      tGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAg                                             <  2:55060/140‑1 (MQ=255)                                                                                              cAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGc                                               <  1:847281/130‑1 (MQ=255)                                                                                              cAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGc                                               >  2:847281/1‑130 (MQ=255)                                                                                              cAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACaa                                     >  2:1081505/1‑140 (MQ=255)                                                                                               aGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAg                                             >  2:1132649/1‑131 (MQ=255)                                                                                               aGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAg                                             <  1:1132649/131‑1 (MQ=255)                                                                                                 gCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCa                                  >  1:1015744/1‑140 (MQ=255)                                                                                                       gTTGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCg                            >  2:1160156/1‑140 (MQ=255)                                                                                                         tGACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgc                          <  2:1015744/140‑1 (MQ=255)                                                                                                          gACGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgct                         <  1:1081505/140‑1 (MQ=255)                                                                                                           aCGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctc                        >  1:724802/1‑140 (MQ=255)                                                                                                           aCGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctc                        >  2:468294/1‑140 (MQ=255)                                                                                                           aCGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgct                         >  2:1159748/1‑139 (MQ=255)                                                                                                           aCGATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgct                         <  1:1159748/139‑1 (MQ=255)                                                                                                             gATATCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcgc                      >  2:837745/1‑140 (MQ=255)                                                                                                               tatCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcgctc                    >  2:403329/1‑140 (MQ=255)                                                                                                                 tCGTTATTTAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCAt                  >  1:1040784/1‑140 (MQ=255)                                                                                                                         tAATGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTa          >  1:788771/1‑140 (MQ=255)                                                                                                                          aaTGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTaa         <  1:468294/140‑1 (MQ=255)                                                                                                                           aTGGGGTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAAc        >  1:600137/1‑140 (MQ=255)                                                                                                                             ggggTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAg      <  2:788771/140‑1 (MQ=255)                                                                                                                             ggggTAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAg      <  1:815240/140‑1 (MQ=255)                                                                                                                                 tAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGcga  >  2:1161883/1‑140 (MQ=255)                                                                                                                                 tAGAAAGTCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGcga  >  2:1105264/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                     GGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGA  >  NC_000913/4320924‑4321191

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: ‑