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

AGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGG  >  NC_000913/3507813‑3508069                                                                                                                                   |                                                                                                                                aGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                        >  2:699192/1‑139 (MQ=255)   tCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCt                                                                                                                      >  1:1254567/1‑139 (MQ=255)     aTAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTg                                                                                                                    >  1:328024/1‑139 (MQ=255)      tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                   >  1:916859/1‑139 (MQ=255)         ccATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                <  1:1295329/139‑1 (MQ=255)             aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                            >  1:1137453/1‑139 (MQ=255)             aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                            >  2:910325/1‑139 (MQ=255)             aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                            >  2:706252/1‑139 (MQ=255)              accaccGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgtg                                                                                                           <  1:706252/139‑1 (MQ=255)                 accGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                        <  2:328024/139‑1 (MQ=255)                 accGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCCATGGTGTCg                                                                                                        >  1:832160/1‑139 (MQ=255)                   cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                              >  1:510383/1‑131 (MQ=255)                   cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                              <  2:510383/131‑1 (MQ=255)                   cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                      >  1:239429/1‑139 (MQ=255)                     ttGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACt                                                                                                    >  1:16935/1‑139 (MQ=255)                        gCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGGCGATGGTGTCGAACTGCa                                                                                                 >  1:871233/1‑139 (MQ=255)                        gCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCa                                                                                                 <  2:239429/139‑1 (MQ=255)                          ttttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACg                                                                                               <  2:1254567/139‑1 (MQ=255)                            ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             <  2:452387/139‑1 (MQ=255)                            ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             <  2:1291745/139‑1 (MQ=255)                                 tGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACgcgc                                                                                        <  2:16935/139‑1 (MQ=255)                                  gggAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGcc                                                                                       >  2:664987/1‑139 (MQ=255)                                    gAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCga                                                                                     >  2:496734/1‑139 (MQ=255)                                       cggcgCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgat                                                                                  <  1:699192/139‑1 (MQ=255)                                          cgcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                <  1:205891/106‑1 (MQ=255)                                          cgcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                >  2:205891/1‑106 (MQ=255)                                            cgtaagagtcATCGAGAGCATGACCCGGTTCAGCAACACACTGTCGCGGAGCACATGAAAGCA‑CCCAATACGCTTTTCGTACCGGGGCTAACTGTTCTTGCCGATAGTGTCAAGATTCACGTACGCGCCGAgatcgagc                                                                             <  2:776136/130‑3 (MQ=255)                                            cgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                             <  1:496734/139‑1 (MQ=255)                                                      aTCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                             <  2:713814/81‑1 (MQ=255)                                                      aTCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                             >  1:713814/1‑81 (MQ=255)                                                                 gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                        <  1:664987/139‑1 (MQ=255)                                                                 gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                        >  2:780590/1‑139 (MQ=255)                                                                   cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                      <  1:910325/139‑1 (MQ=255)                                                                   cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                      >  1:446331/1‑139 (MQ=255)                                                                        ttCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             <  2:1082978/95‑1 (MQ=255)                                                                        ttCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             >  1:1082978/1‑95 (MQ=255)                                                                             cAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGa                                            <  2:871233/139‑1 (MQ=255)                                                                                cccACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCg                                         >  1:309580/1‑139 (MQ=255)                                                                                  cACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                       <  1:155293/139‑1 (MQ=255)                                                                                     ggTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTTGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                        <  1:260213/119‑1 (MQ=255)                                                                                     ggTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTTGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                        >  2:260213/1‑119 (MQ=255)                                                                                          cgTAGCGCATG‑AAGCATCGTAATACGCTTTTCGTCCGTGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                      <  2:1274807/68‑1 (MQ=255)                                                                                          cgTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                      >  1:1274807/1‑68 (MQ=255)                                                                                          cgTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             <  2:402164/77‑1 (MQ=255)                                                                                          cgTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             >  1:402164/1‑77 (MQ=255)                                                                                          cgTAGCGCATG‑AAGCATCGCAATACGCTTTTCGACCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAATTGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCaa                               >  1:283582/1‑139 (MQ=255)                                                                                            tAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACg                             >  1:1118913/1‑139 (MQ=255)                                                                                                 cATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAc                        >  1:154250/1‑139 (MQ=255)                                                                                                  aTG‑AAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAcg                       <  1:538821/139‑1 (MQ=255)                                                                                                      aaGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGGTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAcgcgc                    <  2:916859/139‑1 (MQ=255)                                                                                                        gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                  <  2:309580/139‑1 (MQ=255)                                                                                                           tCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             <  1:347201/61‑1 (MQ=255)                                                                                                           tCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             >  2:347201/1‑61 (MQ=255)                                                                                                           tCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAg               >  2:613469/1‑139 (MQ=255)                                                                                                                 tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                        <  2:159236/44‑1 (MQ=255)                                                                                                                 tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                        >  1:159236/1‑44 (MQ=255)                                                                                                                 tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCtt                                                 <  2:459672/99‑1 (MQ=255)                                                                                                                 tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCtt                                                 >  1:459672/1‑99 (MQ=255)                                                                                                                 tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAAc         >  1:511104/1‑139 (MQ=255)                                                                                                                      ttttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                             >  2:1229815/1‑66 (MQ=255)                                                                                                                      ttttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                             <  1:1229815/66‑1 (MQ=255)                                                                                                                      ttttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                        <  2:505985/87‑1 (MQ=255)                                                                                                                      ttttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                        >  1:505985/1‑87 (MQ=255)                                                                                                                        ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgat                                                                                  <  1:827934/59‑1 (MQ=255)                                                                                                                        ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgat                                                                                  >  2:827934/1‑59 (MQ=255)                                                                                                                        ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTgg  >  1:924120/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                AGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATG‑AAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGG  >  NC_000913/3507813‑3508069

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