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,3550GT96.7% 90.6 / ‑2.6 30A6890E (GCA→GAA) ebhextracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (14/15);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CATCGCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTCAGCAG  >  CP000730/1481219‑1481490
                                                                                                                                        |                                                                                                                                       
cATCGCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCAtt                                                                                                                                     <  1:210552/141‑1 (MQ=255)
    gCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCTTTGttt                                                                                                                                 <  2:366183/141‑1 (MQ=255)
     cTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGtttt                                                                                                                                <  2:303985/141‑1 (MQ=255)
               ttCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTATCATCTGt                                                                                                                      <  1:327225/141‑1 (MQ=255)
                   tGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTCCt                                                                                                                  >  1:9492/1‑141 (MQ=255)
                    gggCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACtt                                                                                                                 <  1:242789/139‑1 (MQ=255)
                      gCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTTATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACtttt                                                                                                               >  1:83254/1‑141 (MQ=255)
                            tGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCa                                                                                                         >  2:20485/1‑141 (MQ=255)
                            tGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCa                                                                                                         <  2:335872/141‑1 (MQ=255)
                                    ttGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCaa                                                                                                 >  2:31275/1‑141 (MQ=255)
                                      gCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACg                                                                                               <  1:27057/141‑1 (MQ=255)
                                           aCTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATCCAACGCTTGt                                                                                          >  1:300278/1‑141 (MQ=255)
                                                    aCAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAgtt                                                                                 <  2:111202/141‑1 (MQ=255)
                                                    aCAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAgtt                                                                                 <  1:325218/141‑1 (MQ=255)
                                                          tgCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACa                                                                           >  1:344095/1‑141 (MQ=255)
                                                           gcgcTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAg                                                                          <  2:121892/141‑1 (MQ=255)
                                                             gcTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGtt                                                                        <  1:31275/141‑1 (MQ=255)
                                                                     gatGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGtgc                                                                <  1:74027/141‑1 (MQ=255)
                                                                          ttCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCtt                                                           >  2:220177/1‑141 (MQ=255)
                                                                            ccAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCa                                                         >  1:110052/1‑141 (MQ=255)
                                                                                  aCAGCACCTTTTTGTGTTTTGTTTCAGTTATCTAACGACGCTAATGTACGGTATTAATATTTTTTAGAATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTAAAAGTGCTGCTTCAACTtgt                                                   >  2:166199/1‑141 (MQ=255)
                                                                                                 gCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCAttt                                    <  2:140455/141‑1 (MQ=255)
                                                                                                  cATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTg                                   >  1:26664/1‑141 (MQ=255)
                                                                                                          aaGTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATGGGTTCGCATTTGTACCAttt                           >  1:105409/1‑141 (MQ=255)
                                                                                                            gTTATCTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGc                         >  2:281522/1‑141 (MQ=255)
                                                                                                                 cTAACGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCATGAt                    >  1:144228/1‑141 (MQ=255)
                                                                                                                    aCGTAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATtaa                 <  1:222951/141‑1 (MQ=255)
                                                                                                                       tAGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATtaataa              <  2:324851/141‑1 (MQ=255)
                                                                                                                        aGATAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTaataat             <  1:271590/141‑1 (MQ=255)
                                                                                                                           tAATGTTTGGTTTTCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATAtt          >  2:350351/1‑141 (MQ=255)
                                                                                                                                   ggttttCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTcagcag  >  1:74033/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CATCGCTGTATTCAGTTCCTGAGCGGTTTGAATGGCTTGCGTTACTTCAGCTACAGTGTGCGCTTGATTGATGTTTCCAGTAACAGCACCTTTTTGTGCATTGTTTAAGTTATCTAACGTAGATAATGTTTGGTTTGCATTGTTTTTAGCATCTGTTACTTTTCTATCACCATTCAACGCTTGTTTCGTAGTTGTTACAGTTGAAAGTGCTGCTTCAACTTGTGATTGGTTCGCATTTGTACCATTTGCTTGATTAATAATATTTTCAGCAG  >  CP000730/1481219‑1481490

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