Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,309,246 T→C I60V (ATT→GTT)  phnJ ← alpha‑D‑ribose 1‑methylphosphonate 5‑phosphate C‑P‑lyase PhnJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,309,2460TC78.8% 101.2 / 27.7 52I60V (ATT→GTT) phnJalpha‑D‑ribose 1‑methylphosphonate 5‑phosphate C‑P‑lyase PhnJ
Reads supporting (aligned to +/- strand):  ref base T (7/4);  new base C (24/17);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-04

ATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCG  >  NZ_CP009273/4309113‑4309372
                                                                                                                                     |                                                                                                                              
aTAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGc                                                                                                                           <  2:757593/139‑1 (MQ=255)
     cGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCg                                                                                                                      <  2:389817/139‑1 (MQ=255)
       tCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCgg                                                                                                                     >  2:741365/1‑138 (MQ=255)
         gCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGt                                                                                                                    >  1:89938/1‑137 (MQ=255)
         gCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGt                                                                                                                    >  2:722212/1‑137 (MQ=255)
                         ccGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCgg                                                                                                   <  1:57138/138‑1 (MQ=255)
                         ccGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGt                                                                                                  >  1:124570/1‑139 (MQ=255)
                         ccGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGt                                                                                                  >  1:682144/1‑139 (MQ=255)
                         ccGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGt                                                                                                  >  2:421573/1‑139 (MQ=255)
                                 tttACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGc                                                                                                                           >  1:131694/1‑106 (MQ=255)
                                 tttACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGc                                                                                                                           <  2:131694/106‑1 (MQ=255)
                                         ggTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTAcgg                                                                                   >  2:644614/1‑138 (MQ=255)
                                                    tGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCAtctc                                                                       >  1:342540/1‑139 (MQ=255)
                                                           gTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTa                                                                                      <  1:514460/117‑1 (MQ=255)
                                                           gTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTa                                                                                      >  2:514460/1‑117 (MQ=255)
                                                           gTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGcgcccg                                                                >  1:686516/1‑139 (MQ=255)
                                                                   tCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                        <  2:724714/139‑1 (MQ=255)
                                                                     gACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCa                                                      <  2:686516/139‑1 (MQ=255)
                                                                              ttGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAAc                                             <  2:682144/139‑1 (MQ=255)
                                                                                      gtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGc                                                                                                                           >  2:697432/1‑53 (MQ=255)
                                                                                      gtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGc                                                                                                                           <  1:697432/53‑1 (MQ=255)
                                                                                      gtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTa                                                                                      <  1:608027/90‑1 (MQ=255)
                                                                                      gtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTa                                                                                      >  2:608027/1‑90 (MQ=255)
                                                                                      gtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATc                                     >  2:325475/1‑139 (MQ=255)
                                                                                        cATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCg                                           <  1:660570/131‑1 (MQ=255)
                                                                                        cATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCg                                           >  2:660570/1‑131 (MQ=255)
                                                                                         aTCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTc                                                                                                >  1:84429/1‑77 (MQ=255)
                                                                                         aTCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTc                                                                                                <  2:84429/77‑1 (MQ=255)
                                                                                          tCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCa                                 <  1:576589/139‑1 (MQ=255)
                                                                                            cgcgCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAcc                               <  2:390515/139‑1 (MQ=255)
                                                                                                          tCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGc                                                       >  1:405740/1‑101 (MQ=255)
                                                                                                          tCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGc                                                       <  2:405740/101‑1 (MQ=255)
                                                                                                          tCACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACC‑‑‑GATCGCCACCGCTTTTAAGATGGcgcg              >  2:684290/1‑139 (MQ=255)
                                                                                                           cACCTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGgcg                >  1:256792/1‑139 (MQ=255)
                                                                                                             ccTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTc                                                                                                >  1:106134/1‑57 (MQ=255)
                                                                                                             ccTTCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTc                                                                                                <  2:106134/57‑1 (MQ=255)
                                                                                                               ttCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCgcggc            >  2:864120/1‑139 (MQ=255)
                                                                                                               ttCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                        >  2:336360/1‑95 (MQ=255)
                                                                                                               ttCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                        <  1:336360/95‑1 (MQ=255)
                                                                                                               ttCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCgcggc            >  2:289539/1‑139 (MQ=255)
                                                                                                               ttCAGCACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGCCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCgcggc            >  1:793212/1‑139 (MQ=255)
                                                                                                                    cACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                <  2:193143/114‑1 (MQ=255)
                                                                                                                    cACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                >  1:193143/1‑114 (MQ=255)
                                                                                                                    cACGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGAt        >  2:706631/1‑138 (MQ=255)
                                                                                                                     aCGTCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCa      <  2:89938/139‑1 (MQ=255)
                                                                                                                        tCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGgc                 >  2:330163/1‑125 (MQ=255)
                                                                                                                        tCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGgc                 <  1:330163/125‑1 (MQ=255)
                                                                                                                        tCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGc   >  1:775210/1‑139 (MQ=255)
                                                                                                                        tCGCTTTCGCCAACCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGc   >  1:32542/1‑139 (MQ=255)
                                                                                                                         cGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                >  2:650536/1‑109 (MQ=255)
                                                                                                                         cGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                <  1:650536/109‑1 (MQ=255)
                                                                                                                         cGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCg  >  1:237151/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                              
ATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCG  >  NZ_CP009273/4309113‑4309372

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