Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,924,567 A→T 8.6% F115Y (TTT→TAT)  PP_1724 ← 2‑aminoethylphosphonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,924,5670AT 8.6% 147.2 / 3.8 48F115Y (TTT→TAT) PP_17242‑aminoethylphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base T (2/2);  total (24/24)
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.84e-01

GCAGCAAGACGCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCTCGACCAGAAGCTGATTTCCAG  >  NC_002947/1924421‑1924715
                                                                                                                                                  |                                                                                                                                                    
gcagcaAGACGCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGAAATGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAgg                                                                                                                                                   >  2:1261607/1‑150 (MQ=255)
   gcaAGACGCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCa                                                                                                                                                >  1:823836/1‑150 (MQ=255)
      aGACGCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCg                                                                                                                                             <  2:1208457/150‑1 (MQ=255)
      aGACGCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCg                                                                                                                                             >  1:721267/1‑150 (MQ=255)
          gCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTAc                                                                                                                                         <  1:675945/150‑1 (MQ=255)
                    aTCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAg                                                                                                                               <  2:533462/150‑1 (MQ=255)
                        ggAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGc                                                                                                                           <  2:309576/150‑1 (MQ=255)
                                gtgtACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATg                                                                                                                   >  2:450936/1‑150 (MQ=255)
                                 tgtACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGAAGGATAACAAAGCCAGAGCCGAGCCTTACGCCGGCAAATTTGCTGGTCAtat                                                                                                                  >  2:727048/1‑148 (MQ=255)
                                    acacCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTgg                                                                                                               <  1:450936/150‑1 (MQ=255)
                                     cacCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGc                                                                                                              <  2:823836/150‑1 (MQ=255)
                                      acCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCg                                                                                                             >  1:887242/1‑150 (MQ=255)
                                        cAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTAAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGaa                                                                                                           >  1:457841/1‑150 (MQ=255)
                                         aGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAg                                                                                                          >  1:1028683/1‑150 (MQ=255)
                                           gATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTc                                                                                                                      <  2:1127089/136‑1 (MQ=255)
                                           gATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTc                                                                                                                      >  1:1127089/1‑136 (MQ=255)
                                           gATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGgc                                                                                                        >  2:1004141/1‑150 (MQ=255)
                                             tCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCTAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCATAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTgg                                                                                                               <  1:684862/141‑1 (MQ=255)
                                             tCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCTAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCATAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTgg                                                                                                               >  2:684862/1‑141 (MQ=255)
                                             tCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCATAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTgg                                                                                                               <  1:684934/141‑1 (MQ=255)
                                             tCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCATAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTgg                                                                                                               >  2:684934/1‑141 (MQ=255)
                                               aaGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTg                                                                                                    >  1:1127715/1‑150 (MQ=255)
                                                 gCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGa                                                                                                   >  2:1367023/1‑149 (MQ=255)
                                                 gCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGa                                                                                                   <  1:1367023/149‑1 (MQ=255)
                                                 gCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGAc                                                                                                  >  1:1262840/1‑150 (MQ=255)
                                                  ccGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACg                                                                                                 >  1:629764/1‑150 (MQ=255)
                                                            aaTAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCa                                                                                       <  2:1421518/150‑1 (MQ=255)
                                                                gATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCtt                                                                                   >  1:1281595/1‑150 (MQ=255)
                                                                    cTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAg                                                                               <  1:1004141/150‑1 (MQ=255)
                                                                     tGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCTAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGACATAAAGACTATGGCGCTGACGAAGTCGCGCAGCTTCGAGt                                                                              >  2:474153/1‑150 (MQ=255)
                                                                      gAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGTGCTCACGAAGTCGCGCAGCTTCGAGTc                                                                             <  1:1261607/150‑1 (MQ=255)
                                                                        aGTCCTCCAGCAGGCCGACCTGCTTCGGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAACTCGCGCAGCTTCCAGTCCa                                                                           >  2:1087567/1‑150 (MQ=255)
                                                                         gTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCAc                                                                          <  2:1281595/150‑1 (MQ=255)
                                                                              ccAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGgc                                                                     <  2:1028683/150‑1 (MQ=255)
                                                                               cagcagGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGgcg                                                                    >  2:546647/1‑150 (MQ=255)
                                                                                agcagGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGgcgc                                                                   <  2:721267/150‑1 (MQ=255)
                                                                                   agGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTCCGTCTAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGt                                                                <  1:727048/150‑1 (MQ=255)
                                                                                   agGCCGACCTGCTTAAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGt                                                                <  2:457841/150‑1 (MQ=255)
                                                                                      ccGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGaa                                                             <  2:1262840/150‑1 (MQ=255)
                                                                                       cGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGGTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGATACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAAt                                                            >  2:100767/1‑150 (MQ=255)
                                                                                                           gTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAgg                                                                                                         >  2:1393654/1‑85 (MQ=255)
                                                                                                           gTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAgg                                                                                                         <  1:1393654/85‑1 (MQ=255)
                                                                                                                            gCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCt                       >  1:946684/1‑150 (MQ=255)
                                                                                                                             cctaGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAAACACCGTCGAGGGTAGCGATGACGATGCCGAACACGCTc                      <  1:474153/147‑1 (MQ=255)
                                                                                                                             ccAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCTc                      <  2:946684/150‑1 (MQ=255)
                                                                                                                              cAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCTCg                     <  2:887242/150‑1 (MQ=255)
                                                                                                                                    ggATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATCTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCTCGACCAGa               <  1:1249205/150‑1 (MQ=255)
                                                                                                                                               gCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCTCGACCAGAAGCTGATTTcc    >  2:503806/1‑150 (MQ=255)
                                                                                                                                                 aaaGGCCAGCGGTACGCCGGCAAAGTGGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGGGCAGCTTCGAGTACACCGGTCGCAGTGGATAGGCGCCGAGGGTAGCGATTACGATGCCGAACAGGCTCGACCAGAGGCTGATTTcccg  >  1:174087/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCAGCAAGACGCCCAGCGGGATCTGGAAGTAGGTGTACACCAGGATCAAGCCGCTTTTGGAATAGATACTGAAGTCCTCCAGCAGGCCGACCTGCTTCAGCAGCAGCGTCAGCGCGCCGTTGAAGCCAAGCAGGATGATGAAGGCAAAGGCCAGCGGTACGCCGGCAAAGTTGCTGGTCATGTTGGCGAAGGCGCTGACGAAGTCGCGCAGCTTCGAGTCCACCTGGCGCAGTGAATAGGCGCCGAGGGTAGCGATGACGATGCCGAACACGCTCGACCAGAAGCTGATTTCCAG  >  NC_002947/1924421‑1924715

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑