Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,481,355 G→T A6890E (GCA→GAA)  ebh ← extracellular matrix binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,481,3550GT92.3% 71.6 / ‑3.8 26A6890E (GCA→GAA) ebhextracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (10/14);  total (11/15)
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.94e-01

ACCCATCGCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTCAGCAGC  >  CP000730/1481216‑1481491
                                                                                                                                           |                                                                                                                                        
aCCCATCGCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGc                                                                                                                                         >  2:82341/1‑141 (MQ=255)
         tGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGttttt                                                                                                                                <  1:284935/141‑1 (MQ=255)
          gTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTa                                                                                                                               >  1:223374/1‑141 (MQ=255)
             ttCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCa                                                                                                                            >  2:256511/1‑141 (MQ=255)
                aGTTCCTGAGCGGTTTGAATGGCTTGCGATACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCt                                                                                                                         <  1:258534/141‑1 (MQ=255)
                                              aCTTGAGCTACAGTGTGCGCTTGATTGATTTTTCGAGTAACGGCACTTTTTTGTGCATTATTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTGTTTGCCTCTGTTACTTCTCTATCACCATTAAACGCTTGt                                                                                           <  2:48730/141‑1 (MQ=255)
                                              aCTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGt                                                                                           <  2:31475/141‑1 (MQ=255)
                                                       aCAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAgtt                                                                                  <  2:178229/141‑1 (MQ=255)
                                                               cgcTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGt                                                                          <  2:11052/141‑1 (MQ=255)
                                                                  ttgattgaTGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGa                                                                       >  2:244289/1‑141 (MQ=255)
                                                                     attgatGTTTCCAGTAACAGCAACTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGATTGTTTAATAGTTGTTACAATTGAAAg                                                                    >  2:262593/1‑141 (MQ=255)
                                                                           ttttCCAGTAACAGCACCGTTTTGAACACTAGTTCACGAAGCTAACAAATATAATGTGTTGGTTGCTTTGGTTTGAGCGTCGTTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGtgctgc                                                              <  1:42152/140‑1 (MQ=255)
                                                                                cAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCaa                                                         >  1:1023/1‑141 (MQ=255)
                                                                                  gTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACt                                                       <  2:33892/141‑1 (MQ=255)
                                                                                  gTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACt                                                       <  2:33888/141‑1 (MQ=255)
                                                                                       aGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGa                                                  <  2:186773/141‑1 (MQ=255)
                                                                                                ttGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATTAGTTACTGTTCTATAACGATTGGACGCGGGTTTCTTTGTGGTTACGGTTGAAAGTGCTGTTTCAACTTGTGATTGGTTCGc                                         >  1:277572/1‑141 (MQ=255)
                                                                                                  tttttttttttAAGTTATATAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCACCGCTTTTTTTTTTTTTTTTACAGTTGTAAGTGCTGCTTCAACTTGTTATTGGttttcat                                       >  2:199933/9‑136 (MQ=255)
                                                                                                  gtgCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCAt                                       >  1:105261/1‑141 (MQ=255)
                                                                                                           ttAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCAt                              <  2:173440/141‑1 (MQ=255)
                                                                                                              aGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTg                           <  1:164908/141‑1 (MQ=255)
                                                                                                                ttATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCt                         <  2:151129/141‑1 (MQ=255)
                                                                                                                  aTCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTg                       >  1:66149/1‑141 (MQ=255)
                                                                                                                           aGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTaataat              <  2:67732/141‑1 (MQ=255)
                                                                                                                             aTAATGTTTTGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTGGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAatat            <  2:1023/141‑1 (MQ=255)
                                                                                                                                aTGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATAtttt         >  1:253716/1‑141 (MQ=255)
                                                                                                                                       gttttCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTcagcagc  >  2:233008/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                        
ACCCATCGCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTCAGCAGC  >  CP000730/1481216‑1481491

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