Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,937 C→A 12.0% Q78H (CAG→CAT)  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,9370CA12.0% 143.7 / 7.4 50Q78H (CAG→CAT) fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base C (17/27);  new base A (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 = 5.80e-01

TGCCTGCAATTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGG  >  NC_002947/243790‑244081
                                                                                                                                                   |                                                                                                                                                 
tgcctgcAATTCGCTGTACAGCACACCTTCCGTTAGTTCTTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCActgcctg                                                                                                                                                 >  1:100446/1‑150 (MQ=255)
 gcctgcAATTCGCTGTACAGCACACCTNCCGTTAGTTCGTNCACCCCNTNACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCActgcctgc                                                                                                                                                >  1:696950/1‑150 (MQ=255)
       aaTTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTCGCCATTAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCGGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑Gt                                                                                                                                         >  2:174802/1‑150 (MQ=255)
            gCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               <  1:1069130/140‑1 (MQ=255)
            gCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               >  2:1069130/1‑140 (MQ=255)
            gCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCtt                                                                                                                                    <  1:1211800/150‑1 (MQ=255)
             cTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCGGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTTGCTGCATGC‑GTGCCTTg                                                                                                                                   <  2:315443/150‑1 (MQ=255)
              tGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGt                                                                                                                                  <  2:603597/150‑1 (MQ=255)
               gTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTc                                                                                                                                 <  1:985558/150‑1 (MQ=255)
               gTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTc                                                                                                                                 <  1:429371/150‑1 (MQ=255)
                tACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCg                                                                                                                                >  2:1090824/1‑150 (MQ=255)
                tACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCg                                                                                                                                >  1:723286/1‑150 (MQ=255)
                 aCAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGa                                                                                                                               <  2:723286/150‑1 (MQ=255)
                         ccTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAgg                                                                                                                       <  2:777850/150‑1 (MQ=255)
                          cTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAg                                                                                                                        >  1:1106288/1‑148 (MQ=255)
                          cTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAg                                                                                                                        <  2:1106288/148‑1 (MQ=255)
                             ccGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGtt                                                                                                                   <  2:310674/150‑1 (MQ=255)
                              cGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCCGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTg                                                                                                                  >  2:392781/1‑150 (MQ=255)
                               gTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGc                                                                                                                 >  1:285756/1‑150 (MQ=255)
                                    ttCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               >  1:747044/1‑116 (MQ=255)
                                    ttCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCa                                                                                                                                               <  2:747044/116‑1 (MQ=255)
                                     tCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGa                                                                                                           <  2:285756/150‑1 (MQ=255)
                                        tGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCTTTGTCGAGTTTCAGGGGTTGCCCCTGATAg                                                                                                        <  2:976002/150‑1 (MQ=255)
                                        tGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAg                                                                                                        <  2:910696/150‑1 (MQ=255)
                                           aCCCCCTACCGCCCCCGTTCGTTCGCCCTGACACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGCTAGGAc                                                                                                     <  2:10322/150‑1 (MQ=255)
                                              ccATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTc                                                                                                                                 >  1:1202430/1‑119 (MQ=255)
                                              ccATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTc                                                                                                                                 <  2:1202430/119‑1 (MQ=255)
                                                  aaCGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑Gt                                                                                                                                         <  2:426087/107‑1 (MQ=255)
                                                  aaCGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑Gt                                                                                                                                         >  1:426087/1‑107 (MQ=255)
                                                         aCGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGa                                                                                                      >  1:145768/1‑135 (MQ=255)
                                                         aCGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGa                                                                                                      <  2:145768/135‑1 (MQ=255)
                                                           gTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCTCCTGATAGGa                                                                                                      >  2:762560/1‑133 (MQ=255)
                                                           gTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCTCCTGATAGGa                                                                                                      <  1:762560/133‑1 (MQ=255)
                                                              cATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCa                                                                                  <  1:1090824/150‑1 (MQ=255)
                                                                   gCCATGAAACGCGTTTCGGCGGCCTGCAGCGTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCa                                                                             >  2:1154941/1‑150 (MQ=255)
                                                                    ccATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCAt                                                                            >  2:261770/1‑150 (MQ=255)
                                                                      aTGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATcc                                                                          >  2:1123369/1‑150 (MQ=255)
                                                                                tttCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGt                                                                                                                    <  2:872515/98‑1 (MQ=255)
                                                                                tttCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGt                                                                                                                    >  1:872515/1‑98 (MQ=255)
                                                                                   cggcggCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTg                                                             <  2:696950/150‑1 (MQ=255)
                                                                                     gcggcCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGc                                                           <  2:42997/150‑1 (MQ=255)
                                                                                       ggcCTGCAGCTTTTGCGCGGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCAt                                                         <  1:174802/150‑1 (MQ=255)
                                                                                            gCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCt                                                    <  1:1123369/150‑1 (MQ=255)
                                                                                                  tttGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCa                                              <  1:392781/150‑1 (MQ=255)
                                                                                                              cgTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCa                                  <  2:100446/150‑1 (MQ=255)
                                                                                                                                      gacaagGCAC‑GCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGAc                                                                                    <  2:502773‑M1/70‑1 (MQ=255)
                                                                                                                                      gacaagGCAC‑GCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGAc                                                                                    >  1:502773‑M1/7‑76 (MQ=255)
                                                                                                                                        caagGCAC‑GCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       >  1:808213‑M1/5‑87 (MQ=255)
                                                                                                                                        caagGCAC‑GCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       >  1:808161‑M1/5‑87 (MQ=255)
                                                                                                                                        caagGCAC‑GCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       <  2:808161‑M1/83‑1 (MQ=255)
                                                                                                                                        caagGCAC‑GCATGCAGGCAGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGc                                                                       <  2:808213‑M1/83‑1 (MQ=255)
                                                                                                                                              ACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTgg  >  1:261770‑M1/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCCTGCAATTCGCTGTACAGCACACCTTCCGTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGC‑GTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGG  >  NC_002947/243790‑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: ‑