Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,388,322 T→G 36.8% L336V (TTG→GTG)  ruvB_1 → Holliday junction ATP‑dependent DNA helicaseRuvB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,388,3220TG36.8% 22.6 / 2.7 30L336V (TTG→GTG) ruvB_1Holliday junction ATP‑dependent DNA helicaseRuvB
Reads supporting (aligned to +/- strand):  ref base T (15/2);  new base G (7/4);  total (23/7)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGCGTCATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCCACAG  >  KT2440/1388192‑1388455
                                                                                                                                  |                                                                                                                                     
agagCGTCATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGGCTACCTGCACTTTTGGCTGAATATTACCGGGGGCTTgggg                                                                                                                                  >  3:213556/1‑136 (MQ=255)
     gTCATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCGTGGGGGAggg                                                                                                                             >  4:11269/1‑136 (MQ=255)
     gTCATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCATGGGGGAgg                                                                                                                              >  6:272523/1‑135 (MQ=255)
                 aaGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCCGCATGGCTTCCTGCCCTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGAttttt                                                                                                                  >  8:59588/1‑135 (MQ=255)
                   gATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCAACCCGCCTGCCTACCTGCACTTTTGGGGGGATATTACCGGGCCCGTGGGGGGTGGGGGCGGTTTTTcc                                                                                                                >  6:18580/1‑135 (MQ=255)
                                     cTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCCGCATGCCTACCTGCACTTTGGGCGGGATATTCCCGGGCGCGTGGGGGAGGGGGGCGCTTTTTACGCGCCCGGCGgtgaatga                                                                                              >  3:88644/1‑135 (MQ=255)
                                     cTGATTCAGCAGGGCTACATCATCAGCACCCCACGAGGCCGAGTGGACACACCGCCTCCCAACCAGCACATCGGGGTGAATATTACCCGGCGCGTGGGGGCAGGTGGCGAATTTTACGACCCACGCGatggatga                                                                                              >  6:159158/1‑135 (MQ=255)
                                               agGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTACGAGCCACGCGATGAATGACAGATAGAAAg                                                                                   >  3:271970/1‑136 (MQ=255)
                                               agGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGCCATGCCTACCTGCACGTTGGGCGGAATAATACCGGGCGCGTGGGGGCGGGGGGCGAATTTTCCGCGCCAGGCGCTGCATGACAGAGCGaaa                                                                                    >  6:300441/1‑135 (MQ=255)
                                                         catGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATTCATGACAGCTCGAAAGCGCCttttt                                                                          >  8:195375/1‑135 (MQ=255)
                                                              gcACGCCACGCGGCCGCGTGGTCCCCCGGCATGCCTACCTGCACTTTGGGCTGAATATTACCGGGCGCTTGGGGGAGGGGGGCGATTTTTACGAGCCAGGCGCTGAATGACAGGTAGAAAGCGCCgtttttttgt                                                                     >  5:172499/1‑130 (MQ=255)
                                                              gcACGCCACGCGGCCGCGTGGTCCCCCGGCATGCCGACCTGCACTTTGGGCTGAATATTGCCGGGCGCGTGGGGGAGGGGGGCGTTTTTTTCGCGCCCGGCGGTGAATGACACCTGGCAAGCGACGTTTTGTGGt                                                                     >  2:17651/1‑135 (MQ=255)
                                                              gcACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCGTGGGGGAGGGGGGCGATTTTTGCGCGCCAGGCGCTGAATGACCGATCGAAAGCGACTTGTTGTGGt                                                                     >  8:200950/1‑135 (MQ=255)
                                                              gcACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTAACGGCACTTTGGGCTGAATATTACCGGGCGCTTGGGGGCGGGGGGCGCTTTTTGCGCGCCCGGCGATGAATGACACCACGGAAACGCCGTTTTGTggg                                                                     >  5:78036/1‑134 (MQ=255)
                                                              gcACGCCACGCGGCCGCGGGGTCACCCGGCCTGCCTAACAGCAATTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGAGTTTTCCCCGCCAGGCCAGTAATGACAAAACGAAAGCGCCGGTTTGTggg                                                                     >  5:58654/1‑134 (MQ=255)
                                                                       gcgcCCGGGGGGTCCCCCGGCATGCCTACCTGCCCTTGGGGCTCACTCTGCCGGGGCGCGGGGGGGCGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTg                                                           <  3:289096/136‑1 (MQ=255)
                                                                               gTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGGGGGGGGCGATTTTTGCCAGCCAGGCCATGAATTCCACAAAGAAACCGCCGTTTTTTGGTCCAACCGCTGGCATGcc                                                    >  6:115611/1‑135 (MQ=255)
                                                                                    cACCCGGCATCCCCCCCGCCCCTTGGGCCGGAATTTTCCCGCGCCCGTGGGGGCGGGGGCCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGt                                              <  4:304146/136‑1 (MQ=255)
                                                                                             gcccccccTGCCCTGTGGGCTGAATCTGCCCGGGCGCGGGGGCGCGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGACAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCa                                     <  4:105039/130‑1 (MQ=255)
                                                                                              gCCGACCTGCACTTTGGGCGAAATCTGCCGGGGCGCTGGGGGCCGGGGGGCCATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCa                                     <  8:196094/135‑1 (MQ=255)
                                                                                                      gCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGCTCGAAAGCGACGTTTTGTGGTCCTACCGCTGGCATGCCCCGGGGGCGCTGGCAATAAGACa                             >  1:141965/1‑135 (MQ=255)
                                                                                                        acgttGGGCTCACTAAGCCCGCGCGGGCGGGGCCGGGGGCCCCTTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACAtt                           <  6:311924/132‑1 (MQ=255)
                                                                                                                     tatTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACATTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGaaaaaa              >  3:79268/1‑135 (MQ=255)
                                                                                                                           cGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGc       <  3:139041/136‑1 (MQ=255)
                                                                                                                           cGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGc       >  4:139041/1‑136 (MQ=255)
                                                                                                                            gggCGCCGGGGGGCGGGGCGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGcc      <  3:11269/136‑1 (MQ=255)
                                                                                                                                gcTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTTATGAAGAAAAAAAAGTTGCCACAg  >  4:36412/1‑136 (MQ=255)
                                                                                                                                gcTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCCACAg  >  4:210148/1‑136 (MQ=255)
                                                                                                                                gcTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTACTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCcaca   >  8:350955/1‑135 (MQ=255)
                                                                                                                                gcTTGGGGGAGGGGGGCGAATTTTCCGACCCAGGCGATGAAAGACACATCGAAAGCGACTTTTTGTTGTGCTACCGCGGGCATGCCCCAGGTGCGCAGGCAATAACTCAATAATGAAGAAAAAACAGTTGcccaa   >  2:40738/1‑132 (MQ=255)
                                                                                                                                  |                                                                                                                                     
AGAGCGTCATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCCACAG  >  KT2440/1388192‑1388455

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