Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,359,996 T→C 7.4% F181F (TTT→TTC proCC → redundant pyrroline‑5‑carboxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,359,9960TC 7.4% 154.6 / 4.0 54F181F (TTT→TTCproCCredundant pyrroline‑5‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base T (19/31);  new base C (2/2);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

AGCACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCTGGAACGT  >  BS168_reference_genome/1359858‑1360134
                                                                                                                                          |                                                                                                                                          
aGCACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCttt                                                                                                                                            >  6:12571/1‑139 (MQ=255)
  cACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga                                                                                                                                          <  5:338876/139‑1 (MQ=255)
     tGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTgaag                                                                                                                                        >  3:202498/1‑138 (MQ=255)
      ggCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTgaagaa                                                                                                                                      <  1:158534/139‑1 (MQ=255)
       gCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCGGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAAc                                                                                                                                     <  7:273445/139‑1 (MQ=255)
       gCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTgaagaa                                                                                                                                      >  7:131883/1‑138 (MQ=255)
           cAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTTCAATCTTTGAAGAACTgg                                                                                                                                  >  3:231971/1‑138 (MQ=255)
           cAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTgg                                                                                                                                  >  6:152429/1‑138 (MQ=255)
            agagAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga                                                                                                                                          <  4:28257/129‑1 (MQ=255)
            agagAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga                                                                                                                                          >  3:28257/1‑129 (MQ=255)
            agagAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCa                                                                                                                                >  6:122235/1‑139 (MQ=255)
                aaaaaaGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTgg                                                                                                                            <  1:368422/139‑2 (MQ=255)
                 aaaaaGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTc                                                                                                                           >  8:305598/1‑139 (MQ=255)
                   acaGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTg                                                                                                                         <  1:123169/137‑1 (MQ=255)
                          gTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGt                                                                                                                     >  8:401135/1‑136 (MQ=255)
                          gTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAg                                                                                                                  <  2:52425/139‑1 (MQ=255)
                          gTCTGAATGAATTGGTATAGGTGTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAg                                                                                                                  <  2:151887/139‑1 (MQ=255)
                                    atTGTTTTCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCGAGCAATTTAACGAGCTCGGCACCGGGATTCATAGGTGCAATCTTTGAAGAACTGGCAATGTCTGCCGTCAGAAACagcagc                                                                                                        <  7:53805/139‑1 (MQ=255)
                                    atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  3:184324/139‑1 (MQ=255)
                                    atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  7:305598/139‑1 (MQ=255)
                                    atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  3:382363/139‑1 (MQ=255)
                                    atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  4:186581/139‑1 (MQ=255)
                                    atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  5:152429/139‑1 (MQ=255)
                                     ttGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCGCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  3:291938/138‑1 (MQ=255)
                                         tATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTc                                                                                                   >  5:72648/1‑139 (MQ=255)
                                         tATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTc                                                                                                   >  5:66677/1‑139 (MQ=255)
                                                  tCAGCCGTGTACGGGAAATCAAAGTATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCGATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAaagaag                                                                                           <  7:11343/138‑1 (MQ=255)
                                                  tCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCa                                                                                                 >  2:49459/1‑132 (MQ=255)
                                                  tCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCa                                                                                                 <  1:49459/132‑1 (MQ=255)
                                                        gtgtACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAAATGGCACTGTCTGCCGTCAGAAACAGCAGCCTCTCCAAACAAGAAGCCtt                                                                                    >  5:354493/1‑139 (MQ=255)
                                                         tgtACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCttt                                                                                   >  3:103215/1‑139 (MQ=255)
                                                               ggAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGAtttt                                                                             <  5:12571/139‑1 (MQ=255)
                                                                   aTCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGa                                                                         >  1:10136/1‑139 (MQ=255)
                                                                           aTCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCg                                                                 <  1:293662/139‑1 (MQ=255)
                                                                           aTCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATGCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCg                                                                 <  3:399676/139‑1 (MQ=255)
                                                                            tCGAATTTAGATGTGGCCATCAAATTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGGCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAATCCTTTGATTTTCTGATTCATTCGc                                                                <  2:54808/139‑1 (MQ=255)
                                                                                aTTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCtata                                                            <  4:103215/139‑1 (MQ=255)
                                                                                     gATGAGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGaa                                                       <  2:10136/139‑1 (MQ=255)
                                                                                           gCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga                                                                                                                                          <  6:85924/50‑1 (MQ=255)
                                                                                           gCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga                                                                                                                                          >  5:85924/1‑50 (MQ=255)
                                                                                                      caCGAGCTCGGCACCGGTATTCATAGAGGCCATCTTTTAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGAt                                       <  4:231971/137‑1 (MQ=255)
                                                                                                         gAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGaa                                   <  5:122235/139‑1 (MQ=255)
                                                                                                           gCTCGGCACCGGGATTCATAGCTGCAATCTATGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGaaa                                  >  7:84375/1‑138 (MQ=255)
                                                                                                             tCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        >  1:286317/1‑66 (MQ=255)
                                                                                                             tCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc                                                                                                        <  2:286317/66‑1 (MQ=255)
                                                                                                                 cACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAatat                           <  6:398901/139‑1 (MQ=255)
                                                                                                                          tCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTc                    >  1:296969/1‑137 (MQ=255)
                                                                                                                            aTAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGTTTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgagg                <  6:66677/139‑1 (MQ=255)
                                                                                                                            aTAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGTTTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgagg                <  6:72648/139‑1 (MQ=255)
                                                                                                                             tAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgagga               <  7:401135/139‑1 (MQ=255)
                                                                                                                               gCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgaggag              <  1:28569/138‑1 (MQ=255)
                                                                                                                                   cAATCTTTGAAGCACTGTCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTTATTTTCTGATTCATTCGCTATACGTAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCt         >  7:38860/1‑139 (MQ=255)
                                                                                                                                   cAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGattcatt                                                                   >  5:387890/1‑81 (MQ=255)
                                                                                                                                   cAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGattcatt                                                                   <  6:387890/81‑1 (MQ=255)
                                                                                                                                   cAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCt         >  7:293278/1‑139 (MQ=255)
                                                                                                                                          tGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCTGGAACgt  >  1:283010/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                          
AGCACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCTGGAACGT  >  BS168_reference_genome/1359858‑1360134

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: