Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,791,318 C→A 9.3% P159P (CCG→CCT)  nucM ← Nuclease NucM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,791,3180CA 9.3% 116.3 / 4.7 43P159P (CCG→CCT) nucMNuclease NucM
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base A (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CTTGCGGCGGGTAGGTCTTGTTCCAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCGCCGTTGACCTCACCGATGCTGG  >  KT2440/2791194‑2791439
                                                                                                                            |                                                                                                                         
cTTGCGGCGGGTAGGTCTTGTTCCAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCt                                                                                                                 <  5:126019/135‑1 (MQ=255)
      gcggGTAGGTCTTGTTCCAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCg                                                                                                           >  3:46901/1‑135 (MQ=255)
            aGGTCTTGTTCCAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGa                                                                                                     >  1:151407/1‑135 (MQ=255)
                cTTGTTCCAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTc                                                                                                 <  4:46823/135‑1 (MQ=255)
                      ccAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTg                                                                                           >  2:160410/1‑135 (MQ=255)
                      ccAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTg                                                                                           >  7:10172/1‑135 (MQ=255)
                       cAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCTCCTTGAAGTCGACCTg                                                                                           <  5:187118/134‑1 (MQ=255)
                        aGGCCTCGAACCACTGGCGAGCCTGTGTGGACAGGCGCAGGTGGTACTGCTTGCTCATGTAGAAGTCGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTg                                                                                           <  5:107378/133‑1 (MQ=255)
                           ccTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCa                                                                                      >  1:24008/1‑135 (MQ=255)
                           ccTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCa                                                                                      <  8:295058/135‑1 (MQ=255)
                           ccTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCa                                                                                      <  6:281247/135‑1 (MQ=255)
                           ccTCCAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCa                                                                                      <  6:55739/135‑1 (MQ=255)
                            cTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAg                                                                                     >  3:134483/1‑135 (MQ=255)
                            cTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGAGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAg                                                                                     >  7:92700/1‑135 (MQ=255)
                            cTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGAGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAg                                                                                     >  7:87071/1‑135 (MQ=255)
                                aacaacTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAg                                                                                 <  7:250315/135‑1 (MQ=255)
                                         cGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATAc                                                                        <  8:321746/135‑1 (MQ=255)
                                         cGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGACTTGAAGTCGACCTGGGTCAGGCAGCGGCCATAc                                                                        >  4:382769/1‑135 (MQ=255)
                                            tCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGAGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTgg                                                                     <  8:92700/135‑1 (MQ=255)
                                            tCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGAGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTgg                                                                     <  8:87071/135‑1 (MQ=255)
                                              cTGTTTGGACAGGCGCAGGTTTTACTGCTTGCGCATGTAGAAGTAGGTGCGGGCGATCATGCCGCTGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGcc                                                                   <  8:211622/135‑1 (MQ=255)
                                               tGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGcc                                                                   <  8:10172/134‑1 (MQ=255)
                                               tGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCcg                                                                  <  2:151407/135‑1 (MQ=255)
                                                  ttGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCCTGCGGGCcgcgg                                                               >  1:173188/1‑134 (MQ=255)
                                                            gcAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTGGTTGTggga                                                     >  5:122183/1‑133 (MQ=255)
                                                                     aCTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGa                                            <  7:370625/135‑1 (MQ=255)
                                                                       tgcttgctCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGTGTCAGGCCGCGGCCATACTGGCCGCGGTTGTGTGGTAGGCAGCCgc                                            >  6:33275/1‑132 (MQ=255)
                                                                       tgcttgctCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGGTGTTGTGGTTGCCAGCCgca                                           >  2:129637/1‑132 (MQ=255)
                                                                        gcttgctCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTGGGGTAGCCCGCCGaa                                           >  2:356759/1‑133 (MQ=255)
                                                                             ctCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGGTGGTGTGGGAGCCACCCGCAAaagaag                                    >  5:122922/1‑135 (MQ=255)
                                                                              tCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGGCAGGCAGCGGCCATACTGGCCGCGTTGTTATGGTAGCCATCCCAAACTGaact                                   >  7:314124/1‑133 (MQ=255)
                                                                                    aGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCGGGGGCTGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTtgctgc                             <  4:379525/135‑1 (MQ=255)
                                                                                         tAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCg                        <  3:382111/135‑1 (MQ=255)
                                                                                            gTGCGGGCGAGCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGGCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCGcc                      <  5:366753/134‑1 (MQ=255)
                                                                                              gCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCGCCGtt                   <  7:222179/135‑1 (MQ=255)
                                                                                               cGGGCGATCATGCCGCGGATCGACTGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGGTGTTCTGGTAGCCCGCCGAAACTGAACTTTCTGCGGTCGCCGTTg                  >  6:264293/1‑135 (MQ=255)
                                                                                               cGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCcg                                                                  <  6:235136/87‑1 (MQ=255)
                                                                                               cGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCcg                                                                  >  5:235136/1‑87 (MQ=255)
                                                                                               cGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTGTGGTAGCGACCCGAAACTGAAGTTGCTGAGGTCGCCGTTg                  >  7:319965/1‑135 (MQ=255)
                                                                                                    gATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGCCCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGGCCGCCGAAACTGAAGTTGCTGCGGTCGCCGTTGGCCTc             >  1:20767/1‑135 (MQ=255)
                                                                                                            cgcgGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCGCCGTTGACCTCACCGAt       <  4:252619/133‑1 (MQ=255)
                                                                                                            cgcgGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGGCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCGCCGTTGACCTCACCGAt       <  1:160410/133‑1 (MQ=255)
                                                                                                               ggATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGCAACTGAAGTTGCTGCGGTCGCCGTTGCCCTCCCCGCTGCggg  >  5:192641/1‑135 (MQ=255)
                                                                                                                            |                                                                                                                         
CTTGCGGCGGGTAGGTCTTGTTCCAGGCCTCGAACAACTGGCGATCCTGTTTGGACAGGCGCAGGTTGTACTGCTTGCTCATGTAGAAGTAGGTGCGGGCGATCATGCCGCGGATCGACGGGCGCGGCATCACCTTCCTCGCCTTGAAGTCGACCTGGGTCAGGCAGCGGCCATACTGGCCGCGTTGTTCTGGTAGCCAGCCGAAACTGAAGTTGCTGCGGTCGCCGTTGACCTCACCGATGCTGG  >  KT2440/2791194‑2791439

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