Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,680,753 G→A 9.3% L16L (CTC→CTT)  PP_1475 ← CAIB/BAIF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,680,7530GA 9.3% 117.3 / 2.7 43L16L (CTC→CTT) PP_1475CAIB/BAIF family protein
Reads supporting (aligned to +/- strand):  ref base G (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 = 9.27e-01

ACGAACCACCACAGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCGGGTACTGCTGCGGGGCATGGAAGGTGCGGATGG  >  NC_002947/1680610‑1680889
                                                                                                                                                |                                                                                                                                        
aCGAACAACCCCAGCTAGGGGCCATCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCCTCGGTCTATTCGACCTTGACCACTTCTTCTCCGAATTCGTCACAGATACGCTAGGCGAACGGCCCGTCGATCAGGTTACCGAGTTc                                                                                                                                     <  5:71895/149‑1 (MQ=255)
  gAACCACCACAGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGa                                                                                                                                   >  1:90998/1‑149 (MQ=255)
  gAACCACCACAGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGa                                                                                                                                   >  4:306359/1‑149 (MQ=255)
      caccacAGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGCGGGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACaa                                                                                                                                >  1:257390/1‑149 (MQ=255)
       accacAGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCAAGTTCGACAACt                                                                                                                              >  7:25641/1‑149 (MQ=255)
          acaGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTa                                                                                                                           <  6:341624/149‑1 (MQ=255)
           caGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTaa                                                                                                                          <  8:205034/149‑1 (MQ=255)
                      ccTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAgg                                                                                                               >  2:263613/1‑149 (MQ=255)
                       cTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTc                                                                                                                                     >  8:77280/1‑126 (MQ=255)
                       cTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAggg                                                                                                              <  2:257390/149‑1 (MQ=255)
                             cAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCAAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTg                                                                                                         <  8:25641/148‑1 (MQ=255)
                               gCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTg                                                                                                      <  1:263613/149‑1 (MQ=255)
                               gCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTg                                                                                                      <  3:306359/149‑1 (MQ=255)
                                cTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTgg                                                                                                     <  4:296378/149‑1 (MQ=255)
                                             cgcAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATgg                                                                                         <  2:101313/148‑1 (MQ=255)
                                             agcAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGTGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATgg                                                                                         <  8:281995/147‑1 (MQ=255)
                                                    gggTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTgg                                                                                                     >  2:65414/1‑130 (MQ=255)
                                                    gggTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTgg                                                                                                     <  1:65414/130‑1 (MQ=255)
                                                      gTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGt                                                                                <  4:299887/149‑1 (MQ=255)
                                                       tCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTg                                                                               >  8:135959/1‑149 (MQ=255)
                                                        cgccgcCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTgg                                                                              <  7:135959/149‑1 (MQ=255)
                                                                          cGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCTACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGt                                                             <  7:220596/148‑1 (MQ=255)
                                                                             ccTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGt                                                             <  3:318026/145‑1 (MQ=255)
                                                                             ccTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGt                                                             >  4:318026/1‑145 (MQ=255)
                                                                              cTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGttt                                                        >  3:47269/1‑149 (MQ=255)
                                                                                       ttCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTg                                                                        <  2:195317/124‑1 (MQ=255)
                                                                                         cGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGAGTCCGTGGCATTTGGGCAGTGCCGTGGGTTTATCACCCGtg                                             >  1:290183/1‑148 (MQ=255)
                                                                                                gAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCg                                      >  6:55582/1‑149 (MQ=255)
                                                                                                   ttCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGCGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCg                                   <  2:90998/149‑1 (MQ=255)
                                                                                                      ggCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCggg                                 >  6:331097/1‑148 (MQ=255)
                                                                                                        cacaGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGtt                                                         <  2:58112/122‑1 (MQ=255)
                                                                                                        cacaGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGtt                                                         >  1:58112/1‑122 (MQ=255)
                                                                                                                 gcgAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATggg                                                                                        >  6:44560/1‑82 (MQ=255)
                                                                                                                   gAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCGGGTACTGCTGCgggg                    >  7:259441/1‑148 (MQ=255)
                                                                                                                   gAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCAAACGTCGGTTTCGCCGGGTACTGCTGCgggg                    >  6:88592/1‑148 (MQ=255)
                                                                                                                       ctaccGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACAGGTCGGTTTCGCCGGGTACTGCTGCGGGGCATg                <  7:222693/144‑1 (MQ=255)
                                                                                                                          aCGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATAACACGTGGGTTTCGCCGGGTACTGCTGCGGGGAATg                >  8:222693/1‑145 (MQ=255)
                                                                                                                            ggCCCGGCGATCAGGGTACCAAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCGGGTACTGCTGCGGGGCATGGAAGGt          >  8:323117/1‑149 (MQ=255)
                                                                                                                              cccGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTcc                                                                                  <  4:76487/75‑1 (MQ=255)
                                                                                                                              cccGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTca                                                                                  >  3:76487/1‑74 (MQ=255)
                                                                                                                                    gATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCGGGTACTGCTGCGGGGCATGGAAGGTGCGGATgg  <  5:331097/149‑1 (MQ=255)
                                                                                                                                    gATCAGGGTACCAAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCGGGTACTGCTGCGGGGCATGGAAGGTGCGGATgg  <  7:323117/149‑1 (MQ=255)
                                                                                                                                               cGAGTTCGACAACTTTAAGGCCGGCGAGGGCTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGGTGGGCAGAGCCGTGGTATTATCCCACGTCGGTTTCGCCGGGTACTGCTGCg                       >  8:326996/1‑117 (MQ=255)
                                                                                                                                                |                                                                                                                                        
ACGAACCACCACAGCGAGGTGCCCTCGTACAGCTTGCGCCACTTGCGCAGC‑GGGTCGCCGCCATCGGGCGATTCGACCTTGACCACTTCGGCGCCGAATTCGGCACAGATACGCGAGGCGAACGGCCCGGCGATCAGGGTACCGAGTTCGACAACTTTAAGGCCGGCGAGGGGTTTGCTGGGCGTCGACATGGGGTGTCCGTGGCAGTTGGGCAGAGCCGTGGTTTTATCACACGTCGGTTTCGCCGGGTACTGCTGCGGGGCATGGAAGGTGCGGATGG  >  NC_002947/1680610‑1680889

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