Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,934 Δ1 bp 12.0% coding (237/663 nt) fkl ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947243,9340T.12.0% 143.2 / 14.8 50coding (237/663 nt)fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (17/27);  new base . (3/3);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GCCTGCCTGCAATTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGG  >  NC_002947/243787‑244081
                                                                                                                                                   |                                                                                                                                                    
gcctgcctGCAATTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCActgc                                                                                                                                                    <  1:568544/150‑1 (MQ=255)
   tgcctgcAATTCGCTGTACAGCACACCTTCCGTTAGTTCTTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCActgcctg                                                                                                                                                 >  1:100446/1‑150 (MQ=255)
    gcctgcAATTCGCTGTACAGCACACCTNCCGTTAGTTCGTNCACCCCNTNACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCActgcctgc                                                                                                                                                >  1:696950/1‑150 (MQ=255)
          aaTTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTCGCCATTAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCGGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGt                                                                                                                                         >  2:174802/1‑150 (MQ=255)
               gCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               <  1:1069130/140‑1 (MQ=255)
               gCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               >  2:1069130/1‑140 (MQ=255)
               gCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCtt                                                                                                                                    <  1:1211800/150‑1 (MQ=255)
                cTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCGGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTTGCTGCATGCGTGCCTTg                                                                                                                                   <  2:315443/150‑1 (MQ=255)
                 tGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGt                                                                                                                                  <  2:603597/150‑1 (MQ=255)
                  gTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTc                                                                                                                                 <  1:429371/150‑1 (MQ=255)
                  gTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTc                                                                                                                                 <  1:985558/150‑1 (MQ=255)
                   tACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCg                                                                                                                                >  1:723286/1‑150 (MQ=255)
                   tACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCg                                                                                                                                >  2:1090824/1‑150 (MQ=255)
                    aCAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGa                                                                                                                               <  2:723286/150‑1 (MQ=255)
                            ccTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAgg                                                                                                                       <  2:777850/150‑1 (MQ=255)
                             cTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAg                                                                                                                        >  1:1106288/1‑148 (MQ=255)
                             cTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAg                                                                                                                        <  2:1106288/148‑1 (MQ=255)
                                ccGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGtt                                                                                                                   <  2:310674/150‑1 (MQ=255)
                                 cGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCCGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTg                                                                                                                  >  2:392781/1‑150 (MQ=255)
                                  gTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGc                                                                                                                 >  1:285756/1‑150 (MQ=255)
                                       ttCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               >  1:747044/1‑116 (MQ=255)
                                       ttCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               <  2:747044/116‑1 (MQ=255)
                                        tCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGa                                                                                                           <  2:285756/150‑1 (MQ=255)
                                           tGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCTTTGTCGAGTTTCAGGGGTTGCCCCTGATAg                                                                                                        <  2:976002/150‑1 (MQ=255)
                                           tGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAg                                                                                                        <  2:910696/150‑1 (MQ=255)
                                              aCCCCCTACCGCCCCCGTTCGTTCGCCCTGACACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGCTAGGAc                                                                                                     <  2:10322/150‑1 (MQ=255)
                                                 ccATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTc                                                                                                                                 <  2:1202430/119‑1 (MQ=255)
                                                 ccATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTc                                                                                                                                 >  1:1202430/1‑119 (MQ=255)
                                                     aaCGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGt                                                                                                                                         <  2:426087/107‑1 (MQ=255)
                                                     aaCGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGt                                                                                                                                         >  1:426087/1‑107 (MQ=255)
                                                            aCGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGa                                                                                                      >  1:145768/1‑135 (MQ=255)
                                                            aCGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGa                                                                                                      <  2:145768/135‑1 (MQ=255)
                                                              gTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCTCCTGATAGGa                                                                                                      >  2:762560/1‑133 (MQ=255)
                                                              gTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCTCCTGATAGGa                                                                                                      <  1:762560/133‑1 (MQ=255)
                                                                 cATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCa                                                                                  <  1:1090824/150‑1 (MQ=255)
                                                                      gCCATGAAACGCGTTTCGGCGGCCTGCAGCGTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCa                                                                             >  2:1154941/1‑150 (MQ=255)
                                                                       ccATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCAt                                                                            >  2:261770/1‑150 (MQ=255)
                                                                         aTGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATcc                                                                          >  2:1123369/1‑150 (MQ=255)
                                                                                   tttCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGt                                                                                                                    >  1:872515/1‑98 (MQ=255)
                                                                                   tttCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGt                                                                                                                    <  2:872515/98‑1 (MQ=255)
                                                                                      cggcggCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTg                                                             <  2:696950/150‑1 (MQ=255)
                                                                                        gcggcCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGc                                                           <  2:42997/150‑1 (MQ=255)
                                                                                          ggcCTGCAGCTTTTGCGCGGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCAt                                                         <  1:174802/150‑1 (MQ=255)
                                                                                               gCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCt                                                    <  1:1123369/150‑1 (MQ=255)
                                                                                                     tttGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCa                                              <  1:392781/150‑1 (MQ=255)
                                                                                                                 cgTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCa                                  <  2:100446/150‑1 (MQ=255)
                                                                                                                                         gacaagGCACGCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGAc                                                                                    <  2:502773‑M1/70‑1 (MQ=255)
                                                                                                                                         gacaagGCACGCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGAc                                                                                    >  1:502773‑M1/7‑76 (MQ=255)
                                                                                                                                           caagGCACGCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       >  1:808161‑M1/5‑87 (MQ=255)
                                                                                                                                           caagGCACGCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       <  2:808161‑M1/83‑1 (MQ=255)
                                                                                                                                           caagGCACGCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       <  2:808213‑M1/83‑1 (MQ=255)
                                                                                                                                           caagGCACGCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       >  1:808213‑M1/5‑87 (MQ=255)
                                                                                                                                                 ACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTgg  >  1:261770‑M1/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCTGCCTGCAATTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGG  >  NC_002947/243787‑244081

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