Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,798,113 C→T 7.4% R29Q (CGG→CAG)  tmk_3 ← Thymidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,798,1130CT 7.4% 178.5 / 3.9 54R29Q (CGG→CAG) tmk_3Thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (17/33);  new base T (2/2);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

AGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCG  >  ACB143/3797982‑3798243
                                                                                                                                   |                                                                                                                                  
aGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTaa                                                                                                                                <  4:365286/136‑1 (MQ=255)
   gTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGc                                                                                                                             >  1:7514/1‑136 (MQ=255)
        cAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAg                                                                                                                        <  5:312394/136‑1 (MQ=255)
        cAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAg                                                                                                                        <  3:120430/136‑1 (MQ=255)
                  aCTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTc                                                                                                              <  8:22117/136‑1 (MQ=255)
                  aCTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTc                                                                                                              <  5:78731/136‑1 (MQ=255)
                       aGGTCCTTGCAAGGGCTTCTGGAAAGCTGGTTGACCAGTGCCATGGTCTCGCCCCGGTGATGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCcagca                                                                                                         <  6:390241/136‑1 (MQ=255)
                         gTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCagcagc                                                                                                       <  4:330623/136‑1 (MQ=255)
                                cAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCt                                                                                                <  4:181545/136‑1 (MQ=255)
                                cAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCt                                                                                                <  3:189270/136‑1 (MQ=255)
                                 aaGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCACGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCtt                                                                                               >  3:87788/1‑136 (MQ=255)
                                   gggCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGc                                                                                             <  5:281547/136‑1 (MQ=255)
                                   gggCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGc                                                                                             <  7:167672/136‑1 (MQ=255)
                                   gggCTTCTGGAAAGCTGGTGGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTTAGTGACGGCCTCCAGCAGCGTGGTCTTGc                                                                                             <  6:104460/136‑1 (MQ=255)
                                       ttCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCgg                                                                                          >  2:327321/1‑135 (MQ=255)
                                         cTGGACAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTAc                                                                                       <  2:359479/136‑1 (MQ=255)
                                            gAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCt                                                                                    >  7:293473/1‑136 (MQ=255)
                                            gAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCt                                                                                    >  6:342350/1‑136 (MQ=255)
                                                  tGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCTGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGcc                                                                              >  4:222099/1‑136 (MQ=255)
                                                  tGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCTGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGcc                                                                              >  4:216850/1‑136 (MQ=255)
                                                     ttGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAATGGCGGATCACCTTTTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCAt                                                                           <  8:328058/136‑1 (MQ=255)
                                                     ttGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCAt                                                                           <  5:157663/136‑1 (MQ=255)
                                                     ttGACCAGGGCCATGGGCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGGTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCAt                                                                           <  8:117296/136‑1 (MQ=255)
                                                      tGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCAt                                                                           <  7:257945/135‑1 (MQ=255)
                                                       gACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATc                                                                          >  7:142152/1‑135 (MQ=255)
                                                       gACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATc                                                                          >  2:111232/1‑135 (MQ=255)
                                                               ccATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGaaa                                                                  >  3:285766/1‑135 (MQ=255)
                                                                           cccGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcgggcg                                                     >  5:302933/1‑136 (MQ=255)
                                                                             cGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGt                                                   <  5:255106/136‑1 (MQ=255)
                                                                                tGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCa                                                <  7:249609/136‑1 (MQ=255)
                                                                                 gAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCAt                                               >  5:93023/1‑136 (MQ=255)
                                                                                       cATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGc                                          >  3:44736/1‑135 (MQ=255)
                                                                                        aTTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCaa                                        <  1:73503/136‑1 (MQ=255)
                                                                                           ctttctttCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATAc                                     >  2:107093/1‑136 (MQ=255)
                                                                                            tttctttcGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATAc                                     >  7:358624/1‑135 (MQ=255)
                                                                                                   cGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTgg                              <  2:402250/135‑1 (MQ=255)
                                                                                                   cGCAAAGGCGGATCACCTTGTTCCCGTACAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGTTATAAACAGCGGGCGGTTCTTCGTGCAATACCAGTTgg                              <  7:100548/135‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCTGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  3:216850/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCTGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  3:222099/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGTGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  4:273354/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  8:69864/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  6:41989/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  3:120849/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  2:336749/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  1:388113/136‑1 (MQ=255)
                                                                                                       aaGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCaa                         <  1:327321/136‑1 (MQ=255)
                                                                                                            ggATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGcgcg                    <  1:328771/136‑1 (MQ=255)
                                                                                                              aTCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCggg                  <  4:313350/136‑1 (MQ=255)
                                                                                                                  ccTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGAt                >  6:139728/1‑134 (MQ=255)
                                                                                                                       ttCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGAcc          <  8:142152/135‑1 (MQ=255)
                                                                                                                         cccGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCAt        <  3:300895/135‑1 (MQ=255)
                                                                                                                           cGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGa     >  1:222365/1‑136 (MQ=255)
                                                                                                                                aGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCg  >  5:4321/1‑134 (MQ=255)
                                                                                                                                aGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCg  >  1:376534/1‑134 (MQ=255)
                                                                                                                                   |                                                                                                                                  
AGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCG  >  ACB143/3797982‑3798243

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