Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 965,051 T→G 14.8% L342R (CTT→CGT)  ssuD → FMNH2‑dependent aliphatic sulfonate monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome965,0510TG14.8% 95.5 / 3.5 43L342R (CTT→CGT) ssuDFMNH2‑dependent aliphatic sulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (19/16);  new base G (5/2);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAGGGGCAGGAACAGCTCTTGTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGCCCCCCGAACA  >  BS168_reference_genome/964914‑965185
                                                                                                                                          |                                                                                                                                       
cggaggGTCAGGAACAGCTCTTGTCGGTGATCAGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCtt                                                                                                                                        <  8:81803/138‑1 (MQ=255)
     gggCAGGAACAGCTCTTGTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCAGGCTATCCTCATTTAGAAGAA‑GCGTATTATTTTGCTGAGCTTGTAtt                                                                                                                                   >  4:330214/1‑139 (MQ=255)
      ggCAGGAACAGCTCTTGTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGCATCGTGTATATTCTCCGGCGTTGCGCATTTAGAAGCA‑GCGTATTATTTTGCAGAGCGTGTAttt                                                                                                                                  >  3:101162/1‑139 (MQ=255)
                 ctctTGTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGcc                                                                                                                      <  3:86629/139‑1 (MQ=255)
                 ctctTGTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTGTATATTTAACCGCTATACGCAATTTAAAGAA‑TCGGAATAATTTTCAAATCATGTATTTACGA‑TTTTGcc                                                                                                                      >  4:86629/1‑139 (MQ=255)
                      gTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTTTCCCCATTTAGAAGAA‑GGGTATTATTTTGCAGAGCTTGTATTTACGC‑GTTTGCCGttt                                                                                                                  >  8:278983/1‑138 (MQ=255)
                           tGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGCGCTTGTATTTCCGC‑TTTTg                                                                                                                        >  8:409936/1‑127 (MQ=255)
                           tGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTg                                                                                                                        <  7:409936/127‑1 (MQ=255)
                              tCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGAt                                                                                                         <  1:194784/139‑1 (MQ=255)
                                 gCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTTGAAGAA‑GCGGATTATTTTGCAGAGCTGGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCgg                                                                                                      >  8:389361/1‑139 (MQ=255)
                                        aTTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGaa                                                                                               <  1:31685/139‑1 (MQ=255)
                                             cGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCtt                                                                                          <  1:26033/139‑1 (MQ=255)
                                                                  ggCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCg                                                                     <  7:278983/139‑1 (MQ=255)
                                                                  ggCACTCGGCATAGAATCGTGTTTATTATGCGGCGATACGCATTTTTAAGAA‑GGGGTTTATTTTTACGAGCGTGTTTTTACGCGATTTGCGGTTGGAAAA‑GAACGGGCCCGCAACCGTGAGACAAAAAGGGGAGAAGCg                                                                     >  2:372811/1‑139 (MQ=255)
                                                                  ggCACTCGGCATAGAATCGTGTATATTGTCCGGCGATCCGCATTTAGAAGAA‑GCGTATTTTTTTGCAGAGCGTGTATTTGCGC‑TTTTGCCGTTTTCAAATGGTCGGACGAGAAAGCGTCAGCACAAAAGGGGCGAAGCg                                                                     >  3:160232/1‑139 (MQ=255)
                                                                  ggCACTCGGCATAGAATCGTGTATATTATCCCGCTATCCCCATTTTGAAGAA‑GCGGATTTTTTTTCAGAGCTTGTATTTGCGC‑TTTTGGGGTTTGTAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGTGATGCg                                                                     >  3:270195/1‑139 (MQ=255)
                                                                       tCGGCATAGAATCGTTTTTATTGTGCGGCGTTCCGCATTTTTAAGAAGGGGTATTTTTTTG‑AGAGGTTGTATTTACGC‑GTTTGGGGGTTGTAAATTTTAGGGGCAGAAAGCGTAAGAACAAAAAGGGTGAAGCGGTAgg                                                                >  4:77646/1‑139 (MQ=255)
                                                                       tCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTTCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAgg                                                                >  8:19369/1‑139 (MQ=255)
                                                                       tCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTAtt                                                                                                                                   >  6:359390/1‑73 (MQ=255)
                                                                       tCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTAtt                                                                                                                                   <  5:359390/73‑1 (MQ=255)
                                                                       tCGGCATAGAATCGTTTATATTATCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCGTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAgg                                                                >  4:364023/1‑139 (MQ=255)
                                                                                     tttATATTCAACGGCTATACGCATTTAGAAGAA‑GCGTATTATTTTTCAGAGCTTGTAATTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAACCATCAGAACAAAAGGGGAGAAGCGGTAGGAAACACATATTTTg                                                  >  3:188567/1‑139 (MQ=255)
                                                                                     tgtttatTCTCCGGGTATCCGCATTTTGAAGGA‑GCGGATTATTTTGCAGCGCGTGTATTTCCGC‑TTTTGCCGGTTGTAAATGAGCGGACGTGGAAGCGTCAGTACAAAAGGGGGGAAGCGGGTGGCAACCCGTATTTTg                                                  >  6:211116/5‑139 (MQ=255)
                                                                                              tcCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAagag                                          >  3:152727/1‑138 (MQ=255)
                                                                                              tcCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAagag                                          >  4:202167/1‑138 (MQ=255)
                                                                                                cgtcTATCCGCAATTAGAAGAA‑GCGTATTATTTTTCATAGCGTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGaa                                        <  6:137542/135‑1 (MQ=255)
                                                                                                     aTCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATg                                   >  6:35859/1‑138 (MQ=255)
                                                                                                      tCCTCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAACGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCa                                 <  1:372811/139‑1 (MQ=255)
                                                                                                          cATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAg                             >  7:95536/1‑139 (MQ=255)
                                                                                                                          tattatTTTGCAGAGCTTGTAATTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACa              <  4:270195/139‑1 (MQ=255)
                                                                                                                            ttattTTGCAGAGCTTTGATTTCCGA‑TTTTGCCTTTTGAAAATGATCGTAAGAGAAAGCTTCAGAACAAACGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAg             <  6:367886/138‑1 (MQ=255)
                                                                                                                            ttattTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAg             <  3:168902/138‑1 (MQ=255)
                                                                                                                            ttattTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGc            >  5:331669/1‑139 (MQ=255)
                                                                                                                               ttttGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTGATTTTGTCAAAGAGAAAAATGCATGAGTTTCTCTAATGTACGGcccc         >  2:287436/1‑139 (MQ=255)
                                                                                                                                 ttGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGcccccc       >  5:367886/1‑139 (MQ=255)
                                                                                                                                  tGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATa                >  7:158181/1‑129 (MQ=255)
                                                                                                                                  tGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATa                <  8:158181/129‑1 (MQ=255)
                                                                                                                                  tGCAGAGCTTGTATTTAAGC‑TTTTGCCGTTTTAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAATCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGccaccc       >  7:412990/1‑138 (MQ=255)
                                                                                                                                    cAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAAc                       <  7:116767/120‑1 (MQ=255)
                                                                                                                                    cAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAAc                       >  8:116767/1‑120 (MQ=255)
                                                                                                                                    cAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGCCCCCCGaa    <  3:202167/139‑1 (MQ=255)
                                                                                                                                      tagCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGCCCCCCGAACa  <  1:287436/138‑1 (MQ=255)
                                                                                                                                      gagCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGCCCCCCGAACa  <  7:19369/139‑1 (MQ=255)
                                                                                                                                      gagCGTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGCCCCCCGAACa  <  3:364023/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
AGGAGGGGCAGGAACAGCTCTTGTCGGTGATCCGCAGACGATTGCCGACCGGATAGCGGAGTACCAGGCACTCGGCATAGAATCGTTTATATTCTCCGGCTATCCGCATTTAGAAGAA‑GCGTATTATTTTGCAGAGCTTGTATTTCCGC‑TTTTGCCGTTTGAAAATGATCGGACGAGAAAGCTTCAGAACAAAAGGGGAGAAGCGGTAGGAAACACTTATTTTGTCAAAGAGAAAAATGCATGAGTTTAACTAAAATACAGCCCCCCGAACA  >  BS168_reference_genome/964914‑965185

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