Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,397,175 G→C 9.9% L295L (CTG→CTC)  PP_3001 → CAIB/BAIF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,397,1750GC 9.9% 142.8 / 5.5 50L295L (CTG→CTC) PP_3001CAIB/BAIF family protein
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base C (3/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CGGCGACCTGGTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACATGGTCATCG  >  NC_002947/3397028‑3397320
                                                                                                                                                   |                                                                                                                                                 
cGGCGACCTGGTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTgg                                                                                                                                                  <  5:153096/149‑1 (MQ=255)
   cGACCTGGTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCgg                                                                                                                                               >  1:113291/1‑149 (MQ=255)
    gACCTGGTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCggt                                                                                                                                              <  2:145717/149‑1 (MQ=255)
      ccTGGTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCggtgg                                                                                                                                            >  7:113900/1‑149 (MQ=255)
          gTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCCAGGATCGCCTCGCggtgg                                                                                                                                            <  3:163936/145‑1 (MQ=255)
          gTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCCAGGATCGCCTCGCggtgg                                                                                                                                            >  4:163936/1‑145 (MQ=255)
                gccgccCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGAtt                                                                                                                                   <  1:29142/148‑1 (MQ=255)
                 ccgccCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCcgc                                                                                                                                       <  7:263345/143‑1 (MQ=255)
                 ccgccCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTgg                                                                                                                                 <  7:151340/149‑1 (MQ=255)
                 ccgccCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCCAGGCGATCCAGGCGGGGGTCcgc                                                                                                                                       >  8:263345/1‑143 (MQ=255)
                     ccAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTg                                                                                                                             <  6:213071/149‑1 (MQ=255)
                        ggTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCTACGCCCACCGCGAGGCGATCCTGTCGGTGGTCCGCGATTGGGTTGgcg                                                                                                                          <  5:267555/149‑1 (MQ=255)
                         gTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGgcgc                                                                                                                         <  1:292739/149‑1 (MQ=255)
                           cgacgaTGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGgcgcgc                                                                                                                       >  5:62787/1‑149 (MQ=255)
                           cgacgaTGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGgcgcgc                                                                                                                       >  3:289459/1‑149 (MQ=255)
                                       gCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGcc                                                                                                            <  2:113291/148‑1 (MQ=255)
                                        cAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCga                                                                                                          <  5:2642/149‑1 (MQ=255)
                                                 gcgGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTgg                                                                                                  <  5:328082/148‑1 (MQ=255)
                                                 gcgGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTgg                                                                                                  <  4:129959/148‑1 (MQ=255)
                                                   ggCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGt                                                                                               <  3:257878/149‑1 (MQ=255)
                                                     ccATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTg                                                                                             <  3:156421/149‑1 (MQ=255)
                                                       aTGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCt                                                                                           <  6:62787/149‑1 (MQ=255)
                                                         gCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCggtg                                                                                                                                             >  8:96889/1‑97 (MQ=255)
                                                         gCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCggtg                                                                                                                                             <  7:96889/97‑1 (MQ=255)
                                                           tCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGAt                                                                                       >  3:259890/1‑149 (MQ=255)
                                                           tCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACAACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGGTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGAt                                                                                       >  5:163034/1‑149 (MQ=255)
                                                             gAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATg                                                                                      <  2:256967/148‑1 (MQ=255)
                                                              aaGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATgg                                                                                     >  4:195121/1‑148 (MQ=255)
                                                                            ccgcccgGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGtt                                                                                              <  4:223998/125‑1 (MQ=255)
                                                                            ccgcccgGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGtt                                                                                              >  3:223998/1‑125 (MQ=255)
                                                                            ccgcccgGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAg                                                                      >  1:2831/1‑149 (MQ=255)
                                                                               cccgGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGGCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGc                                                                   >  3:222162/1‑149 (MQ=255)
                                                                                ccgGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGc                                                                   >  8:17190/1‑148 (MQ=255)
                                                                                               aCCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATc                                                   >  5:165517/1‑149 (MQ=255)
                                                                                                ccgccgCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGt                                                             <  1:113919/138‑1 (MQ=255)
                                                                                                ccgccgCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGt                                                             >  2:113919/1‑138 (MQ=255)
                                                                                                ccgccgCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCg                                                  <  4:222162/149‑1 (MQ=255)
                                                                                                      cTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATg                                               >  7:292528/1‑146 (MQ=255)
                                                                                                      cTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATg                                               <  8:292528/146‑1 (MQ=255)
                                                                                                         ccACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGt                                           <  2:2831/147‑1 (MQ=255)
                                                                                                                     cGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGc                             <  4:289459/149‑1 (MQ=255)
                                                                                                                            aaCGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATc                        >  3:12073/1‑147 (MQ=255)
                                                                                                                             aCGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAg                     <  7:17190/149‑1 (MQ=255)
                                                                                                                               gCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCa                                                                 <  4:253953/103‑1 (MQ=255)
                                                                                                                               gCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCa                                                                 >  3:253953/1‑103 (MQ=255)
                                                                                                                                      gcgAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCaa             >  6:10674/1‑148 (MQ=255)
                                                                                                                                       cgAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGGTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACa           <  8:113900/149‑1 (MQ=255)
                                                                                                                                            cGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACATGGTc      <  6:165517/149‑1 (MQ=255)
                                                                                                                                             tcgccTCGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACATgg        >  6:142535/4‑146 (MQ=255)
                                                                                                                                             tcgccTCGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACATgg        <  5:142535/143‑1 (MQ=255)
                                                                                                                                                ccTCGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACATGGTCATCg  >  5:325568/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGCGACCTGGTGATCGCCGCCCAGGTCGACGATGCCTGGCAGCGCTTCGCGGCCATGCTCGAAGCCAACGGCGGCCCGCCCGGGTTCGGCAGCGACCGCCGCTACCACGGCCTCAACGGCCGCAACGCCCACCGCGAGGCGATCCTGGCGGTGGTCCGCGATTGGGTTGGCGCGCGTCGGGTGGCCGAGATCCTGGCGTTGCTCGATGGCATCGACATCCCCAGTGCCAAGGTTCAACGCATCGATGAAGTGCTGGCCGACCCGCAGATCCAGGCACGCAACATGGTCATCG  >  NC_002947/3397028‑3397320

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