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_0009133,507,9420AG98.4% 210.8 / NA 61Y217H (TAC→CAC) phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/33);  minor base T (0/1);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGTAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGG  >  NC_000913/3507806‑3508069
                                                                                                                                        |                                                                                                                               
aGTAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGggtggt                                                                                                                               >  1:647597/1‑139 (MQ=255)
    aaTAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAAc                                                                                                                            <  2:189817/138‑1 (MQ=255)
    aaTAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAAc                                                                                                                            >  1:189817/1‑138 (MQ=255)
       aGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTCGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                        >  2:998615/1‑139 (MQ=255)
       aGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                        <  2:1116980/139‑1 (MQ=255)
       aGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                        <  1:686764/139‑1 (MQ=255)
           aTAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTg                                                                                                                    <  2:647597/139‑1 (MQ=255)
            tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                   >  1:1024572/1‑139 (MQ=255)
            tAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                   >  2:1088558/1‑139 (MQ=255)
               ccATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                <  2:1332311/139‑1 (MQ=255)
                 aTAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                              <  1:1088558/139‑1 (MQ=255)
                   aaCCACCGTTGGCTTTTAAATGGGAGTGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                            >  2:303580/1‑139 (MQ=255)
                        ccGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGa                                                                                                       <  1:1141641/139‑1 (MQ=255)
                         cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                              <  2:879292/131‑1 (MQ=255)
                         cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                      >  1:562835/1‑139 (MQ=255)
                         cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                      >  1:563826/1‑139 (MQ=255)
                         cGTGGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATg                                                                                                              >  1:879292/1‑131 (MQ=255)
                                         ggAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCg                                                                                      >  2:70833/1‑139 (MQ=255)
                                           aGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCgag                                                                                    >  2:478888/1‑139 (MQ=255)
                                                 gcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgat                                                                              <  2:563826/139‑1 (MQ=255)
                                                    tAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGatcat                                                                           >  1:1189360/1‑139 (MQ=255)
                                                    tAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGatcat                                                                           <  1:70833/139‑1 (MQ=255)
                                                             tCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAt                                                                  <  2:1189360/139‑1 (MQ=255)
                                                               gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                >  1:314535/1‑91 (MQ=255)
                                                               gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                <  2:314535/91‑1 (MQ=255)
                                                               gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTAcgcg                                                                                         <  2:302582/114‑1 (MQ=255)
                                                               gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTAcgcg                                                                                         >  1:302582/1‑114 (MQ=255)
                                                                   gCATGACGCGGTTCAGCAACCCACGGTCGCTTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCTGAGATCGATCATCTTCaaaa                                                                   <  2:87362/132‑1 (MQ=255)
                                                                   gCATGACGCGGTTCAGCAACCCACGGTCGCTTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaaa                                                                   >  1:87362/1‑132 (MQ=255)
                                                                    cATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                           <  2:562835/139‑1 (MQ=255)
                                                                        aCGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGt                                                       <  1:704051/139‑1 (MQ=255)
                                                                         ggcgGTTGAGTAACCGACGGTCGGGTAGCGCATGAATCATCGCAAGACGCTTTTCGTCCGGGTTGTAACTGTTCTTGCCGATCGTGTCTAACGGTACGGACGCGCCGAGATCGATCAGCTTCAAAATGTTGTCGAGGtt                                                      <  1:998615/138‑1 (MQ=255)
                                                                         cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                        >  1:429394/1‑89 (MQ=255)
                                                                         cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                        <  2:429394/89‑1 (MQ=255)
                                                                         cgcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                      <  2:180641/139‑1 (MQ=255)
                                                                          gcgGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTg                                                     <  2:1024572/139‑1 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                        >  2:584325/1‑58 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                        <  1:584325/58‑1 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                <  1:1265003/66‑1 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                >  2:1265003/1‑66 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             >  1:1184060/1‑85 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                             <  2:1184060/85‑1 (MQ=255)
                                                                                        cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                       <  2:370803/139‑1 (MQ=255)
                                                                                           ggTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTg                                    <  1:478888/139‑1 (MQ=255)
                                                                                            gTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGa                                   >  1:410968/1‑139 (MQ=255)
                                                                                                cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCaa                               >  2:146903/1‑139 (MQ=255)
                                                                                                      gcATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGa                         <  2:119221/139‑1 (MQ=255)
                                                                                                       cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAc                        >  2:440819/1‑139 (MQ=255)
                                                                                                        aTGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAcg                       <  1:1154705/139‑1 (MQ=255)
                                                                                                             gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                  <  2:410968/139‑1 (MQ=255)
                                                                                                             gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                  <  2:30758/139‑1 (MQ=255)
                                                                                                             gCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                  <  1:1144199/139‑1 (MQ=255)
                                                                                                              cATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGaa                 >  1:169208/1‑139 (MQ=255)
                                                                                                               aacGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAcgcgcg                   >  1:715829/3‑136 (MQ=255)
                                                                                                               aacGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAcgcgcg                   <  2:715829/134‑1 (MQ=255)
                                                                                                               aTCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACg                                                                                               >  1:597823/1‑60 (MQ=255)
                                                                                                               aTCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACg                                                                                               <  2:597823/60‑1 (MQ=255)
                                                                                                                tCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                  >  1:1339156/1‑136 (MQ=255)
                                                                                                                tCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGa                  <  2:1339156/136‑1 (MQ=255)
                                                                                                                             ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTgg  >  2:41441/1‑139 (MQ=255)
                                                                                                                              tCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                       >  1:406540/1‑101 (MQ=255)
                                                                                                                              tCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                       <  2:406540/101‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
AGTAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGG  >  NC_000913/3507806‑3508069

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