Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,507,942	A→G	Y217H (TAC→CAC)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,507,942	0	A	G	100.0%	115.5 / NA	33	Y217H (TAC→CAC)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/14);  total (19/14)

ATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTA  >  NC_000913/3507817‑3508076                                                                                                                              |                                                                                                                                       aTAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCCATACGCTTTTCGTCCGGGTGGTAACTGTTCTTg                                                                                                                           >  2:714818/1‑139 (MQ=255)  tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                          >  1:648054/1‑139 (MQ=255)    gCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCg                                                                                                                        >  2:133910/1‑139 (MQ=255)     ccATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                       <  2:648054/139‑1 (MQ=255)         aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                                   >  1:70240/1‑139 (MQ=255)             accGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                               <  2:464054/139‑1 (MQ=255)              ccGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGa                                                                                                              >  2:756882/1‑139 (MQ=255)               cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                             >  2:452045/1‑139 (MQ=255)                   ggCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACTCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGc                                                                                                         >  2:257280/1‑139 (MQ=255)                   ggCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGc                                                                                                         >  1:603832/1‑139 (MQ=255)                    gCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCa                                                                                                        >  2:925227/1‑139 (MQ=255)                                                   tCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAt                                                                         >  1:776836/1‑139 (MQ=255)                                                       cAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgttg                                                                     >  2:124321/1‑139 (MQ=255)                                                          cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                                  >  1:372775/1‑139 (MQ=255)                                                               cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                             <  1:133910/139‑1 (MQ=255)                                                                gcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGatcat                                                                                  <  2:823427/117‑1 (MQ=255)                                                                gcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGatcat                                                                                  >  1:823427/1‑117 (MQ=255)                                                                gcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTg                                                            <  1:714818/139‑1 (MQ=255)                                                                    ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCtt                                                        <  2:882149/139‑1 (MQ=255)                                                                        gCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAg                                                    <  2:603832/139‑1 (MQ=255)                                                                                      cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                       >  1:194527/1‑58 (MQ=255)                                                                                      cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                       <  2:194527/58‑1 (MQ=255)                                                                                      cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCaa                                      <  2:927605/139‑1 (MQ=255)                                                                                             cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCCCGTAcg                                                                                                  >  1:466650/1‑72 (MQ=255)                                                                                             cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTAcg                                                                                                  <  2:466650/72‑1 (MQ=255)                                                                                             cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAc                               >  1:827981/1‑139 (MQ=255)                                                                                                   gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                         <  1:756882/139‑1 (MQ=255)                                                                                                              cGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAAccc              >  2:662033/1‑139 (MQ=255)                                                                                                               gCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                                    >  2:228852/1‑52 (MQ=255)                                                                                                               gCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                                    <  1:228852/52‑1 (MQ=255)                                                                                                                       tCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCtt     <  2:827981/139‑1 (MQ=255)                                                                                                                         cGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGt   >  1:666994/1‑139 (MQ=255)                                                                                                                          gggTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTa  <  2:114151/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                       ATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTA  >  NC_000913/3507817‑3508076

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