Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,291,923 T→G 40.4% F194V (TTT→GTT)  cypE → cytochrome P450 CYP109B1, monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,291,9230TG40.4% 20.9 / 19.5 44F194V (TTT→GTT) cypEcytochrome P450 CYP109B1, monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base G (16/2);  total (29/15)
Fisher's exact test for biased strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGAGTGCCTTCAGCGCATATGGAACAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTCCGGTGAAGAGCTGATTCCGTTTTGCACGCTGCT  >  BS168_reference_genome/1291785‑1292050
                                                                                                                                          |                                                                                                                               
ctgGGAGTGCCTTCAGCGCATATGGAACAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGt                                                                                                                                 >  6:392327/1‑139 (MQ=255)
ctgGGAGTGCCTTCAGCGCATATGGAACAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGACTTTTTGTAAGAACGAGATAAGTGTGATGAAGAACTGGCCGCGt                                                                                                                                 >  8:133345/1‑139 (MQ=255)
             cAGCTCATATGGAACAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGca                                                                                                                     <  1:185869/138‑1 (MQ=255)
             cAGCGCATATGGAACAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGc                                                                                                                      <  8:238788/137‑1 (MQ=255)
                         aaCAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTTCCGGCCTCATAGAAGaaaa                                                                                                        >  4:151828/1‑139 (MQ=255)
                          aCAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTTGCGGCCTCAAAgaagaa                                                                                                          >  4:9558/1‑136 (MQ=255)
                              tttAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGGCGGGCTCATTGAAGTAAAGCGa                                                                                                    >  5:390489/1‑138 (MQ=255)
                              tttAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATAAGAGAAGAAAAGCGa                                                                                                    >  7:81220/1‑138 (MQ=255)
                              tttAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGGCGGCAGCATAGAAGAAAAGCGa                                                                                                    >  4:80504/1‑138 (MQ=255)
                                tAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGaaa                                                                                                  >  5:299455/1‑138 (MQ=255)
                                   aGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGCCGGCATCATAGAAGAAAAGCGcaacaaa                                                                                              >  7:373955/1‑139 (MQ=255)
                                    gCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGaaacaaac                                                                                             <  1:74929/139‑1 (MQ=255)
                                    gCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGCCGGCATCATAGAAGAAAAGCGaaacaaac                                                                                             <  6:299455/139‑1 (MQ=255)
                                        ggTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTTGCGGGATCAAAGAAGAAAAAAGAAAAAAACCgg                                                                                          >  1:150929/1‑138 (MQ=255)
                                        ggTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGGCGGCCACATAGAAGAAAAGCGAAACAAACCgg                                                                                          >  6:149165/1‑138 (MQ=255)
                                           cTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACa                                                                                      <  6:24774/139‑1 (MQ=255)
                                              aTCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTTCCGGCATCATATAAGAAAAACGAAACAAACCCGCACAggc                                                                                   >  2:74929/1‑138 (MQ=255)
                                              aTCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTGTTTTTGCGGCCTCAAAAAATAAAAACGACACAAACCCGGACAg                                                                                     >  5:226947/1‑137 (MQ=255)
                                              aTCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTTCCGGCATCCTTGAAGAAAAGCGCAACAAACCGGAACAGGa                                                                                   >  3:87434/1‑139 (MQ=255)
                                                   ctGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGata                                                                                 <  6:326650/136‑1 (MQ=255)
                                                    tGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGGTTTTGGCGGCATGCTAGTAGAAAAGCGAAACAAACCGGAACAGGAtatttt                                                                             >  4:374907/1‑137 (MQ=255)
                                                      gTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGAtattat                                                                             <  5:149165/137‑1 (MQ=255)
                                                          gTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTGTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTAt                                                                       <  3:305747/139‑1 (MQ=255)
                                                           tACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTAtt                                                                      >  8:335271/1‑139 (MQ=255)
                                                                cGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTTCCGGCATCATATAAGAAAACCCAAAACAAACGGAACAGGATATTATTTCTAttttttt                                                                 >  4:349641/1‑136 (MQ=255)
                                                                cGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGGCGGCATCAAAGAAGAAAAGCGCACAAAACCGGCACAGGAAATTATTTGTAtttttat                                                                 >  6:238992/1‑136 (MQ=255)
                                                                cGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGCCGGGCTCATTGAAGAAAAGCGACACAAACCGGGACAGGATATTATTTCTATTTTAGt                                                                 >  7:10618/1‑139 (MQ=255)
                                                                  aaGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGCCGGCAGCCTAGCAGACAAGCGAAACAAACCGCAACACGCTATTATGTCTATTTTAGttg                                                               >  2:62354/1‑137 (MQ=255)
                                                                  aaGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGGTTTTGCCGGCATGCTAGAAGAAAAGCGAACCAAACCGGAACAGGATATTATTTCGATTTTAGTgg                                                               >  2:216033/1‑139 (MQ=255)
                                                                     gATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATagaag                                                                                                            <  8:271972/91‑1 (MQ=255)
                                                                     gATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATagaag                                                                                                            >  7:271972/1‑91 (MQ=255)
                                                                       tAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTTGCGGCATCATACAAGAAAAGCCAACCAAACCGGAAAAGGTTATTTTTTCTATTTTAGGGGAAGCg                                                          >  7:32598/1‑139 (MQ=255)
                                                                        aaaaGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATATAAGCAAAGCGAAAAAAACCGGAACAGGATATTATTTATATTTTAATGGAAGcc                                                          >  8:117797/1‑137 (MQ=255)
                                                                          aaGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGGTTTTTGCCGGCATCATAGAATAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGGAGCgg                                                         >  7:1823/1‑137 (MQ=255)
                                                                               aagaagCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGaaa                                                   <  7:335271/138‑1 (MQ=255)
                                                                               aagaagCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGaaa                                                   <  3:349641/138‑1 (MQ=255)
                                                                               aagaagCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAAc                                                  <  3:80504/139‑1 (MQ=255)
                                                                                         aaaaGCCTTTTTGGAAGAACTCGATAAGTGTGAGGAAGAACGGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAg                                        <  6:226947/139‑1 (MQ=255)
                                                                                                              agaTAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCTTCATAGAAGAAAATCGAAACAAACCGTAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTCCGGTGAAGAGCTGAtt                   <  7:133345/139‑1 (MQ=255)
                                                                                                              agaTAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTCCGGTGAAGAGCTGAtt                   <  3:9558/139‑1 (MQ=255)
                                                                                                                         aGGAAGAACTGGCCGCGTTTTTTGCCGGGAGCATAGAAGAAAAGCGAAAAAAACCCGAACAGGATATTATTTATATTTTATTTGAAGCGGGAGAAACACGCGCGAAGCGGTCCGGTGAAGCGCTGATTACGGTTTGCAc        >  6:371442/1‑139 (MQ=255)
                                                                                                                         aGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATATAAGAAAAGCGCAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTCCGGTGAAGAGCTGATTCCGTTTTGCAc        >  7:278782/1‑139 (MQ=255)
                                                                                                                         aGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACACACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTCCGGGGAAGAGCTGATTCCGGTTTGCAc        >  8:319626/1‑139 (MQ=255)
                                                                                                                         aGGAAGAACTGGCCGCGGTTTTTGCCGGCATAATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTGCGGTGAAGAGCGGATTGCGGTTTGCAc        >  5:150995/1‑139 (MQ=255)
                                                                                                                               aaCTGGCCGCGGTTTTTGCCGGCATCATAGAAGAAAAGCGCAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGCAGCTGTCCGGGGAAGCGCTGATTCCGGTTTGCACgctgct  >  8:323723/1‑139 (MQ=255)
                                                                                                                               aaCTGGCCGCGGTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCCGAACAGGATATTATTTCTATTTTAGTTGAAGCGGAAGCAACATGCGAGAAGCTGTGCGGTGAAGAGCTGAATCCGGTTTTCACgctgct  >  2:350312/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                               
CTGGGAGTGCCTTCAGCGCATATGGAACAGTTTAAAGCATGGTCTGATCTTCTGGTCAGTACACCGAAGGATAAAAGTGAAGAAGCTGAAAAAGCCTTTTTGGAAGAACGAGATAAGTGTGAGGAAGAACTGGCCGCGTTTTTTGCCGGCATCATAGAAGAAAAGCGAAACAAACCGGAACAGGATATTATTTCTATTTTAGTGGAAGCGGAAGAAACAGGCGAGAAGCTGTCCGGTGAAGAGCTGATTCCGTTTTGCACGCTGCT  >  BS168_reference_genome/1291785‑1292050

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: ‑