Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,152	A→G	I5T (ATC→ACC)	phnC ←	phosphonate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,152	0	A	G	100.0%	125.7 / NA	36	I5T (ATC→ACC)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/18);  total (18/18)

TAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACA  >  NC_000913/4325025‑4325270                                                                                                                                |                                                                                                                        tAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                              <  2:317307/139‑1 (MQ=255)         tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                  <  2:478461/127‑1 (MQ=255)         tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                  >  1:478461/1‑127 (MQ=255)          ggATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAg                                                                                                     >  2:340899/1‑139 (MQ=255)           gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                    <  1:388585/139‑1 (MQ=255)                                     cAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTg                                                                         <  1:841976/139‑1 (MQ=255)                                     cAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTg                                                                         <  1:723032/139‑1 (MQ=255)                                         gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCgg                                                                     >  2:720826/1‑139 (MQ=255)                                             ccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                              <  1:425805/95‑1 (MQ=255)                                             ccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                              >  2:425805/1‑95 (MQ=255)                                                tttCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                       >  1:117621/1‑83 (MQ=255)                                                tttCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                       <  2:117621/83‑1 (MQ=255)                                                        tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTACTCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAtt                                                      >  1:448644/1‑139 (MQ=255)                                                           tgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                   <  1:340899/139‑1 (MQ=255)                                                                    aGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                              <  1:279164/72‑1 (MQ=255)                                                                    aGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                              >  2:279164/1‑72 (MQ=255)                                                                      aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑aaaaaa                                                                                       <  1:834000/92‑1 (MQ=255)                                                                      aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑aaaaaa                                                                                       >  2:834000/1‑92 (MQ=255)                                                                           ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCAt                                                                             >  1:308021/1‑97 (MQ=255)                                                                           ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCAt                                                                             <  2:308021/97‑1 (MQ=255)                                                                           ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                   <  1:488751/123‑1 (MQ=255)                                                                           ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                   >  2:488751/1‑123 (MQ=255)                                                                           ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                   >  2:892990/1‑139 (MQ=255)                                                                                     cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                         <  1:441556/139‑1 (MQ=255)                                                                                     cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                         >  1:268786/1‑139 (MQ=255)                                                                                     cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTAGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                         >  1:646031/1‑139 (MQ=255)                                                                                               tGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAt                                                       <  2:216437/99‑1 (MQ=255)                                                                                               tGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAt                                                       >  1:216437/1‑99 (MQ=255)                                                                                                   tGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg              >  1:511971/1‑136 (MQ=255)                                                                                                   tGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg              <  2:511971/136‑1 (MQ=255)                                                                                                     aaGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTc         >  2:214196/1‑139 (MQ=255)                                                                                                      aGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                         <  1:142075/123‑1 (MQ=255)                                                                                                      aGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAAGTGGCTCCATCGTGGCGGATTATGGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                         >  2:142075/1‑123 (MQ=255)                                                                                                            tGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACa  <  2:268786/139‑1 (MQ=255)                                                                                                                  gCTTCTCGACACGGGTAATTGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                   >  2:33898/1‑100 (MQ=255)                                                                                                                  gCTTCTCGACACGGGTAATTGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                   <  1:33898/100‑1 (MQ=255)                                                                                                                                |                                                                                                                        TAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACA  >  NC_000913/4325025‑4325270

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