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%	179.3 / NA	50	I5T (ATC→ACC)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/24);  total (26/24)

AGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGA  >  NC_000913/4325018‑4325268                                                                                                                                       |                                                                                                                     aGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAAt                                                                                                                  <  2:1205070/139‑1 (MQ=255)       gTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTg                                                                                                            >  2:60440/1‑139 (MQ=255)         aaaaGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCa                                                                                                          <  2:1338882/139‑1 (MQ=255)                tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                   >  2:220002/1‑139 (MQ=255)                  gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                 <  2:202988/139‑1 (MQ=255)                   aTTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGcc                                                                                                >  1:946427/1‑139 (MQ=255)                                  aCGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTgg                                                                                 >  2:516918/1‑139 (MQ=255)                                        cAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                                                                           <  1:831370/139‑1 (MQ=255)                                               agCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg                                                                    >  1:56159/1‑139 (MQ=255)                                                 ccaccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                                                                           >  1:1254022/1‑130 (MQ=255)                                                 ccaccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                                                                           <  2:1254022/130‑1 (MQ=255)                                                         tCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg                                                          <  1:202259/139‑1 (MQ=255)                                                               tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt                                                           >  1:1129315/1‑132 (MQ=255)                                                               tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt                                                           <  2:1129315/132‑1 (MQ=255)                                                               tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAtt                                                    >  2:7911/1‑139 (MQ=255)                                                                 atgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTg                                                  >  1:707442/1‑139 (MQ=255)                                                                  tgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                 <  2:707442/139‑1 (MQ=255)                                                                    aaTGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGtt                                               >  1:923535/1‑139 (MQ=255)                                                                       gTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaat                                            >  1:1208199/1‑139 (MQ=255)                                                                           aGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCCCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTaaa                                        >  2:920950/1‑139 (MQ=255)                                                                            gATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTaaaa                                       <  1:220002/139‑1 (MQ=255)                                                                                 aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                               >  1:883907/1‑78 (MQ=255)                                                                                 aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                               <  2:883907/78‑1 (MQ=255)                                                                                 aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAg                                                                                    >  2:98436/1‑89 (MQ=255)                                                                                 aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAg                                                                                    <  1:98436/89‑1 (MQ=255)                                                                                 aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                 >  1:452550/1‑124 (MQ=255)                                                                                 aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                 <  2:452550/124‑1 (MQ=255)                                                                                  ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                 >  2:190459/1‑139 (MQ=255)                                                                                  ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                 >  2:369381/1‑139 (MQ=255)                                                                                          aGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTGCTCGACACGGGTAATCGTTTGCATGATTCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACGAACAGATg                         <  2:436759/139‑1 (MQ=255)                                                                                            cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                       <  1:369381/139‑1 (MQ=255)                                                                                                 gATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                  >  1:69613/1‑91 (MQ=255)                                                                                                 gATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                  <  2:69613/91‑1 (MQ=255)                                                                                                  aTGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTa                 >  1:847228/1‑139 (MQ=255)                                                                                                       gATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg            <  1:190459/139‑1 (MQ=255)                                                                                                       gATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg            <  1:60440/139‑1 (MQ=255)                                                                                                        aTTGCAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                 >  2:783721/1‑53 (MQ=255)                                                                                                        aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                 <  1:783721/53‑1 (MQ=255)                                                                                                           gAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAAtt        >  1:162136/1‑139 (MQ=255)                                                                                                            aaGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTc       >  2:870361/1‑139 (MQ=255)                                                                                                               gTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGAt    <  2:923535/139‑1 (MQ=255)                                                                                                                 tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa  <  1:870361/139‑1 (MQ=255)                                                                                                                  ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                 <  2:350889/43‑1 (MQ=255)                                                                                                                  ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                 >  1:350889/1‑43 (MQ=255)                                                                                                                  ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg            <  2:661281/128‑1 (MQ=255)                                                                                                                  ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg            >  1:345891/1‑128 (MQ=255)                                                                                                                  ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg            <  2:345891/128‑1 (MQ=255)                                                                                                                  ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg            >  1:661281/1‑128 (MQ=255)                                                                                                                          cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                 >  1:943053/1‑83 (MQ=255)                                                                                                                          cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                 <  2:943053/83‑1 (MQ=255)                                                                                                                                       |                                                                                                                     AGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGA  >  NC_000913/4325018‑4325268

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