Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,388,322 T→G 46.5% 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
*ACB1221,388,3220TG46.5% 9.4 / 5.8 31L336V (TTG→GTG) ruvB_1Holliday junction ATP‑dependent DNA helicaseRuvB
Reads supporting (aligned to +/- strand):  ref base T (14/2);  new base G (13/2);  total (27/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCCACA  >  ACB122/1388199‑1388454
                                                                                                                            |                                                                                                                                    
cATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCCGCCTGCCTACCGGCACTTTGGGCGGAATATTGCCGGGCGCGTGGGGGCgggg                                                                                                                           >  3:185782/1‑135 (MQ=255)
     cATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCCCCCGGCCTGCCGACCTGCCCTTTGGGCGGCAAATTACCGGGCGCTGGGGGGAGGGGGGCGa                                                                                                                      >  1:919/1‑135 (MQ=255)
          aaGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCCCGCGGCCGCGTGGGCCCCCGGCCTGCCTACCTGCCCTTTGGGCTGAATATTCCCGGGCGCGTGGGGGCGGGGGGCGAgtttt                                                                                                                 >  7:99161/1‑135 (MQ=255)
          aaGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCCCTTTGGGCTGAATATTGCCGGGCGCGTGGGGGCGGGGGGCGAttttt                                                                                                                 >  4:410258/1‑135 (MQ=255)
          aaGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCCGCATGGCTACCTGCCCTTTTGGCTGAATATTCCCGGGCGCGTGGGGGAGGGGGGCGAttttt                                                                                                                 >  7:64068/1‑135 (MQ=255)
          aaGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCCGCATGCCTACCTGCACTTTGGGCTGAATATTACCGGGCGCGTGGGGGGGGGGGGCGAttttt                                                                                                                 >  7:9475/1‑135 (MQ=255)
          aaGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGGGGGCACCCCCCAAGCCTACCTGGCCGATTGGGCGTATATTACCGGGCGCGTGGGGGGGGGGGGCGCctttt                                                                                                                 >  5:43007/1‑135 (MQ=255)
           aGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGAtttttg                                                                                                                >  7:460137/1‑134 (MQ=255)
                   ttGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGGGGTCCCCCCGCCTGCCAACCTGCACTTTTGGCGGGATATTGCCGGGGGGGGGGGGGGGGGGGGGGCTTTTTTCggcgccc                                                                                                        >  1:218601/1‑129 (MQ=255)
                              cTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCCCCCGGCCTGCCTACCTGCCCTTTGGGCTGAATATTCCCGGGCGCGTGGGGGAGGGGGGCGCGTTTTAAGAGCCCGGCGatgaagt                                                                                              >  3:317736/1‑132 (MQ=255)
                                    cagcagGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCCTGCCGACCTGCCCTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGCTTTTTCCGCGCCAGGCGCTGGATGGCAGCTc                                                                                       >  8:234851/1‑135 (MQ=255)
                                        agGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGCGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGaaa                                                                                   >  1:54410/1‑135 (MQ=255)
                                        agGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCCGCCTGCCAACCTGCACTTTGGGCTGAATATTACCGGGCCCTTGGGGGTGGGGGGCGATTTTTACGCGCCAGGCGAATAATGAAATAACGaaa                                                                                   >  4:261352/1‑135 (MQ=255)
                                                       gcACGCCACGCGGCCGCGTGGTCCCCCGGCATGCCTACCTGCACTTTGGGCTGAATATACCCGGGCGCGTGGGGGTGGGGGGCGCTTTTTGCGCGCCCGGCGCTGAATGCCAGATCGAAAGCGACGTTTTTTGGt                                                                    >  2:44543/1‑135 (MQ=255)
                                                       gcACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTTGGGGAGGGGGGCGATTTTTGCGAGCCCGGCGATGAATTACAGATCGAAAGCGACTTTTTGTGGt                                                                    >  3:325988/1‑135 (MQ=255)
                                                       gcACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTACCGGGCGCGTTGGGGGGGGGGGCGATTTTTGCGCGCCAGGCGAAGTATGAAAGATCGAAAGCGGCATTTTGTGGt                                                                    >  2:363271/1‑135 (MQ=255)
                                                       gcACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCAGCACTTTGGGCGGAATATTCACGGGGGCGTGGGGGAGGGGGGCGAATTTTACCAGCCCGGCGAAGAAAGAAAGATAGAAAGCCACTTTTtttgtt                                                                    >  6:175764/1‑133 (MQ=255)
                                                       gcACGCCACGCGGCCGCGTGGTCACCCGGCATGCCGACCTGCACTTTGGGGTGAATATTCCCGGGCGCGTGGGGGAGGGGGGCGCTTTTTCCGAGCCCGGCGATGCATGACAGATCGCAAGCGCCCTTTTGTGGt                                                                    >  4:461280/1‑135 (MQ=255)
                                                           cccccgcgGCCCGGGGGCCACCCGGCACGCCAACCGGCCAGTTGGCGAGAACATGCCGGGCGGCGTGGGGCAGGGCGCCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTa                                                                <  2:367038/131‑1 (MQ=255)
                                                             cacgcGGCCGCGTGGTCCCCCGGCATGCCTACCTGCCCTTTGGGCTGAATATTCCCGGGCGCTGGGGGGGGGGGGGCGTGTTTTGCGCGGCCGGGGCGGTATTGCCGCTCGAAAGCGACGTTTTGTGGTccgtcc                                                              >  6:273832/1‑131 (MQ=255)
                                                                 cgGCCCCGGGGCCACCCGGCACCCCTCCCGCCACTTCGGGCGCAAGACCCCCGGCCGCGGGGGGGCGGCGCGCGAGTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTg                                                          <  8:9475/135‑1 (MQ=255)
                                                                   gCCGCGTGTTCACGCGGCATGCCTACCGCCAAGTGCGGCGTAATCTCCCCGGGGGGTTGGGGGAGGGGTGGGATGTTTCGGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTg                                                          <  4:80010/133‑1 (MQ=255)
                                                                        gTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCGTGGGGGGGGGGGGCGATTTTTGCGCGCCAGGGGATGAATGACCGATCGAAAGCGCCGTTTTGTGGTCCTACCCCTGGCCTGcc                                                   >  3:117704/1‑135 (MQ=255)
                                                                                 cGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGCCTTTTTGTGGTCCTACCGCTGGCCTGCCCCTGGTggg                                          >  2:248356/1‑133 (MQ=255)
                                                                                     aTGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTg                                       >  2:138173/1‑134 (MQ=255)
                                                                                                 cACTTTGGGCTGAATATTCCCGGGCGCTTGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTa                         >  1:295372/1‑135 (MQ=255)
                                                                                                  aCTTGGGGCTGACTCTTCCCGGGCGGTGGGGGGCGGGCGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACa                            <  7:225630/133‑1 (MQ=255)
                                                                                                      tGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTGTGTGGTACTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAAt                       >  7:383627/1‑134 (MQ=255)
                                                                                                          cTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATgaag                   >  8:204118/1‑134 (MQ=255)
                                                                                                                          gcTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCcaca  >  7:458555/1‑135 (MQ=255)
                                                                                                                          gcTTGGGGGAGGGGGGCGATTTTTACGAGCCAGGCGATGAATGACAGATCGAAAGCGAATTTTTGTGGTCCTACCGCTGGCATGCCCCTGGGGGGCTGGCAATAAGCCATTAATTAAGAAAAAACAGTTGCcaga  >  7:397004/1‑133 (MQ=255)
                                                                                                                            |                                                                                                                                    
CATACCATCGAAGATGTGCTTGAGCCGTACCTGATTCAGCAGGGCTACATCATGCGCACGCCACGCGGCCGCGTGGTCACCCGGCATGCCTACCTGCACTTTGGGCTGAATATTCCCGGGCGCTTGGGGGAGGGGGGCGATTTTTCCGAGCCAGGCGATGAATGACAGATCGAAAGCGACTTTTTGTGGTCCTACCGCTGGCATGCCCCTGGTGCGCTGGCAATAAGACATTAATGAAGAAAAAACAGTTGCCACA  >  ACB122/1388199‑1388454

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: