Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,489,621 C→A 16.1% A4135S (GCT→TCT)  ebh ← extracellular matrix binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,489,6210CA16.1% 88.0 / 7.3 38A4135S (GCT→TCT) ebhextracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base C (20/11);  new base A (3/3);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

TGACTTGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTTACTTGATTAGCTGCAGCTGTTACTTCTG  >  CP000730/1489482‑1489755
                                                                                                                                           |                                                                                                                                      
tGACTTGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCt                                                                                                                                       >  1:284101/1‑141 (MQ=255)
  tcTTGTTTTGACTGTTTTTTCGTTAGCAATACTATCTCTCAAGCCTTCCCTTGCATTGTTCAATGCTTGTCCATTTGTTTAAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCttt                                                                                                                                     <  2:359894/140‑1 (MQ=255)
   cTTGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCtttt                                                                                                                                    >  1:116275/1‑141 (MQ=255)
   cTTGTTTTGACTGTTGTTTCGTTAGCAATACTATATCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTGTTAATTTTGCgttt                                                                                                                                    >  1:106939/1‑141 (MQ=255)
    ttGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGAttttt                                                                                                                                   >  2:393259/1‑141 (MQ=255)
    ttGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGAttttt                                                                                                                                   >  2:204570/1‑141 (MQ=255)
    ttGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGAttttt                                                                                                                                   >  2:204544/1‑141 (MQ=255)
     tGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTTTCCCACTTGTTCTTTTATTTTAGCTTTTTg                                                                                                                                  >  2:207716/1‑141 (MQ=255)
     tGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTAATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTACTCTATTGGCTTGTCACACCTGTACTTTTATTTTAGGTTTTTg                                                                                                                                  >  2:41102/1‑141 (MQ=255)
                    ttccttatataTAGTATCTTCTATGCTTTTTATTGTATTGTTCAATGTTTGTCCATTTGTTTGAACAGTTTGTAAGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAAGTGTATTTAATtg                                                                                                                   <  2:67543/131‑1 (MQ=255)
                      cGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGCACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATtgtg                                                                                                                 <  1:288047/141‑1 (MQ=255)
                       gTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTgtgt                                                                                                                >  2:43010/1‑141 (MQ=255)
                       gTTAGCAATACTATCTCTTAAGCCTTTCATTGCAGTGTTCCATGCTTGTCCATCCGTTTTAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTgtgt                                                                                                                >  2:380426/1‑141 (MQ=255)
                           gCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAAt                                                                                                            <  2:307256/141‑1 (MQ=255)
                                aCTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGATTTTTGAGGTGTATTTAATTGTGTTAATGCATc                                                                                                       <  1:393259/141‑1 (MQ=255)
                                aCTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGATTTTTGAGGTGTATTTAATTGTGTTAATGCATc                                                                                                       <  1:204570/141‑1 (MQ=255)
                                aCTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGATTTTTGAGGTGCATTTAATTGTGTTAATGCATc                                                                                                       <  1:204544/141‑1 (MQ=255)
                                aCTATCCCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATc                                                                                                       >  2:358522/1‑141 (MQ=255)
                                 cTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCa                                                                                                      >  2:324699/1‑141 (MQ=255)
                                             ccTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTg                                                                                          >  2:151004/1‑141 (MQ=255)
                                              cTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTgg                                                                                         >  1:259254/1‑141 (MQ=255)
                                                      gcattgTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCgttttgt                                                                                 <  2:359002/141‑1 (MQ=255)
                                                         tttttCCATGCTTGTCCCTTTGTTTGATCAGTTTGTACGTCTTCGAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTCGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCgttttgtttt                                                                              <  2:227900/138‑1 (MQ=255)
                                                             tCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATCGGCGTTTTGTTTTGCtt                                                                          >  1:298550/1‑141 (MQ=255)
                                                                     tGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAAt                                                                  <  2:284101/141‑1 (MQ=255)
                                                                       tCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCt                                                                >  1:199067/1‑141 (MQ=255)
                                                                                   aaCAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCAtt                                                    <  2:302554/141‑1 (MQ=255)
                                                                                       gTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTTTTTTTTTTTTAATTATGTTAATGCATCAATAGCAGTCTTTGCGTTTTGTTTTGCTTCCCGTAATCTTTCGTCACCATTTTAt                                                >  2:153288/1‑141 (MQ=255)
                                                                                       gTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAAt                                                >  2:346685/1‑141 (MQ=255)
                                                                                       gTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGGATTTAATTGTGTTAAGGCATCAATAGCAGGATTGGCGTTGTGTTTGTCTTCTCGTGGTCTTTTGTAACCAtgtttt                                                >  1:378828/1‑136 (MQ=255)
                                                                                                    ttAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCgcgc                                   <  1:213301/140‑1 (MQ=255)
                                                                                                           aTTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTTAc                            >  1:147404/1‑141 (MQ=255)
                                                                                                                cTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTTACTTGAt                       >  1:66729/1‑141 (MQ=255)
                                                                                                                   gTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTCACTTGTTTAg                    >  2:368448/1‑141 (MQ=255)
                                                                                                                     cccACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTGATTCTTGTTTCGCGTTGTCTACTTGATTAGCt                  >  2:64894/1‑141 (MQ=255)
                                                                                                                     cccACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAAATATTGTTTCGCGCTGTTTACTTGTTTAGCt                  >  2:277049/1‑141 (MQ=255)
                                                                                                                                     aTTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTTACTTGATTAGCTGCAGCTGTTACTTCTg  <  2:332445/141‑1 (MQ=255)
                                                                                                                                     aTTTAGCTTTTTGAGGTGTATTTAATTGCGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTTACTTGATTAGCTGCAGCTGTTACTTCTg  <  1:274821/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
TGACTTGTTTTGACTGTTGTTTCGTTAGCAATACTATCTCTTAAGCCTTTCATTGCATTGTTCAATGCTTGTCCATTTGTTTGAACAGTTTGTACGTCTTCTAATCTATTGGCTTGTCCCACTTGTTCTTTTAATTTAGCTTTTTGAGGTGTATTTAATTGTGTTAATGCATCAATAGCAGTATTGGCGTTTTGTTTTGCTTCACGTAATCTTTCGTCACCATTTAATTCTTGTTTCGCGCTGTTTACTTGATTAGCTGCAGCTGTTACTTCTG  >  CP000730/1489482‑1489755

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: