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,3550GT90.0% 96.4 / 5.8 40A6890E (GCA→GAA) ebhextracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (18/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTCAGCA  >  CP000730/1481226‑1481489
                                                                                                                                 |                                                                                                                                      
gTATTCAGTTCCTGAGCGGTTTGAATGGCCTGCGTTACTTCAGCTACCGTGTGCGCTTGATTGATGCTTCCCGTAACAGCACCTTTTTGTGCCTTGTTTAAGCTATCTAACGTAGATGATTTTTGGTTTTCATTGTTTTTa                                                                                                                             >  2:334269/1‑141 (MQ=255)
    tCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCAt                                                                                                                         >  2:115308/1‑141 (MQ=255)
       gTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGCTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTGATTGTTTTTGGCATCTg                                                                                                                      >  1:18843/1‑141 (MQ=255)
              aGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACttt                                                                                                               >  1:363294/1‑141 (MQ=255)
               gCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACtttt                                                                                                              <  2:94168/141‑1 (MQ=255)
                     tGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCa                                                                                                        <  1:9437/141‑1 (MQ=255)
                         tGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCAt                                                                                                    <  2:263648/141‑1 (MQ=255)
                             ttGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCaa                                                                                                <  1:302418/141‑1 (MQ=255)
                                cGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGc                                                                                             <  2:238541/141‑1 (MQ=255)
                                      ttCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGttt                                                                                       >  2:180550/1‑141 (MQ=255)
                                        cAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCACCGCTTGTTTCg                                                                                     <  2:155189/141‑1 (MQ=255)
                                        cAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCg                                                                                     <  1:194812/141‑1 (MQ=255)
                                        cAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCg                                                                                     <  2:317116/141‑1 (MQ=255)
                                        cAGCTACAGAGTGCGCTTGATTGATGATTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCg                                                                                     >  2:24723/1‑141 (MQ=255)
                                         aGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGt                                                                                    <  1:220440/141‑1 (MQ=255)
                                             aCAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAgtt                                                                                <  2:130135/141‑1 (MQ=255)
                                               aGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAgttgt                                                                              >  2:302466/1‑141 (MQ=255)
                                                         tgattgatGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGaa                                                                    <  1:339322/141‑1 (MQ=255)
                                                         tgattgatGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGaa                                                                    <  1:51748/141‑1 (MQ=255)
                                                         tgattgatGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGaa                                                                    >  1:208304/1‑141 (MQ=255)
                                                            ttgatGTTTCCAGTAACAGCACCTTTTTTTTTTTTTTTTAAGTTATCTAACCTAGATAATGTTTGTTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGt                                                                 >  2:141674/1‑141 (MQ=255)
                                                            ttgatGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGt                                                                 >  2:80750/1‑141 (MQ=255)
                                                                    tCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTc                                                         >  2:20403/1‑141 (MQ=255)
                                                                     ccAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCa                                                        <  2:202159/141‑1 (MQ=255)
                                                                     ccAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCa                                                        >  1:202159/1‑141 (MQ=255)
                                                                           aCAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTtgt                                                  >  1:295968/1‑141 (MQ=255)
                                                                                     tttGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCg                                        <  1:253872/141‑1 (MQ=255)
                                                                                      ttGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGc                                       <  1:330129/141‑1 (MQ=255)
                                                                                        gtgCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGAGTCTGTGACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTAAAAGCGCTGCTTCAACTTGTGATTTGTTCGCAt                                     >  1:362551/1‑141 (MQ=255)
                                                                                        gtgCATTGCTTAAGTTATCTAACGGAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTCGTTTCGTAGTTGTTACAGCTGAAAGTGCTGCTTCAACTTGTGATCGGTTTGCAt                                     >  2:242227/1‑141 (MQ=255)
                                                                                         tgCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAAc                                                      <  1:171775/123‑1 (MQ=255)
                                                                                         tgCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAAc                                                      >  2:171775/1‑123 (MQ=255)
                                                                                           cATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTg                                  >  1:105862/1‑141 (MQ=255)
                                                                                                tttAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCa                             <  1:368107/141‑1 (MQ=255)
                                                                                                   aaGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCAttt                          <  1:192376/141‑1 (MQ=255)
                                                                                                    aGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCACTTg                         >  2:323692/1‑141 (MQ=255)
                                                                                                    aGTTATCTAACGTAGATAATGTTTGGCCTTCATTGTTTTTAGCATGTGTTACTATTCTATCATCATTCTACGCTTGTTTTGTGGTTGTTAGAGTGGAGTGTGCGGCTTCAGCTTGTGATTGGGTGGCATTTGTagtgtttg                         >  2:228485/1‑134 (MQ=255)
                                                                                                      ttATCTAACGTAGATAATGTTTGGTTCTCATTGTTTTTAGCATCTGTTACTTTTCTCTCACCATCCAACGCTTGTTTCGTAGTTGTCACAGTTGACAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCt                       <  2:2048/141‑1 (MQ=255)
                                                                                                           tacgttAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGAtt                  >  1:267634/6‑141 (MQ=255)
                                                                                                             aCGTAGATAATGTTTGGTTTTCATTGTTTTTTTTTTTTTTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATtaa                >  2:98589/1‑141 (MQ=255)
                                                                                                                          ttggttttCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTcagc   >  2:172139/1‑141 (MQ=255)
                                                                                                                           tggttttCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTcagca  <  1:115308/141‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                      
GTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTCAGCA  >  CP000730/1481226‑1481489

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