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%	170.7 / NA	48	Y217H (TAC→CAC)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

TAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGA  >  NC_000913/3507808‑3508057                                                                                                                                       |                                                                                                                    tAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                 >  1:45118/1‑139 (MQ=255)          aTAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTg                                                                                                        >  1:706589/1‑139 (MQ=255)           tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                       >  1:190505/1‑139 (MQ=255)           tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                       >  1:62171/1‑139 (MQ=255)                  aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                >  2:12427/1‑139 (MQ=255)                      accGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                            >  1:933995/1‑139 (MQ=255)                       ccGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGa                                                                                           >  2:699116/1‑139 (MQ=255)                        cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                            <  2:612585/121‑1 (MQ=255)                        cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                            >  1:612585/1‑121 (MQ=255)                        cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                          >  2:397887/1‑139 (MQ=255)                        cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                          >  2:16745/1‑139 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATTGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTc                                                                                             >  1:611931/1‑127 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTc                                                                                             <  2:611931/127‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                 <  2:933995/139‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                 <  2:706589/139‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                 <  2:62171/139‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                 <  1:397887/139‑1 (MQ=255)                                         gAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCga                                                                         <  2:45118/139‑1 (MQ=255)                                         gAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCga                                                                         <  2:160022/139‑1 (MQ=255)                                               cgcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcga                                                                   <  1:835397/139‑1 (MQ=255)                                                   tAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCgag                                                                        <  1:638301/130‑1 (MQ=255)                                                   tAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCgag                                                                        >  2:638301/1‑130 (MQ=255)                                                        tCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                    >  2:374654/1‑97 (MQ=255)                                                        tCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                    <  1:374654/97‑1 (MQ=255)                                                              gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgt                                                    >  1:328474/1‑139 (MQ=255)                                                                   cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                          <  2:805665/96‑1 (MQ=255)                                                                   cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGCCGaa                                                                                          >  1:805665/1‑96 (MQ=255)                                                                    aTGACGCGGTTCAGCAACCCACAGTCGCGTAGCGCATGAAGTATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGa                                              <  2:524559/139‑1 (MQ=255)                                                                      gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                            <  2:328474/139‑1 (MQ=255)                                                                      gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                            <  1:12427/139‑1 (MQ=255)                                                                        cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                          <  1:699116/139‑1 (MQ=255)                                                                             ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                            >  1:543681/1‑132 (MQ=255)                                                                             ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                            <  2:543681/132‑1 (MQ=255)                                                                                  cAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGa                                <  1:394059/139‑1 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTACGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                           >  1:390741/1‑139 (MQ=255)                                                                                               cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                          >  2:345763/1‑68 (MQ=255)                                                                                               cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                          <  1:345763/68‑1 (MQ=255)                                                                                                      cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                               <  1:595635/104‑1 (MQ=255)                                                                                                      cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                               >  2:595635/1‑104 (MQ=255)                                                                                                      cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGa                                >  2:504319/1‑119 (MQ=255)                                                                                                      cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGa                                <  1:504319/119‑1 (MQ=255)                                                                                                              aTCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                            <  2:589636/51‑1 (MQ=255)                                                                                                              aTCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                            >  1:589636/1‑51 (MQ=255)                                                                                                                cGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGa  >  1:178244/1‑139 (MQ=255)                                                                                                                     tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTg                        <  2:685756/112‑1 (MQ=255)                                                                                                                     tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTg                        >  1:685756/1‑112 (MQ=255)                                                                                                                       cGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGAc                      >  2:147373/1‑112 (MQ=255)                                                                                                                       cGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGAc                      <  1:147373/112‑1 (MQ=255)                                                                                                                                       |                                                                                                                    TAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGA  >  NC_000913/3507808‑3508057

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