Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 575,740 T→G 23.9% F186V (TTT→GTT)  selB → Selenocysteine‑specific elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440575,7400TG23.9% 67.9 / 3.4 48F186V (TTT→GTT) selBSelenocysteine‑specific elongation factor
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base G (11/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 2.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGCAGTTTCCGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCTGCATGC  >  KT2440/575609‑575861
                                                                                                                                   |                                                                                                                         
cgcgGCAGTTTTCGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGGTTg                                                                                                                        >  4:220295/1‑135 (MQ=255)
cgcgGCAGTTTCCGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGGTTg                                                                                                                        >  2:229860/1‑135 (MQ=255)
 gcgGCAGTTTCCGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGGGCGGTTGc                                                                                                                       >  2:113958/1‑135 (MQ=255)
     cAGTTTCCGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGtt                                                                                                                   >  5:142523/1‑135 (MQ=255)
           ccGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGGTTTCGCCTGGCCGTCGCCCGTGCGGTTGCGGTTACCGgc                                                                                                             >  7:251509/1‑135 (MQ=255)
             tcTGTCGAGCATCAGCGGGGAGGGCTGGTACGCGTTGCGTAAGGCCCTGCTGCAGGACGAAGCAGGTAAACGGCAGCGTAGTGTCGTCGGTGGTTTTCGACTGGCCGTCGTCCGTGCTTTTGCTGTTACCGgcg                                                                                                            <  6:20312/133‑1 (MQ=255)
             gCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGgcg                                                                                                            <  4:294847/134‑1 (MQ=255)
                         cAGCGGGGAGGGCGGGGACGCGTGGCGTCTGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTggtggt                                                                                               <  5:110416/135‑1 (MQ=255)
                             ggggCGGGGGGGGACGGGTTGCTTCTGGCCGGGCTGGCGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGGGCGTTTGCGGTTACCGGCGCCGGGGTGGTGGTCAcc                                                                                           <  1:107128/135‑1 (MQ=255)
                                   ggggTGGACGCGTTGCGTCAGGCCCGGGTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGGCGGGGGGTGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGGGTGGGGGTCACcggcacg                                                                                     <  4:86179/132‑1 (MQ=255)
                                   ggCGTGCACGCGTGGCGTCAGGCCCTGCGGGAGGCCGAAGCAGGTAGACGGCAGCGTAGGGGCGGCGGTGGTTTGCGCCTGGCCGGCGACCGTGCGGTTGCGGTTACCGGCGCCGGGGTGGTGGTCACcggcacg                                                                                     <  5:175916/135‑1 (MQ=255)
                                    gCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGGTTTCGCCTGGCCGTCGACCCTGGGTTTGCGGTTACCCGCCCCGGTGTGGGGGTGACccgcaccg                                                                                    >  2:139898/1‑133 (MQ=255)
                                      gTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTGGGGGGGTTTGGGCTGGCCGTCGCCCGGGGGGTTGGGGGTTCCGGCCCCGGGGGGGGGGGACCcggcccggca                                                                                  >  2:367218/1‑135 (MQ=255)
                                         gccgcgTTGCGTCAGGCCCTGCGGGAGGCCGACGCAGGTACACGGCAGCGTAGTGTCGGCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGGGTGGGGGTCACCGGCACGGCACTg                                                                               <  8:59270/133‑1 (MQ=255)
                                                tGCGTCAGGCCCTGCTGGCGGCCGAAGCATGTACACGGCCGCGTAGTGTCGTCGGTGGTTTTAGCCTGGGCGTAGCCCGTGCGTTTGGGGTTTCCGGCGCCCGTGTTGTTGTACCCCGCACGGCCCGTGCcgggc                                                                        >  2:1370/1‑135 (MQ=255)
                                                tGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGGGGTTTTAGCCTGGCCGTCGCCCCTGGGGTTGCGGTTTCCGGCGCCGGTGTGGTGGTCACCGGCACCGCCCAGGCcgggc                                                                        >  7:370476/1‑135 (MQ=255)
                                                 gcgcCAGGCCCTGCTGGCGGCCGAAGCAGGTCCACGGCAGCGTAGTGGCGTCGGTGGTTTTCGCCTGGCCGGCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGGCACCGGCACGGCACTGGCcgggcg                                                                       <  6:381331/135‑1 (MQ=255)
                                                      aGGCCCGGCGGGAGCCCGACGCAGGAACACGGCAGCGTAGGGTCGTCGGGGGTTTTCGCCGGGCCGGCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGGGGGGGGCACCGGCACGGCACTGGCCGGGCGGGTCa                                                                  <  3:325194/135‑1 (MQ=255)
                                                               tgGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCCTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCCGTGCGGGTg                                                         >  3:77736/1‑135 (MQ=255)
                                                               tgGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGGGGTTTTCGCCTGGCCGTCGACCCGGCGGGTGCGGTTACCGGCGCCGGGGTGGTGGTCCCCGGCACGGCACTGGCCGGGCGGGTCCGGGCGGGTg                                                         >  1:88206/1‑135 (MQ=255)
                                                               tgGAGGCCGAAGCAGGTACACGCCAGCGTAGTGTCGTCGGTGGGTTTCGCCTGGCCGTCGCCCGGGCGGTTGCGGGTAACGGCGCCCGGCTGGGGGTCACCGGCACCGCCCTGGCCCGGCGGGTCAGGTCGGGt                                                          >  8:239441/1‑134 (MQ=255)
                                                                 gAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGAt                                                       <  8:68766/135‑1 (MQ=255)
                                                                  cggCCGACGCAGGTACACGGCAGGGTAGGGGCGTCGGTGGTTTTCGCCGGGCCGGCGACCGTGCGTTTGCGGTTACGGGCGCCGGTGTGGGGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATa                                                      <  5:13911/134‑1 (MQ=255)
                                                                  aGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATa                                                      <  6:142523/135‑1 (MQ=255)
                                                                       gaagaaGGTACACGGCAGCGTTGTGTCGGCGGTGGTGTTCGCCTGGCCGTCGCCCGTGCGTTTGCGGTTACCGGCGCCGGTGGGGTGGTCACGGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACgttg                                                 <  2:103571/135‑1 (MQ=255)
                                                                       gAAGCAGGTACACGGCATCTTAGTGTCGTCGGTGGTTTTCGCCTGGCCGGCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACgttg                                                 <  5:50874/135‑1 (MQ=255)
                                                                         aGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACgttgtt                                               >  2:97263/1‑135 (MQ=255)
                                                                           cAGGTACACGGCAGCGTAGTGTCGTCGGTTGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTg                                                         <  2:251931/123‑1 (MQ=255)
                                                                           cAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTGAGTGCGGGTGATACGTTGTTGc                                             >  7:392668/1‑135 (MQ=255)
                                                                           cAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTg                                                         >  1:251931/1‑123 (MQ=255)
                                                                           cAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGGTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGc                                             >  1:208516/1‑135 (MQ=255)
                                                                               tccacGGCAGCGTTGTGTCGTCGGTGGTTTTCGCCTGGCCGGCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGCCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTggg                                         <  2:1049/133‑1 (MQ=255)
                                                                                  cgggCAGCGTAGGGTCGGCGGTGGTTTTCGCCTGGCGGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGGCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCaa                                      <  8:162050/133‑1 (MQ=255)
                                                                                    ggCAGCGTAGTGGCGTCGGTGGTTTTCGCCTGGCCGGCGACGGGGCGTGGGGGGTTACCGGCGCCGGTGGGGTGGTCACCGGCACGGCACTGGCCGGGCGGGGCAGTGCGGGTGATACGTTGTTGCTGGGCAAgg                                    <  7:378414/135‑1 (MQ=255)
                                                                                         cGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGGTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGGGGGTGTTTCGTTGTTGCAGGGCAAGGCCgg                                >  1:270130/1‑134 (MQ=255)
                                                                                            aGTGTCGTCGGTGGTTTTAGCCTGGGCGTCGACCGGGGGGTTTCGGTTTACGGCGCCGGGGTGGTGGGAACCCGGCCCGCACAGGCCGGGGGGGGGGGTGCGGGGGATTCGGTTTTTCTGGGGAAGGGCGGCCAg                            >  2:388679/1‑135 (MQ=255)
                                                                                               gtcggcGGTGGTTGGCGCCTGGCCGCCGACCGGCCGTTGGGGGTTCCCGGCGCCGGGGTGGTGGCCCCCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAggtg                         <  1:379637/135‑1 (MQ=255)
                                                                                                 cgtcgGTGGTTTTCGCCTGGCCATCGACCGTGCGTTTGCGGTTACCGGCGCCGGGGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAggtggt                       >  5:301542/1‑135 (MQ=255)
                                                                                                      gtggtTTTCGCCTGGCCGTCGCACGTGCGTTTGCGGTTACCGGCGCCGGGGTGGTGGTCACCGGCACGGCACTGGCAGGGCGGGTCAGGGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGcgcg                  >  4:269122/1‑135 (MQ=255)
                                                                                                      gtggtTTTCGCCTGGCCGTCGACCGTGCCTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGcgcg                  >  1:11154/1‑135 (MQ=255)
                                                                                                       tgggTTTCGACTGGCCGTCGAGCGGGCGGTTGAGTGTACCGGCGCCTGTGGGGTGGTCACCGGCACGGCACTTGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGAAAGGCCGGCAAGGTGGTTCgcgt                 <  8:326879/135‑1 (MQ=255)
                                                                                                       gggtTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGGGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCgcgt                 <  1:139898/134‑1 (MQ=255)
                                                                                                               gCCGGGCGGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGGGGTGGGCCCCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTc         <  4:12309/135‑1 (MQ=255)
                                                                                                                ccTGGCCGTCTACGGTGCGTTTGGGGTAACCGGGGCCGGTGTGTTGCTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCt        <  1:360376/135‑1 (MQ=255)
                                                                                                                    gCCGGCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGGCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCTGCAt    <  1:113958/135‑1 (MQ=255)
                                                                                                                       gTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCTGCATgc  >  6:247557/1‑134 (MQ=255)
                                                                                                                                 cGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCTGCAt    <  6:366564/122‑1 (MQ=255)
                                                                                                                                 cGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCTGCAt    >  5:366564/1‑122 (MQ=255)
                                                                                                                                   |                                                                                                                         
CGCGGCAGTTTCCGCTGTCGAGCATCAGCGGGGAGGGCGTGGACGCGTTGCGTCAGGCCCTGCTGGAGGCCGAAGCAGGTACACGGCAGCGTAGTGTCGTCGGTGGTTTTCGCCTGGCCGTCGACCGTGCGTTTGCGGTTACCGGCGCCGGTGTGGTGGTCACCGGCACGGCACTGGCCGGGCGGGTCAGTGCGGGTGATACGTTGTTGCTGGGCAAGGCCGGCAAGGTGGTGCGCGTGCGTGGTCTGCATGC  >  KT2440/575609‑575861

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