Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,507,955	G→A	G212G (GGC→GGT)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,507,955	0	G	A	100.0%	193.1 / NA	50	G212G (GGC→GGT)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/27);  total (23/27)

AAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGC  >  NC_000913/3507809‑3508096                                                                                                                                                   |                                                                                                                                              aaaTAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCg                                                                                                                                            >  2:144427/1‑150 (MQ=255) aaaTAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCg                                                                                                                                            >  1:240712/1‑150 (MQ=255)   aTAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATg                                                                                                                                         >  2:24000/1‑151 (MQ=255)     aGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATgg                                                                                                                                        <  2:477944/150‑1 (MQ=255)     aGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGgt                                                                                                                                       >  2:351093/1‑151 (MQ=255)      gTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGgtg                                                                                                                                      <  2:427487/151‑1 (MQ=255)               aTAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGc                                                                                                                             >  1:125679/1‑151 (MQ=255)                tAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGc                                                                                                                             >  1:435660/1‑150 (MQ=255)                 aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCAc                                                                                                                           <  2:331282/151‑1 (MQ=255)                 aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCAc                                                                                                                           <  2:416847/151‑1 (MQ=255)                          tGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACgcgc                                                                                                                   <  2:134118/150‑1 (MQ=255)                            gCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCg                                                                                                                 >  1:681681/1‑150 (MQ=255)                               tttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCgaga                                                                                                              <  1:254112/150‑1 (MQ=255)                                   aaTGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAgatc                                                                                                            >  1:394206/1‑148 (MQ=255)                                    aTGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                                                        >  2:274854/1‑151 (MQ=255)                                    aTGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                                                        <  1:134166/151‑1 (MQ=255)                                         aGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCtt                                                                                                   >  1:4378/1‑151 (MQ=255)                                         aGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCtt                                                                                                   >  1:19531/1‑151 (MQ=255)                                             gcgcGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaa                                                                                               <  1:274780/151‑1 (MQ=255)                                               gcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaa                                                                                               <  2:311760/149‑1 (MQ=255)                                                 gTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgt                                                                                           >  2:359267/1‑151 (MQ=255)                                                   aaTATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgttg                                                                                         <  2:644366/151‑1 (MQ=255)                                                   aaTATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtgttg                                                                                         <  2:148681/151‑1 (MQ=255)                                                       tCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                                                      >  1:637763/1‑150 (MQ=255)                                                        ccATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                                                      >  1:65487/1‑149 (MQ=255)                                                        ccATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAg                                                                                    <  1:351082/151‑1 (MQ=255)                                                            cGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                                                 <  2:226516/150‑1 (MQ=255)                                                             gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGt                                                                               <  1:589510/151‑1 (MQ=255)                                                                  cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTc                                                                          >  2:674796/1‑151 (MQ=255)                                                                   aTGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCa                                                                         <  2:514434/151‑1 (MQ=255)                                                                     gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAg                                                                        <  2:346824/150‑1 (MQ=255)                                                                     gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGa                                                                       <  2:329609/151‑1 (MQ=255)                                                                      aCGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGAt                                                                      <  2:562565/151‑1 (MQ=255)                                                                        gcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCg                                                                    <  2:125681/151‑1 (MQ=255)                                                                         cgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGc                                                                   <  2:445126/151‑1 (MQ=255)                                                                                gCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGAcc                                                            >  2:432424/1‑151 (MQ=255)                                                                                   aCCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAAc                                                         <  2:589424/151‑1 (MQ=255)                                                                                     ccACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACgg                                                       <  1:183856/151‑1 (MQ=255)                                                                                       aCGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTg                                                     <  1:191635/151‑1 (MQ=255)                                                                                         ggTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGCGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGCACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCACCGGt                                                      >  1:658375/1‑148 (MQ=255)                                                                                                tAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                                             >  1:180637/1‑150 (MQ=255)                                                                                                 aGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGaaa                                           <  1:50411/151‑1 (MQ=255)                                                                                                      aTGAAGCATCGCAATACGCTTTTCGTCCGGGAAGTAACTGTTCTTACAGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATc                                       >  1:492604/1‑150 (MQ=255)                                                                                                       tGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCa                                      >  2:341071/1‑150 (MQ=255)                                                                                                           gCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAAcccc                                 <  1:179790/151‑1 (MQ=255)                                                                                                                         ttttCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGc                    >  1:396892/1‑150 (MQ=255)                                                                                                                            tCGTCCGGGTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAAc                 <  2:103124/150‑1 (MQ=255)                                                                                                                                  gggTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGtt          <  1:496310/151‑1 (MQ=255)                                                                                                                                  gggTAGTAACTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGCGCTTGTAGCAACGCCAGtt          >  1:495490/1‑151 (MQ=255)                                                                                                                                          aCTGTTCTTACCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGc  >  2:506601/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                              AAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGC  >  NC_000913/3507809‑3508096

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