Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 648,025 T→G 14.5% F29C (TTT→TGT)  tatCY → component of the twin‑arginine pre‑protein translocation pathway

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome648,0250TG14.5% 125.6 / 6.9 58F29C (TTT→TGT) tatCYcomponent of the twin‑arginine pre‑protein translocation pathway
Reads supporting (aligned to +/- strand):  ref base T (16/32);  new base G (7/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 2.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACAGACCCGCTTTATGTGTTT  >  BS168_reference_genome/647889‑648149
                                                                                                                                        |                                                                                                                             
ggATTAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTg                                                                                                                             <  1:238318/139‑1 (MQ=255)
 gATTAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGt                                                                                                                            >  8:292969/1‑139 (MQ=255)
 gATTAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTATga                                                                                                                            >  7:53183/1‑138 (MQ=255)
   ttAACTAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGtag                                                                                                                          >  2:181380/1‑137 (MQ=255)
   ttAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACATGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCg                                                                                                                          >  8:303304/1‑139 (MQ=255)
   ttAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGtgtgggg                                                                                                                          >  8:103578/1‑136 (MQ=255)
     aaCGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGtt                                                                                                                        >  3:393271/1‑139 (MQ=255)
     aaCGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTGTGTCGtt                                                                                                                        >  4:20490/1‑139 (MQ=255)
     aaCGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTGTGGCGtg                                                                                                                        >  5:303745/1‑138 (MQ=255)
      aCGAGTGCTGAAGAGGAAACACAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGtt                                                                                                                        <  6:303745/138‑1 (MQ=255)
      aCGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGtt                                                                                                                        <  1:232517/138‑1 (MQ=255)
                    ggAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGa                                                                                                         <  8:413536/139‑1 (MQ=255)
                       aaaaaaGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGAttt                                                                                                     >  8:240953/1‑139 (MQ=255)
                       aaaaaaGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTGTGTGGTTGTTTGTATTTGTGGCttt                                                                                                     >  1:135227/1‑139 (MQ=255)
                                agaagcTCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCa                                                                                            <  4:398310/139‑1 (MQ=255)
                                agaagaTCAGTAAGTTAGGATGACATGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCa                                                                                            <  7:303304/139‑1 (MQ=255)
                                agaagaTCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCa                                                                                            <  7:204361/139‑1 (MQ=255)
                                  aagaTCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCaaa                                                                                          <  1:315816/139‑1 (MQ=255)
                                   agaTCAGTTAGTTAGGATGAGACGAATGTAAGTGTATCTAATGTCGGTGCTGGAGGATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTATCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCaaa                                                                                          <  3:398995/138‑1 (MQ=255)
                                        aGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGat                                                                                    <  4:36452/139‑1 (MQ=255)
                                               tAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTTTCGTGTTCTTTATTGCTGGATTTTTTTTTGCAAAGCCGATTATtgtg                                                                             >  1:38471/1‑139 (MQ=255)
                                               tAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATtgtg                                                                             >  6:310689/1‑139 (MQ=255)
                                                 ggATGACACGAATGACAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTgtgt                                                                            <  2:328724/138‑1 (MQ=255)
                                                 ggATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTgtgt                                                                            <  8:169857/138‑1 (MQ=255)
                                                 ggATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTGTGTCGTTTTGTTTTTTTGTGGTTTTTTTTTATGAAAAACCATTTTTgtgtt                                                                           >  3:374984/1‑138 (MQ=255)
                                                  gATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTAt                                                                          <  3:20490/139‑1 (MQ=255)
                                                      acacGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTAGCAAAGCCGATTATTGTGTATCTGCa                                                                     <  2:71714/139‑1 (MQ=255)
                                                          gaatgaatGTGAATCAAATGTCGCTGATGGAGCATATTGCTGAGCTTCGAATACGGTTGCTGATTTTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGa                                                                  <  8:194656/139‑1 (MQ=255)
                                                                 aGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTTCTGGTTTTTTTTTTGCAAAGCCGATTATTTTGTATCAGCAAGAAACAGAt                                                           >  8:2589/1‑139 (MQ=255)
                                                                         aaaTGTCGCTGCTGGAGCTTATTGCTGAGCTTCGAAAACGGTGGCTGAGTGTAGCGCTGGCGTTTGTCGTATTATTTATTGCTGGATTTTTTTTAGCATAGCCGATTATTGTGTATCTGCAAGAAACAGATgaagcgaa                                                   <  8:235640/139‑1 (MQ=255)
                                                                         aaaTGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCGGATTATTGTGTATCTGCAAGAAACAGATgaagcgaa                                                   <  7:302266/139‑1 (MQ=255)
                                                                         aaaTGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACAGTTGCTGATTGTAGCGCTGGCGTGTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATgaagcgaa                                                   <  3:146676/139‑1 (MQ=255)
                                                                          aaTGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATgaagcgaag                                                  <  3:347901/139‑1 (MQ=255)
                                                                          aaTGTCGCGGCTGGAGCATATTGCTGAGCTTCGAAAACGGGTGCTGATTGTAGCGCTGGCGTTTGTCGGTTTCTTTATTTCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGTTgaagcgaag                                                  <  7:359350/139‑1 (MQ=255)
                                                                          aaTGCCGCTGCTGGAGCATAATGCTGAGCTTCGAAAACGGTTGCTGCATGTAGCGCCTGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTTAGCAAAGCCGATAATTGTGTATCTGCAAGAAACAGATgaagcgaa                                                   <  4:391808/139‑1 (MQ=255)
                                                                               cGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGTCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCt                                             <  7:412713/139‑1 (MQ=255)
                                                                               cGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTGTGTCGGGTTGTTTTTTGGTGGAGTTTTTTTTGAAAAGCCGATTTTTTTGTTTATGGAAGAAACAGATGAAGCGAagcagcg                                             >  8:173882/1‑138 (MQ=255)
                                                                                gctgctGGAGCATATTGCTGTGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGTAACAGATGAAGCGAAGCAGCt                                             <  5:253920/138‑1 (MQ=255)
                                                                                 ctgctgGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTa                                           <  4:393271/139‑1 (MQ=255)
                                                                                 ctgctgGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTGTGTCGTTTTGTTTTTTGGTGGTTTTTTTTTTGCAAATCCGATTATTTGTTTTAGGCAAGAAAAACATTAAGCCGAGCAGCtt                                            >  4:297629/1‑138 (MQ=255)
                                                                                        aGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTaa                                    <  2:127892/139‑1 (MQ=255)
                                                                                        aGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTaa                                    <  7:240953/139‑1 (MQ=255)
                                                                                        aGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGAGTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTaa                                    <  7:2589/139‑1 (MQ=255)
                                                                                            tatTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTGGGTTTGTTTTTTTGTGGATTTTTTTTTTCAAAGCCGATTTTTTTGTTTATGCAAGAAACAGATTAATCGAAGCCGCTTAAGCGTAAagcg                                >  6:368084/1‑139 (MQ=255)
                                                                                            tatTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTGTTTTTTTGGGGTTTTTTTTTTGGAAAGCCGAATTTTGTGTATATGCAAGAAACAGATGAAGCGCATCAGCATACGCTTTAcgcg                                >  3:291707/1‑139 (MQ=255)
                                                                                              atGCTGAGCTTCGAAAACGGGTGCTGACTGTAGAGCTGGCGTGTGTCGTTTTCTTAATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGtt                              <  4:194696/138‑1 (MQ=255)
                                                                                                gCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGTCGTTTGTCGTTTTCTTTATTTCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGttt                             <  8:53183/138‑1 (MQ=255)
                                                                                                 cTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCGGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTaa                           <  2:135227/139‑1 (MQ=255)
                                                                                                    aGCTTCGAAAACGGTTGCTGATTGTGGCGCTGGCGTTTGTCGTTTTAATTATTTCAGGATTTTTTTTTCCAAACCCCATTATTTTTTATAAGCAAAAAACAAATAAACCGAACCACCTTACCCTTAACCCGTTTAAca                         >  6:367671/1‑137 (MQ=255)
                                                                                                      cTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGa                      <  5:126487/139‑1 (MQ=255)
                                                                                                         cGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTGATTTTTTCTTGATTTTTTTTAGCAAAGCCCATTTTTTTTTATGTGCAAGAAACAAATGAAGCGAAGCACCATACCCTTAAACCGTTTAAACtta                      >  6:126487/1‑134 (MQ=255)
                                                                                                         cGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACa                    <  2:22603/138‑1 (MQ=255)
                                                                                                         cGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTATTTATTTCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGa                      >  2:13057/1‑136 (MQ=255)
                                                                                                          gAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACAGa                  <  4:291707/139‑1 (MQ=255)
                                                                                                           aaaaCGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGGTGGATTTTTTTTTGCAAACCCGATTATTGTGTTTCTTCAAGAAACAGATGTAGCCAAGCAGCTAACGCTTAACGCCGTTAACCTGACAGAc                 >  7:350754/1‑139 (MQ=255)
                                                                                                            aaaCGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACAGAcc                <  1:167473/139‑1 (MQ=255)
                                                                                                                ggTTGCTGATTGTAGCGCTGGCGTGTGTCGTTTTCTTTTTTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACAGACCCGCt            >  5:159352/1‑139 (MQ=255)
                                                                                                                          tGTAGCGCTGGCGTTTGTCGTTTTGTTTATTTGTGGATTTTTTTTAGCAAAGCCGATTATTTTGTATCTGCAAGAAACAGATTAAGCGAAGCCGCTTACCCTTAACGCGTTTAACCTGACAGCCCCGCGTTATGTGttt  >  3:238875/1‑139 (MQ=255)
                                                                                                                          tGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACAGACCCGCTTTATGTGttt  >  7:146484/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                             
GGATTAACGAGTGATGAAGAGGAAAAAAAGAAAGAAGATCAGTAAGTTAGGATGACACGAATGAAAGTGAATCAAATGTCGCTGCTGGAGCATATTGCTGAGCTTCGAAAACGGTTGCTGATTGTAGCGCTGGCGTTTGTCGTTTTCTTTATTGCTGGATTTTTTTTAGCAAAGCCGATTATTGTGTATCTGCAAGAAACAGATGAAGCGAAGCAGCTTACGCTTAACGCGTTTAACCTGACAGACCCGCTTTATGTGTTT  >  BS168_reference_genome/647889‑648149

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