Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,844,193 C→T 7.7% D163N (GAC→AAC)  PP_3393 ← CAIB/BAIF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,844,1930CT 7.7% 151.1 / 2.2 52D163N (GAC→AAC) PP_3393CAIB/BAIF family protein
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base T (2/2);  total (26/26)
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.23e-01

ACCACGTACTGCTGCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCTGAGAAACCCCTTGTGCGAAGCGTA  >  NC_002947/3844048‑3844331
                                                                                                                                                 |                                                                                                                                             
aCCACGTACTGCTGCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑Gt                                                                                                                                            <  8:14863/148‑1 (MQ=255)
  cACGTACTGCTGCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCa                                                                                                                                         >  4:157046/1‑149 (MQ=255)
  cACGTACTGCTGCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCa                                                                                                                                         >  2:191772/1‑149 (MQ=255)
        ctgctgCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAg                                                                                                                                    <  3:209103/148‑1 (MQ=255)
        ctgctgCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGc                                                                                                                                   <  2:131764/149‑1 (MQ=255)
           ctgcATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGctgc                                                                                                                                >  1:293969/1‑149 (MQ=255)
              cATATGCTGGGCAGCCAGCAGCCCGCAATTCTCGTACAGGGCACTTTGAACCTCACGGCCCCTGCCGCTGCTGTGGCGCGCGTTGAGTGCGGCCCGCACGCCAATGGCGCCGAACCTGACACCCATGATGTC‑GTTCACCGAGctgcctg                                                                                                                             >  3:255467/1‑149 (MQ=255)
                       ggCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGc                                                                                                                     <  7:265930/148‑1 (MQ=255)
                       ggCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGc                                                                                                                     <  3:157046/148‑1 (MQ=255)
                        gCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCgg                                                                                                                   <  7:278496/149‑1 (MQ=255)
                        gCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCgg                                                                                                                   <  3:55249/149‑1 (MQ=255)
                         cagccagcAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCgg‑‑g                                                                                                                <  1:302977/149‑1 (MQ=255)
                            ccagcAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑gcgc                                                                                                             >  3:111310/1‑149 (MQ=255)
                            ccagcAGCACGCAATTCTCGTACAGGGCACGTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GATCACCGAGCTGCCTGCGCGCAGCGG‑‑gc                                                                                                               <  8:165756/147‑1 (MQ=255)
                               gcagcaCGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGcccc                                                                                                          >  1:79461/1‑149 (MQ=255)
                                        aaTTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCt                                                                                                 <  1:191772/149‑1 (MQ=255)
                                           tctcGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTT‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGt                                                                                               >  4:50698/1‑148 (MQ=255)
                                           tctcGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTT‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGt                                                                                               >  8:314167/1‑148 (MQ=255)
                                           tctcGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGt                                                                                               >  3:122104/1‑148 (MQ=255)
                                            ctcGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTc                                                                                              <  2:142991/148‑1 (MQ=255)
                                               gTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCAt                                                                                            <  7:250721/147‑1 (MQ=255)
                                               gTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCAt                                                                                            <  4:111310/147‑1 (MQ=255)
                                                   aGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGc                                                                                      >  6:350692/1‑149 (MQ=255)
                                                    gggCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACCTGCCCCCCATGATGTC‑GTGCCCCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCCCCGGCGCGGTCATGTAGGcc                                                                                     >  3:255700/1‑149 (MQ=255)
                                                    gggCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCgg                                                                                                                   >  6:292089/1‑121 (MQ=255)
                                                    gggCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCgg                                                                                                                   <  5:292089/121‑1 (MQ=255)
                                                     ggCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGc                                                                                      >  4:330097/1‑147 (MQ=255)
                                                                 cTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCcatca                                                                        >  6:46679/1‑149 (MQ=255)
                                                                  tCACGGCCCCTGCCGGTGCTGTGGCGCACGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCGAGACCGGCcatcat                                                                       >  7:215259/1‑149 (MQ=255)
                                                                   cACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGGGTGGGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCcagc                                                                         >  1:115294/1‑146 (MQ=255)
                                                                    aCGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCcatcat                                                                       <  3:330097/147‑1 (MQ=255)
                                                                    aCGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCt                                                                     >  4:138663/1‑149 (MQ=255)
                                                                        ccccTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATAATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCCCCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTAc                                                                 >  5:199145/1‑149 (MQ=255)
                                                                               cGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCa                                                          >  3:122159/1‑149 (MQ=255)
                                                                                  tGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATcc                                                       <  2:293969/149‑1 (MQ=255)
                                                                                     tgtgGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCAAGACCGGCCATCATCTGTACCACTTCATCCAgc                                                    >  5:330760/1‑149 (MQ=255)
                                                                                           gcTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTCGGTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAg                                               >  1:267214/1‑149 (MQ=255)
                                                                                           gcTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGcc                                                                            >  1:216770/1‑119 (MQ=255)
                                                                                           gcTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGcc                                                                            <  2:216770/119‑1 (MQ=255)
                                                                                               gTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCg                                            >  3:227102/1‑147 (MQ=255)
                                                                                               gTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCg                                            <  4:227102/147‑1 (MQ=255)
                                                                                               gTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCg                                            <  2:79461/147‑1 (MQ=255)
                                                                                                   aGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTAt                                       <  6:330760/148‑1 (MQ=255)
                                                                                                     tGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCa                                     >  8:139640/1‑148 (MQ=255)
                                                                                                       cGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCatat                                  >  5:137133/1‑149 (MQ=255)
                                                                                                          ccAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGc                               >  6:184698/1‑149 (MQ=255)
                                                                                                              cACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGc                           <  8:253914/149‑1 (MQ=255)
                                                                                                               aCGCCAATGGCGCCGAACATGCCACCCATGATGTT‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCt                          <  3:50698/149‑1 (MQ=255)
                                                                                                               aCGCCAATGGCGCCGAACATGCCACCCATGATGTT‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCt                          <  7:314167/149‑1 (MQ=255)
                                                                                                                              aaCATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCTGAGAAACCCCTtgtg           <  4:122159/149‑1 (MQ=255)
                                                                                                                                    ccACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCTGAGAAACCCCTTGTGCGAAg      <  5:184698/148‑1 (MQ=255)
                                                                                                                                       cccATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCTGAGAAACCCCTTGTGCGAAGCGTa  <  4:122104/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ACCACGTACTGCTGCATATGCTGGGCAGCCAGCAGCACGCAATTCTCGTACAGGGCACTTTGCACCTCACGGCCCCTGCCGGTGCTGTGGCGCTCGTTGAGTGCGGCCAGCACGCCAATGGCGCCGAACATGCCACCCATGATGTC‑GTTCACCGAGCTGCCTGCGCGCAGCGG‑‑GCGCCCCACCGGCCCTGTCATGTAGGCCAGACCGGCCATCATCTGTACCACTTCATCCAGCGCCAGGCGGTTATCATATGGCCCGCTGAGAAACCCCTTGTGCGAAGCGTA  >  NC_002947/3844048‑3844331

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