Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 88,548 C→T 11.7% C305C (TGC→TGT)  dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome88,5480CT11.7% 94.6 / 4.0 34C305C (TGC→TGT) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (14/16);  new base T (2/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TATGCAAACTGCATCTTGACCGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCACCTCTATTTGCTGCTCACTGCTAGTTTTACGCTGGCAG  >  BS168_reference_genome/88412‑88679
                                                                                                                                        |                                                                                                                                   
tATGCAAACTGCATCTTGACCGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGcga                                                                                                                                   <  1:46434/139‑1 (MQ=255)
tATGCAAACTGCATCTTGACCGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGcga                                                                                                                                   <  3:171953/139‑1 (MQ=255)
  tGCAAACTGCATCTTGACCGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGcgaga                                                                                                                                 <  5:178379/139‑1 (MQ=255)
           cATCTTGACCGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACgagaga                                                                                                                          <  1:46844/137‑1 (MQ=255)
                   ccGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTa                                                                                                                <  8:88593/139‑1 (MQ=255)
                   ccGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTa                                                                                                                <  7:187625/139‑1 (MQ=255)
                   ccGCCTGATTGACTTAAAAGGTGAAAACGTAGAGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTa                                                                                                                <  8:168702/139‑1 (MQ=255)
                     gCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGcga                                                                                                                                   >  5:160833/1‑118 (MQ=255)
                     gCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGcga                                                                                                                                   <  6:160833/118‑1 (MQ=255)
                                           aaaCGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGt                                                                                         <  8:29118/138‑1 (MQ=255)
                                                  gCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGtt                                                                                  >  2:143815/1‑138 (MQ=255)
                                                        agagAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGtt                                                                                        >  2:26176/1‑126 (MQ=255)
                                                        agagAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGtt                                                                                        <  1:26176/126‑1 (MQ=255)
                                                                           cTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACgagaga                                                                                                                          >  6:62025/1‑73 (MQ=255)
                                                                           cTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACgagaga                                                                                                                          <  5:62025/73‑1 (MQ=255)
                                                                             gCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCATAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      >  7:695/1‑139 (MQ=255)
                                                                              cATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      <  8:57714/138‑1 (MQ=255)
                                                                                  gTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCCTTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAACAGAGCTTCAAAATGCAAAAGTAGGa                                                 <  3:76202/139‑1 (MQ=255)
                                                                                  gTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGa                                                 <  1:100112/139‑1 (MQ=255)
                                                                                   tATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      >  2:201859/1‑133 (MQ=255)
                                                                                   tATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      <  1:201859/133‑1 (MQ=255)
                                                                                     tCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      >  6:45081/1‑131 (MQ=255)
                                                                                     tCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      <  5:45081/131‑1 (MQ=255)
                                                                                      cTGAAAGGCGTGAAAGGCAAAGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGataa                                              >  8:119065/1‑138 (MQ=255)
                                                                                        gAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATAc                                           >  6:195771/1‑139 (MQ=255)
                                                                                         aaaGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACt                                          <  5:195771/139‑1 (MQ=255)
                                                                                             gCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCAcc                                      <  7:79932/139‑1 (MQ=255)
                                                                                             gCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCAcc                                      >  1:119223/1‑139 (MQ=255)
                                                                                             gCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCAGTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCAcc                                      >  4:90394/1‑139 (MQ=255)
                                                                                                tGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACgagaga                                                                                                                          <  1:46262/52‑1 (MQ=255)
                                                                                                tGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACgagaga                                                                                                                          >  2:46262/1‑52 (MQ=255)
                                                                                                          tGCGAACGTTCGTAATGATATCAATCATTGTGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      >  3:223109/1‑110 (MQ=255)
                                                                                                          tGCGAACGTTCGTAATGAAATCAATCATTGTGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      <  4:226544/110‑1 (MQ=255)
                                                                                                          tGCGAACGTTCGTAATGAAATCAATCATTGTGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      <  4:223109/110‑1 (MQ=255)
                                                                                                          tGCGAACGTTCGTAATGAAATCAATCATTGTGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAg                                                      >  3:226544/1‑110 (MQ=255)
                                                                                                          tGCGAACGTTCGTAATGAAATAAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCACCTCTATTTgctgct                         <  7:77624/139‑1 (MQ=255)
                                                                                                                    cGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCACCTCTATTTGCTGCTCACTGCTAGt               >  6:116452/1‑139 (MQ=255)
                                                                                                                    cGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCACCTCTATTTGCTGCTCACTGCTAGt               >  1:80852/1‑139 (MQ=255)
                                                                                                                                 atcaTTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCACCTCTATTTGCTGCTCACTGCTAGTTTTACGCTGGCAg  >  6:22929/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TATGCAAACTGCATCTTGACCGCCTGATTGACTTAAAAGGTGAAAACGTAGCGGTAAGAGAAATGAGAAAGCACGCTGCATGGTATCTGAAAGGCGTGAAAGGCAATGCGAACGTTCGTAATGAAATCAATCATTGCGAGACGAGAGAAGAGTTTGTACAGCTTTTAGATGCATTTACGGTTGAAGTTGAGGCAAAAGAGCTTCAAAATGCAAAAGTAGGATAATACTCACCTCTATTTGCTGCTCACTGCTAGTTTTACGCTGGCAG  >  BS168_reference_genome/88412‑88679

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