Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,508,454 Δ1 bp coding (137/879 nt) php ← phosphotriesterase homology protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,4510C.98.3% 239.7 / ‑2.4 58coding (140/879 nt)phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (28/29);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

ATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508315‑3508583
                                                                                                                                        |                                                                                                                                    
aTGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    <  1:256602/139‑1 (MQ=255)
 tGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                   <  1:424995/139‑2 (MQ=255)
  gTTCCGGGAAAAGCGCGTCTTGTTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGTATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCTCGcccgg                                                                                                                                  <  2:727626/139‑3 (MQ=255)
                          aataaCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                         <  1:92894/139‑1 (MQ=255)
                              aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  1:574392/139‑1 (MQ=255)
                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:596292/139‑1 (MQ=255)
                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:586556/139‑1 (MQ=255)
                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:294976/139‑1 (MQ=255)
                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  1:305010/139‑1 (MQ=255)
                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  2:32934/1‑139 (MQ=255)
                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           <  1:32934/139‑1 (MQ=255)
                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:27060/1‑124 (MQ=255)
                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:27060/124‑1 (MQ=255)
                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      >  2:795441/1‑139 (MQ=255)
                                                  tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:1017558/119‑1 (MQ=255)
                                                  tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:1017558/1‑119 (MQ=255)
                                                  tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                 >  2:211251/1‑139 (MQ=255)
                                                  tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                 >  1:830780/1‑139 (MQ=255)
                                                       cccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccggg                                                                                                                                 <  1:957479/87‑4 (MQ=255)
                                                       cccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCggg                                                                                                                                 >  2:957479/1‑87 (MQ=255)
                                                              cgcgCATTACATCAAGCATAAATTTCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAAGGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                     <  2:803115/139‑1 (MQ=255)
                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  1:1078066/92‑1 (MQ=255)
                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  2:1078066/1‑92 (MQ=255)
                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      >  2:288110/1‑107 (MQ=255)
                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      <  1:288110/107‑1 (MQ=255)
                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                      >  1:223308/1‑139 (MQ=255)
                                                                              cATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCgg                                                     <  1:37317/139‑1 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  1:1037546/93‑1 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  2:1037546/1‑93 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  1:196579/1‑139 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  2:196579/139‑1 (MQ=255)
                                                                                      gcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                      >  1:629497/1‑130 (MQ=255)
                                                                                      gcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                      <  2:629497/130‑1 (MQ=255)
                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcc                                                                                                                                     >  1:657822/1‑51 (MQ=255)
                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcc                                                                                                                                     <  2:657822/51‑1 (MQ=255)
                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:454550/136‑1 (MQ=255)
                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               >  2:454550/1‑136 (MQ=255)
                                                                                        gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCTGGTCATCAGGt                                                                                                                       <  2:644455/63‑1 (MQ=255)
                                                                                        gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGt                                                                                                                       >  1:644455/1‑63 (MQ=255)
                                                                                        gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:294580/135‑1 (MQ=255)
                                                                                        gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               >  2:294580/1‑135 (MQ=255)
                                                                                              cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGtcat                                                                                                                            >  1:802904/1‑52 (MQ=255)
                                                                                              cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGtcat                                                                                                                            <  2:802904/52‑1 (MQ=255)
                                                                                              cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTcc                                                                        >  2:1073122/1‑104 (MQ=255)
                                                                                              cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTcc                                                                        <  1:1073122/104‑1 (MQ=255)
                                                                                              cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAg                                     <  2:223308/139‑1 (MQ=255)
                                                                                                 cccATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATg                                  <  2:830780/139‑1 (MQ=255)
                                                                                                 cccATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATg                                  >  2:609479/1‑139 (MQ=255)
                                                                                                             ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAg                      >  1:72357/1‑139 (MQ=255)
                                                                                                             ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAg                      >  2:273234/1‑139 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCAt                             <  2:1010979/130‑1 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCAt                             >  1:1010979/1‑130 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        >  1:1042810/1‑135 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:215733/135‑1 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        >  2:215733/1‑135 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  2:1042810/135‑1 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  1:740474/1‑139 (MQ=255)
                                                                                                                          tCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                              >  2:496290/1‑86 (MQ=255)
                                                                                                                          tCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                              <  1:496290/86‑1 (MQ=255)
                                                                                                                          tCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                              <  2:77663/86‑1 (MQ=255)
                                                                                                                          tCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                              >  1:77663/1‑86 (MQ=255)
                                                                                                                             caTTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgag                                                   >  2:303392/1‑94 (MQ=255)
                                                                                                                             caTTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGAAAATGAACGCATACTGATCAAGGCGGCAGTACACGTTGTTTTTAAAGCCGgag                                                   <  1:303392/94‑1 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc    >  1:763395/1‑139 (MQ=255)
                                                                                                                                 acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  1:211251/139‑1 (MQ=255)
                                                                                                                                 acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:763395/139‑1 (MQ=255)
                                                                                                                                 acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACGGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:740474/139‑1 (MQ=255)
                                                                                                                                  cgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   >  2:3410/1‑41 (MQ=255)
                                                                                                                                  cgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   <  1:3410/41‑1 (MQ=255)
                                                                                                                                  cgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  2:731885/1‑49 (MQ=255)
                                                                                                                                  cgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           <  1:731885/49‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
ATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508315‑3508583

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