Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,559,702 A→C 35.7% V146G (GTC→GGC)  mutL_2 ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,559,7020AC35.7% 34.4 / 4.0 39V146G (GTC→GGC) mutL_2DNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base C (2/13);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCC  >  ACB122/5559568‑5559833
                                                                                                                                      |                                                                                                                                   
cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccga                                                                                                                                     >  6:105051/1‑135 (MQ=255)
cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccga                                                                                                                                     >  2:17970/1‑135 (MQ=255)
cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTACCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGCGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccga                                                                                                                                     >  2:375650/1‑135 (MQ=255)
 gACATCCAACCGGGCCCGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccgac                                                                                                                                    <  5:254924/135‑1 (MQ=255)
   cATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGGTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccgaccg                                                                                                                                  <  6:26156/135‑1 (MQ=255)
     tCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGCCCggg                                                                                                                                >  5:44258/1‑135 (MQ=255)
         accGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTgcg                                                                                                                            <  5:364541/133‑1 (MQ=255)
                 agcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGc                                                                                                                    >  4:272987/1‑135 (MQ=255)
                 agcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGc                                                                                                                    >  2:38641/1‑135 (MQ=255)
                      cagcCGCCGGATCCCTTCCTGCAAGTGATCGAATTCGGTTATCTCGGCTTTGCGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAc                                                                                                               >  1:108792/1‑135 (MQ=255)
                       agcCGCCGGATCACTTCCTGCACGTGATCGAATGCGGTTTTCGCGGCTTGGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAccc                                                                                                             <  7:234426/136‑1 (MQ=255)
                        gccgccGGATCACTTCCTGCAGGTCATCTAATGCGGTTTTCTCGGCTTTGAGCAACTGGCGCCTGGCGGGGGTATTCAACAACAAACCGCCCACTCCCCCCGCGGTCCCGCCCGGGGGCGCCGCCGGCTCCAccc                                                                                                             <  4:21830/135‑1 (MQ=255)
                          cgccgGATCACTTCCTGCATGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGCTCAAATCGCGCACTTCCACCGAGGTACCGCCCGGGTGCGCCGCCGGCTGCCCCcgcg                                                                                                          >  8:246229/1‑136 (MQ=255)
                           gccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCGGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcgg                                                                                                         >  8:180956/1‑136 (MQ=255)
                           gccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAACTCGCGCACTGCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCCCCcgcg                                                                                                          >  3:23754/1‑135 (MQ=255)
                           gccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCCCCcgcg                                                                                                          >  1:173843/1‑135 (MQ=255)
                           gccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCcgcg                                                                                                          >  3:268973/1‑135 (MQ=255)
                           gccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcgg                                                                                                         >  8:176240/1‑136 (MQ=255)
                                    cTTCCTCCAAGGGATCGAATTCGGTTTTCTCGGCTTGGAGGAACTTGCGCCTGGCCGGGGTATTGAATAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTc                                                                                                <  8:193969/136‑1 (MQ=255)
                                                     aaTTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCGTGGGTTTc                                                                                 >  1:103263/1‑134 (MQ=255)
                                                             tttCTCGGCTTTGAGTACCTGGCGCGGGGCGGGGGATTTGAATAAAAAATCGCGCACTCCCACCGAGGTACCGCCCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCa                                                                        <  1:13955/135‑1 (MQ=255)
                                                             tttCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGCCCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCa                                                                        <  4:137161/135‑1 (MQ=255)
                                                                  cGGCTTTGCGCAACTTGCGCCTGGCGGGGGTATTGAACAACAACTCCCGCACTCCCACCGCGGTACCCACCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTc                                                                  <  7:176240/136‑1 (MQ=255)
                                                                  cGGCTTTGAGCAACTTGCGCCGGGCGGGGGTATTGACCAACAACTCGCGCACTCCCCCCGCGGTCCCGACCGGGGGCGCCGCCGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTc                                                                  <  7:180956/136‑1 (MQ=255)
                                                                  cGGCTTGGAGCAACTGGCGCCGGGCGGGGGAATGCAACACCAAAGCGCGCACTCCCCCCGAGGAACCGCCCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTc                                                                  <  7:171750/136‑1 (MQ=255)
                                                                   ggCTTTGCGTAACTGGCGCCTGCCGGGGGTATTGAAGAACAAATCGCCCACTTCCCCGGAGTTCCCGCCCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTc                                                                  <  5:209774/135‑1 (MQ=255)
                                                                   ggCTTTGCGGAACTGGCGCCGGGCGGGGGAAGTGAACAACAACCCGCGCACTTCCACCGCGGACCCGCCCGGGTGCGCCGCCGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTc                                                                  <  1:330148/135‑1 (MQ=255)
                                                                   ggCTTTGAGGAAGTTGCGCCGGGCCGGGGTATTGACCAACAAACCGCGCACTCCCCCCGAGGAACCGACCGGGGGCGCCGCCGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTc                                                                  <  1:340268/135‑1 (MQ=255)
                                                                             aactgccgcCGGGCCGGGGTATTAAAGAACAAATCGCCACCTTCCCCCGAGGACCCCCCCGGGGCCGCCGCCGGCGGCCCCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTcactggcactgg                                                       <  7:299743/130‑1 (MQ=255)
                                                                                  gcgagtGGCCGGGGTATTCAAGAACAAAGCGCGCACTTCCACCGAGGTACCGCCCGGGGGCGCCGCGGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGt                                                   <  8:86274/130‑1 (MQ=255)
                                                                                   ggcgtGGCGGGGTTTTTTAACAACAAATCCCGAAATACCCCGGAGGTACCGCCCGGGTCCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTa                                                  <  8:295234/131‑1 (MQ=255)
                                                                                    gcgtGGCGGGGGTATTCAACAACAACCCCCCAACTTCCCCCGTGGTACCGCCCGGGGGCGCCGCGGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTAc                                                 <  5:391242/132‑1 (MQ=255)
                                                                                        ggCGGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAg                                            <  8:87621/136‑1 (MQ=255)
                                                                                                     acgacaaaaGCGCGCACTCCCACGGCGGTACCGACCGGGGGCGCCGCCGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAggtcagggtcag                                 <  4:17024/128‑1 (MQ=255)
                                                                                                     aataacaacTCCCGAACTTCCACTTTGGAACCCCCCGGTGGCGCCGGGGCCTGCACCAGCGGCGTCATGTCGCGGCCTTCTGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAggtcagggtcag                                 <  2:87482/131‑1 (MQ=255)
                                                                                                      agaaCAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGgcg                              >  8:86815/1‑136 (MQ=255)
                                                                                                         ccAACTCGCCCATTACCCCCCAGGAACCGACGGGGTGCGCCGCCGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGgcgcg                            <  6:203828/134‑1 (MQ=255)
                                                                                                           aaaTCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGcc                          >  2:66990/1‑135 (MQ=255)
                                                                                                             aTCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCAc                        >  2:236645/1‑135 (MQ=255)
                                                                                                              tCGCGCACTTCCACCGAGGTACCGACCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCAcc                       <  4:83089/135‑1 (MQ=255)
                                                                                                              tCGCCCACTTCCACCGAGGTACCGCCCGGGGGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCAcc                       <  1:44468/135‑1 (MQ=255)
                                                                                                               cgcgCACTCCCCCGGCGCTCCCGAGCGGGGGGCCCGCCGGGTCCACCCGCGGCGTCATGTCGCGGCCTTCGGTTGCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCg                      <  5:66177/135‑1 (MQ=255)
                                                                                                                gcgcACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCGTCGGTTTCTACATGCCACGCCGCCCAGGCCCTGGCCGTAACCGCGGTCACGGTCACGCGCGCCCCCgc                     >  5:74462/1‑134 (MQ=255)
                                                                                                                    aCTGCCCCCGAGGGACCGCCGGGGGGGGCCGCCGGCGGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTg                 <  3:272987/135‑1 (MQ=255)
                                                                                                                                   ccgaccgGGTGCGCCGCCGGCTGCCCCCGCGGCGTCATGTCGCGGCCGTGGGTTTGTACTTGCCAGGCCTCACGGGCCCTGGCGGTACGCGAGGGCAGGGTCAGGCGCGCCCCCGCACTGATCGCGGCCAACGcc  >  6:70602/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCC  >  ACB122/5559568‑5559833

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