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

TAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGT  >  NC_000913/3507818‑3508067                                                                                                                             |                                                                                                                              tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                 >  1:505722/1‑139 (MQ=255) tAGCCATAACCACCGTTAGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                 >  1:1056964/1‑139 (MQ=255)        aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGcc                                                                                                                >  2:992782/1‑133 (MQ=255)        aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGcc                                                                                                                <  1:992782/133‑1 (MQ=255)        aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                          >  1:915249/1‑139 (MQ=255)            accGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                      <  2:915249/139‑1 (MQ=255)            accGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                      >  2:115298/1‑139 (MQ=255)              cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                    >  2:1386585/1‑139 (MQ=255)              cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                    >  2:592107/1‑139 (MQ=255)              cGTTGGCTTTTAAATGGGAGCGGCGCGTAAT‑TCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAAc                                                                                                   >  2:401389/1‑139 (MQ=255)               gTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAAc                                                                                                   >  1:1387598/1‑139 (MQ=255)                  ggCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGc                                                                                                <  2:505722/139‑1 (MQ=255)                  ggCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGGTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGc                                                                                                >  2:1477659/1‑139 (MQ=255)                       ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTATCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTATCGAACTGCACGTa                                                                                           <  1:1028208/139‑1 (MQ=255)                       ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGtt                                                                                           <  2:1056964/139‑2 (MQ=255)                       ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           <  2:563796/139‑1 (MQ=255)                       ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           <  2:1125142/139‑1 (MQ=255)                               gAGCGGCGCGTAATATCGATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAATCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCga                                                                                   >  1:1057026/1‑139 (MQ=255)                                 gcggcgCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCgaga                                                                                 >  2:251067/1‑139 (MQ=255)                                  cggcgCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgat                                                                                >  2:1413993/1‑139 (MQ=255)                                      gcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgat                                                                            <  1:1386585/139‑1 (MQ=255)                                       cgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                           <  1:115298/139‑1 (MQ=255)                                                 aTCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaaa                                                                 >  2:1203290/1‑139 (MQ=255)                                                 aTCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaaa                                                                 >  2:1654044/1‑139 (MQ=255)                                                  tCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                 >  2:1406945/1‑90 (MQ=255)                                                  tCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                 <  1:1406945/90‑1 (MQ=255)                                                    gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgt                                                              <  1:401389/139‑1 (MQ=255)                                                      cAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           >  1:410146/1‑108 (MQ=255)                                                      cAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           <  2:410146/108‑1 (MQ=255)                                                       aGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgttgt                                                           >  1:865396/1‑139 (MQ=255)                                                        gCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           <  2:827627/106‑1 (MQ=255)                                                        gCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           >  1:827627/1‑106 (MQ=255)                                                         cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                      >  2:1155167/1‑78 (MQ=255)                                                         cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                      <  1:1155167/78‑1 (MQ=255)                                                         cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                         >  1:311635/1‑139 (MQ=255)                                                              cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                    <  1:1654044/139‑1 (MQ=255)                                                                   ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                            <  1:612170/78‑1 (MQ=255)                                                                   ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                            >  2:612170/1‑78 (MQ=255)                                                                   ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgtgt                                                                                                        >  2:1362729/1‑82 (MQ=255)                                                                   ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgtgt                                                                                                        <  1:1362729/82‑1 (MQ=255)                                                                        cAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                            >  1:268886/1‑73 (MQ=255)                                                                        cAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                            <  2:268886/73‑1 (MQ=255)                                                                         aaCCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGAt                                         <  2:1387598/139‑1 (MQ=255)                                                                             cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                     <  1:1203290/139‑1 (MQ=255)                                                                                ggTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                      >  2:177733/1‑71 (MQ=255)                                                                                ggTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                      <  1:177733/71‑1 (MQ=255)                                                                                ggTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGtc                                  <  2:1057026/139‑2 (MQ=255)                                                                                ggTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTg                                  <  1:592107/139‑1 (MQ=255)                                                                                 gTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGa                                 <  2:865396/139‑1 (MQ=255)                                                                                     cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                              >  2:954164/1‑58 (MQ=255)                                                                                     cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                              <  1:954164/58‑1 (MQ=255)                                                                                     cgTAGCGCATGAAGCATCGCAATACGCTTCTCGTCCGGGTGGTAATTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCaa                             >  2:977995/1‑139 (MQ=255)                                                                                      gTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTg                                                                                                 >  2:1639988/1‑70 (MQ=255)                                                                                      gTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTg                                                                                                 <  1:1639988/70‑1 (MQ=255)                                                                                       tAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           >  2:1622553/1‑75 (MQ=255)                                                                                       tAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                           <  1:1622553/75‑1 (MQ=255)                                                                                                  gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                <  1:251067/139‑1 (MQ=255)                                                                                                  gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                <  1:378753/139‑1 (MQ=255)                                                                                                    aTCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCtt                                               >  2:45966/1‑106 (MQ=255)                                                                                                    aTCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCtt                                               <  1:45966/106‑1 (MQ=255)                                                                                                      cGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGatca                                                                          >  1:710671/1‑77 (MQ=255)                                                                                                      cGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGatca                                                                          <  2:710671/77‑1 (MQ=255)                                                                                                           tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTAcgcg                                                                                       <  2:1128737/59‑1 (MQ=255)                                                                                                           tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTAcgcg                                                                                       >  1:1128737/1‑59 (MQ=255)                                                                                                           tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAAc       >  2:397312/1‑139 (MQ=255)                                                                                                                ttttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGt  >  1:1378125/1‑139 (MQ=255)                                                                                                                   tCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                      <  1:1487125/36‑1 (MQ=255)                                                                                                                   tCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                      >  2:1487125/1‑36 (MQ=255)                                                                                                                             |                                                                                                                              TAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGT  >  NC_000913/3507818‑3508067

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