Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,324,474	A→G	V231A (GTC→GCC)	phnC ←	phosphonate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,324,474	0	A	G	100.0%	120.1 / NA	34	V231A (GTC→GCC)	phnC	phosphonate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/19);  total (15/19)

GCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCGCGCAGGGTGTCCATCACGATGCG  >  NC_000913/4324345‑4324601                                                                                                                                  |                                                                                                                                gCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTg                                                                                                                       >  2:170796/1‑140 (MQ=255)      cATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGgcg                                                                                                                    >  1:253803/1‑138 (MQ=255)      cATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGgcg                                                                                                                    <  2:253803/138‑1 (MQ=255)      cATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCa                                                                                                                  >  2:953508/1‑140 (MQ=255)             tttCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGTTGCCGTCTTAGTAGACGTGCCACTGGCGCAGGGcgac                                                                                                           >  1:424724/1‑140 (MQ=255)              ttCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGcgacg                                                                                                          >  2:1020789/1‑140 (MQ=255)                 cTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATg                                                                                                       >  2:83385/1‑140 (MQ=255)                     atgatgGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGtt                                                                                                   <  2:151525/140‑1 (MQ=255)                         tgGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCa                                                                                               <  1:925373/140‑1 (MQ=255)                           gggATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGt                                                                                             >  2:759476/1‑140 (MQ=255)                            ggATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAAATGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGAGGCGTTGGCAGTa                                                                                            >  2:633764/1‑140 (MQ=255)                             gATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGA‑‑‑GGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcgc                                                                                        <  2:206844/140‑1 (MQ=255)                              aTGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgc                                                                                          >  2:240840/1‑140 (MQ=255)                               tGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGAGCCCCTGGCGCAGGGCGACGAGGCGTTCGCAGTAgcg                                                                                         >  1:323801/1‑140 (MQ=255)                                gTCAGGCAGCTTTCGTGTTCTCCTCCACGAGGTTAATACTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcgc                                                                                        <  1:326384/140‑1 (MQ=255)                                gTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcgc                                                                                        <  2:886379/140‑1 (MQ=255)                                gTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcgc                                                                                        >  2:927019/1‑140 (MQ=255)                                 tCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGTCCCTGGCTCAGGGCCACGATGCGTTCGCAGTAGCGCa                                                                                       >  1:452788/1‑140 (MQ=255)                                 tCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCa                                                                                       <  1:1020789/140‑1 (MQ=255)                                   aGGCAGCTATCGCGTTTTCTTCAACGCGGTTTAAGATGCTGTAGAGATAGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAgg                                                                                     <  2:452788/140‑1 (MQ=255)                                    ggCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAggg                                                                                    <  2:512921/140‑1 (MQ=255)                                    ggCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAggg                                                                                    <  1:170796/140‑1 (MQ=255)                                             cgcgTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcgc                                                                                        <  2:1000806/127‑1 (MQ=255)                                             cgcgTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcgc                                                                                        >  1:1000806/1‑127 (MQ=255)                                                          aCGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGt                                                              <  1:876770/140‑1 (MQ=255)                                                          aCGCGGTTAATGCTGCGCTAGAAATGGTCAAAACGTTCGGTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGt                                                              <  2:426509/140‑1 (MQ=255)                                                             cgGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTgac                                                           <  1:83385/140‑1 (MQ=255)                                                                ttAATGCTGCGGTAGAGTTGGTCAAAACGTTGGTTCTCATACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTgacgac                                                        <  2:315267/140‑1 (MQ=255)                                                                          ggTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGGGTAATCCGCCTTATCCAGCGTAACGACCACGGTGATg                                              >  1:484376/1‑140 (MQ=255)                                                                                   ggTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGTGCGACGATGCGCTCGCCAAAGCGCAGGGCGTAATCCACCTGATGCGGCGTGACGACCACGGTGATGCCGTCGTTc                                     >  2:979280/1‑140 (MQ=255)                                                                                                 ttATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTcgcg                       <  2:75716/140‑1 (MQ=255)                                                                                                  tATCAAACTGTTGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCgcgc                      <  1:953508/140‑1 (MQ=255)                                                                                                            ttGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCGCGCAgg                   <  1:851997/133‑1 (MQ=255)                                                                                                            ttGGCTGCTGCCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCGCGCAgg                   >  2:851997/1‑133 (MQ=255)                                                                                                                 tgctgcCGTCGTAGAAGGCGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCGCGCAGGGTGTCCATCACg       >  2:562853/1‑140 (MQ=255)                                                                                                                      ccGTCGTAGAAGGCGTGCCCCTAGCGCAGGGTGACGATGCGTTCGCAGGAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCGCGCAGGGTGTCCATCACGATgcg  <  2:424724/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                GCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTCGTTCTGGTTGATGTCGCGCAGGGTGTCCATCACGATGCG  >  NC_000913/4324345‑4324601

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