Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,499,075 G→A 9.2% N983N (AAC→AAT)  ebh ← extracellular matrix binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,499,0750GA 9.2% 124.7 / 4.2 44N983N (AAC→AAT) ebhextracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (2/2);  total (22/22)
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.58e-01

TTCTCAGTATTTCCTAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTCTAT  >  CP000730/1498938‑1499208
                                                                                                                                         |                                                                                                                                     
ttCTCAGTATTTCCTAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTa                                                                                                                                    <  1:30411/141‑1 (MQ=255)
 tctcAGTATTTCCTAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTACTAGTGACGTTGTCGGTCGCTGACACTTTAAATGTAATATCACTGTTaa                                                                                                                                   >  1:35590/1‑141 (MQ=255)
 tctcAGTATTTCCCAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTaa                                                                                                                                   <  2:449563/141‑1 (MQ=255)
     aGTATTTCCTAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGc                                                                                                                               >  1:387113/1‑141 (MQ=255)
             cTAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACCAGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCa                                                                                                                       >  1:108158/1‑141 (MQ=255)
               aaTACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATa                                                                                                                     <  1:345422/141‑1 (MQ=255)
                    aaTCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATactcac                                                                                                                <  1:328616/141‑1 (MQ=255)
                       cGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTATTAAATGCCTGATTCATACTCActct                                                                                                             >  2:359325/1‑141 (MQ=255)
                       cGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTATTAAATGCCTGATTCATACTCActct                                                                                                             >  2:67289/1‑141 (MQ=255)
                        ggATGAGCATTTTCACTTATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTa                                                                                                            <  2:23531/141‑1 (MQ=255)
                         gATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTAc                                                                                                           >  1:451228/1‑141 (MQ=255)
                           tgcgCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCt                                                                                                         <  1:82423/138‑1 (MQ=255)
                           tGAGCATCTTCACTAATTTTACCTGCATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTActt                                                                                                         <  2:286538/141‑3 (MQ=255)
                             aGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGt                                                                                                       >  1:88424/1‑141 (MQ=255)
                             aGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACCCTACCCGt                                                                                                       >  1:116200/1‑141 (MQ=255)
                                  cTTCACTAATTTTCCCTTCATTAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAg                                                                                                  <  1:172571/141‑1 (MQ=255)
                                  cTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAg                                                                                                  <  2:355041/141‑1 (MQ=255)
                                  cTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAg                                                                                                  <  2:77498/141‑1 (MQ=255)
                                  cTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAg                                                                                                  <  2:355038/141‑1 (MQ=255)
                                               aCCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCatt                                                                                     >  2:120511/1‑141 (MQ=255)
                                                    cATGAATTGAAACATGTTTAGATTGAGTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAACCCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCt                                                                                >  2:276344/1‑141 (MQ=255)
                                                       gAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTATTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTAt                                                                             <  1:67289/141‑1 (MQ=255)
                                                       gAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTATTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTAt                                                                             <  1:359325/141‑1 (MQ=255)
                                                        aaTTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATc                                                                            <  1:347409/141‑1 (MQ=255)
                                                         aTTGAAACATGTGTAGATTGACTATCATTTGTCGTCTTATTGACACTGACTGGCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATcc                                                                           <  1:356964/141‑1 (MQ=255)
                                                           tGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGc                                                                         <  2:451228/141‑1 (MQ=255)
                                                             aaaCATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGGTAAATACTTGATGCATACACACTCTACCTGTAATAGATAATGAGCCCTTATTCTTATCCGCtt                                                                       >  2:443518/1‑141 (MQ=255)
                                                                  tGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGtt                                                                  >  1:378903/1‑141 (MQ=255)
                                                                        gATTGTCTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAAttt                                                            <  1:207584/141‑1 (MQ=255)
                                                                         aTTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTg                                                           <  1:377959/141‑1 (MQ=255)
                                                                                             ttagtGACATTGTCGGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAAt                                       <  2:174629/141‑1 (MQ=255)
                                                                                                   acgttGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTa                                 >  1:231616/4‑141 (MQ=255)
                                                                                                       tgtctgtcGCTGACACTTTCAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCAAACCATTCGGCAATGTCGTAATAg                             <  1:294344/141‑1 (MQ=255)
                                                                                                       tgtctgtcGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTACTAg                             >  1:38420/1‑141 (MQ=255)
                                                                                                        gtctgtcGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCACTGTCGTAATAGt                            >  2:252772/1‑141 (MQ=255)
                                                                                                         tctgtcGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTCATAGATAATGAGCCATGATTCTTATCCGCTTTTATTAAATTTGATGTCCAACCATTCGGCACAGTCGTAATGGtt                           >  2:149514/1‑141 (MQ=255)
                                                                                                                     acTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTtata               >  2:27372/1‑141 (MQ=255)
                                                                                                                              gTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTATAATAGCTAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATATCCAACCATTCGGCAATGTCTGAATAATTGTATTTTTTATACCGGAATGt      >  1:135446/1‑141 (MQ=255)
                                                                                                                               tAATATCACTGTTAAATGCCTGATTCACACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTc     <  1:374041/141‑1 (MQ=255)
                                                                                                                                aaTATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATTTCGTAATAGTTGTATTTTTTATACCAGAATGTCt    >  2:162096/1‑141 (MQ=255)
                                                                                                                                aaTATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTCt    >  1:399266/1‑141 (MQ=255)
                                                                                                                                aaTATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTCt    >  2:6852/1‑141 (MQ=255)
                                                                                                                                aaTATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTCt    >  2:9344/1‑141 (MQ=255)
                                                                                                                                  gatCACTTTTAAATGCCTGATTCATTCTCCCTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAATTTTGATGTCCACCCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTCTAt  <  2:284416/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TTCTCAGTATTTCCTAATACAATCGGATGAGCATCTTCACTAATTTTACCTACATGAATTGAAACATGTTTAGATTGACTATCATTTGTCGTATTATTGACATTGTCTGTCGCTGACACTTTAAATGTAATATCACTGTTAAATGCCTGATTCATACTCACTCTACCTGTAATAGATAATGAGCCATTATTCTTATCCGCTTTTGTTAAATTTGATGTCCAACCATTCGGCAATGTCGTAATAGTTGTATTTTTTATACCAGAATGTCTAT  >  CP000730/1498938‑1499208

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