Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,318,211	A→G	L185P (CTC→CCC)	phnK ←	carbon‑phosphorus lyase complex subunit, putative ATP transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,318,211	0	A	G	100.0%	200.8 / NA	53	L185P (CTC→CCC)	phnK	carbon‑phosphorus lyase complex subunit, putative ATP transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/23);  total (30/23)

ACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGT  >  NC_000913/4318080‑4318343                                                                                                                                    |                                                                                                                                     aCCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCg                                                                                                                               >  2:1717554/1‑139 (MQ=255)                 ttGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGAcaca                                                                                                               <  2:1294889/139‑1 (MQ=255)                   gCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATc                                                                                                             >  1:626452/1‑139 (MQ=255)                          cTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCccg                                                                                                      >  1:1359434/1‑139 (MQ=255)                                       aaaCGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTcgt                                                                                         <  1:1717554/139‑3 (MQ=255)                                             tCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCa                                                                                      >  2:1399763/1‑136 (MQ=255)                                             tCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCa                                                                                      <  1:1399763/136‑1 (MQ=255)                                              ccGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATaaa                                                                                  >  1:1476388/1‑139 (MQ=255)                                                gCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAaca                                                                                <  1:1593358/139‑1 (MQ=255)                                                gCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCT‑‑‑CCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCAt                                                                                         >  2:852739/1‑127 (MQ=255)                                                gCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCT‑‑‑CCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCAt                                                                                         <  1:852739/127‑1 (MQ=255)                                                       ggcgggcgACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAgg                                                                                                                                 <  1:1074428/83‑1 (MQ=255)                                                       ggcgggcgACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAgg                                                                                                                                 >  2:1074428/1‑83 (MQ=255)                                                       ggcgggcgACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGtt                                                                                            >  2:26539/1‑120 (MQ=255)                                                       ggcgggcgACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGtt                                                                                            <  1:26539/120‑1 (MQ=255)                                                       ggcgggcgACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCtt                                                                         >  2:194642/1‑139 (MQ=255)                                                               aCGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCg                                                                 >  2:461512/1‑139 (MQ=255)                                                                      aaTCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAgg                                                                                                                                 <  1:1518266/68‑1 (MQ=255)                                                                      aaTCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAgg                                                                                                                                 >  2:1518266/1‑68 (MQ=255)                                                                            gggTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGc                                                                                                                            >  1:1233875/1‑67 (MQ=255)                                                                            gggTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGc                                                                                                                            <  2:1233875/67‑1 (MQ=255)                                                                             ggTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAacac                                                                               <  1:482328/111‑1 (MQ=255)                                                                             ggTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAacac                                                                               >  2:482328/1‑111 (MQ=255)                                                                             ggTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCgg                                                   >  2:714663/1‑139 (MQ=255)                                                                                 aCAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCaa                                               <  2:1446477/139‑1 (MQ=255)                                                                                    aTCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGtctgtctcttata                                            <  2:1359434/139‑13 (MQ=255)                                                                                       accaccGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACA‑CAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCaa                                        >  2:1778484/1‑139 (MQ=255)                                                                                           ccGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACg                                     <  1:461512/139‑1 (MQ=255)                                                                                           ccGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACg                                     <  2:626452/139‑1 (MQ=255)                                                                                             gcgAGGTTCAGCTCCACCACCAGGTCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgct                                   >  2:298194/1‑139 (MQ=255)                                                                                             gcgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATg                                                                   >  2:1756406/1‑107 (MQ=255)                                                                                             gcgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATg                                                                   <  1:1756406/107‑1 (MQ=255)                                                                                              cgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAgg                                                                                                                                 >  1:838733/1‑44 (MQ=255)                                                                                              cgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAgg                                                                                                                                 <  2:838733/44‑1 (MQ=255)                                                                                              cgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgctg                                  >  2:537996/1‑139 (MQ=255)                                                                                                aGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgctgct                                >  2:1299748/1‑139 (MQ=255)                                                                                                   ttCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCa                             >  1:1166565/1‑139 (MQ=255)                                                                                                             ccaccaGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAaca                                                                                <  1:930773/78‑1 (MQ=255)                                                                                                             ccaccaGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAaca                                                                                >  2:930773/1‑78 (MQ=255)                                                                                                                ccaGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACg                                     >  2:1461570/1‑118 (MQ=255)                                                                                                                ccaGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACg                                     <  1:1461570/118‑1 (MQ=255)                                                                                                                ccaGGCCGCGCAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaa                >  2:143331/1‑139 (MQ=255)                                                                                                                        cgcAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATAcc                         >  2:540063/1‑122 (MQ=255)                                                                                                                        cgcAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATAcc                         <  1:540063/122‑1 (MQ=255)                                                                                                                        cgcAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCg        >  1:821841/1‑139 (MQ=255)                                                                                                                        cgcAGCAGGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCACAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCg        >  1:664709/1‑139 (MQ=255)                                                                                                                           agcagGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCa     <  2:1476388/139‑1 (MQ=255)                                                                                                                              agGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCt                                            <  2:308953/97‑1 (MQ=255)                                                                                                                              agGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCt                                            >  1:308953/1‑97 (MQ=255)                                                                                                                              agGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGgt  <  1:719900/139‑1 (MQ=255)                                                                                                                              agGTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGgt  >  2:719900/1‑139 (MQ=255)                                                                                                                               ggTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACg                                     >  2:482373/1‑103 (MQ=255)                                                                                                                               ggTCGGGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACg                                     <  1:482373/103‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     ACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGT  >  NC_000913/4318080‑4318343

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