Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 613,495 T→C D783D (GAT→GAC)  SAUSA300_RS02915 → serine‑aspartate repeat‑containing protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793613,4950TC89.7% 38.8 / 1.3 68D783D (GAT→GAC) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base T (6/1);  new base C (60/1);  total (66/2)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCAGATAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTC  >  NC_007793/613398‑613581
                                                                                                 |                                                                                      
tcgcagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTATGATGCGGACAg                                                                                      >  1:1339884/4‑100 (MQ=0)
ccgcagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAg                                                                                      >  1:4406746/4‑100 (MQ=2)
ccgcagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAg                                                                                      >  1:4439231/4‑100 (MQ=2)
ccgcagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAg                                                                                      >  1:1104961/4‑100 (MQ=2)
acccagaTAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATACCGATTCTGATGCGGACAg                                                                                      >  1:5419052/4‑100 (MQ=0)
  gcagaTAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                                                    >  1:1136642/2‑100 (MQ=11)
  acagaTAGCGATTCAGATGCGGACAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGAATCAGACGCAGATAGCGACTCTGATGCGGATAGCg                                                                                    >  1:229652/2‑100 (MQ=11)
   cagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCga                                                                                   >  1:903414/1‑100 (MQ=11)
   cagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCga                                                                                   >  1:493837/1‑100 (MQ=11)
   cagaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCga                                                                                   >  1:2617177/1‑100 (MQ=11)
   cagaTAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCga                                                                                   >  1:3310860/1‑100 (MQ=255)
   cagaTAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCga                                                                                   >  1:4260092/1‑100 (MQ=11)
   cagaTAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCga                                                                                   >  1:6130787/1‑100 (MQ=11)
   cagaTAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCGGATTCGGATAGCga                                                                                   >  1:4930532/1‑100 (MQ=255)
    agaTAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATGCGGACAGCgat                                                                                  >  1:596927/1‑99 (MQ=11)
    aaaTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCgac                                                                                  >  1:2460825/3‑100 (MQ=2)
      aTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCgactc                                                                                >  1:2850462/1‑100 (MQ=11)
      aTAGCGATTCAGATTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCgactc                                                                                <  1:5789987/100‑1 (MQ=11)
               cagatgcggaCAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGActcagacgcaga                                                                       >  1:4095994/1‑100 (MQ=11)
                       gaCAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgatt                                                               >  1:1474529/1‑100 (MQ=11)
                       gaCAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgatt                                                               >  1:2212082/1‑100 (MQ=11)
                       cccAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgatt                                                               >  1:6619162/3‑100 (MQ=11)
                             tactcagattcagACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgattctgatg                                                         >  1:2471071/2‑99 (MQ=11)
                                tcagatgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCTGATGCGGACAGCGAGTCAGACGCAGATAGCGattcagattcag                                                      >  1:5869738/10‑100 (MQ=2)
                                acagattcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGAGAGCgattctgatgcgg                                                      >  1:3350997/2‑96 (MQ=2)
                                 cagatgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCTGATGCGGACAGCGCTTCAGACGCAGATAGCGAttcagattcaga                                                     >  1:3872915/9‑100 (MQ=2)
                                 cagatgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCTGATGCGGACAGCGCTTCAGACGCAGATAGCGAttcagattcaga                                                     >  1:6053280/9‑100 (MQ=2)
                                   ggtgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCAGATGCGGATAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                   >  1:3049587/7‑98 (MQ=11)
                                   ggtgcggacAGCGATTCAGACGCAGATAGCGATTCAGACTCAGACAGCGATTCTGATGCGGACAGCGATTCACACGCAGATAGCGattcagatttagcca                                                   >  1:6573055/7‑97 (MQ=2)
                                   gattcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGAGTCTGCTGCGGACAGCGACTCAGACGCAGATAGCgattctgatgcggaca                                                   >  1:4325731/4‑93 (MQ=2)
                                   gattcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgattctgatgcggaca                                                   >  1:5984590/4‑93 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:1650908/1‑100 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:2815356/1‑100 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:4007986/1‑100 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:5369032/1‑100 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:1199722/1‑100 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:1848814/1‑100 (MQ=11)
                                      tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCAGACAGCg                                                >  1:2965739/1‑100 (MQ=11)
                                      gcgaacAGCGATTCAGACGCGGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCg                                                >  1:970308/5‑100 (MQ=11)
                                       cagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGACGCGGACAGCGACTCAGACGCAGATAGCGATTCAGATTCAGATAGCGa                                               >  1:5568457/1‑100 (MQ=2)
                                       cagacAGCGATTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGa                                               >  1:823330/1‑100 (MQ=255)
                                          acagcgattcagacgcagacagcgattcagacgcagATAGCGATTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTc                                            >  1:4372147/1‑100 (MQ=11)
                                               gactcagacgcagatagcgactcagacgcagatagcgacTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcg                                       >  1:1957731/4‑99 (MQ=2)
                                               gactcagacgcagatagcgactcagacgcagatagcgacTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcg                                       >  1:257146/4‑99 (MQ=2)
                                               ctctcagacgcagatagcgactcagacgcagatagcgacTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcg                                       >  1:4254332/4‑99 (MQ=2)
                                                   cAGACGCAGATAGCGATTCAGACTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACTCAGATAGCGATTCTGACGCGGACAGCGATTCTGAcgcgga                                   >  1:2868377/1‑95 (MQ=2)
                                                    aGACCCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGatgcggac                                  >  1:4238419/1‑93 (MQ=2)
                                                     tacgcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGGTTCAGACGCAGata                                 >  1:4804173/7‑98 (MQ=2)
                                                     gacgcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGCTTCTGATGCGGACAGCGATTCAGACGCAGata                                 >  1:2847570/7‑98 (MQ=2)
                                                      attcagacaGCGATTCAGACGCAGATAGCGATTCTGATGCGGCCAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGATTCTAACGCGGACAg                                >  1:6291264/9‑100 (MQ=2)
                                                        gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:4458054/4‑98 (MQ=2)
                                                        gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:1914692/4‑98 (MQ=2)
                                                        gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:410507/4‑98 (MQ=2)
                                                        gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:1938854/4‑98 (MQ=2)
                                                        ccagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTCATGCGGACAGCGTTTCAGACGCAGATAgcg                              >  1:3809547/4‑98 (MQ=2)
                                                        acagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:4023369/4‑98 (MQ=2)
                                                         cagATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTga                             >  1:5469767/3‑100 (MQ=255)
                                                            aTAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCGGATTCAGatcggaagag                          >  1:1805449/1‑91 (MQ=11)
                                                             accagatgcagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattca                         >  1:3078718/9‑100 (MQ=0)
                                                                tctggcagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagac                      >  1:5223022/6‑99 (MQ=2)
                                                                  actcagacgcagATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGCgattcagattc                    >  1:5905857/1‑100 (MQ=11)
                                                                   tgcagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGGTAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgca                   >  1:4572881/3‑96 (MQ=2)
                                                                   gacagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACGGCgattcagacgca                   >  1:4823954/3‑96 (MQ=2)
                                                                   cccagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgca                   >  1:5518651/3‑96 (MQ=2)
                                                                    ccagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgcag                  >  1:5598211/2‑95 (MQ=2)
                                                                     cagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgcgga                 >  1:1724519/1‑94 (MQ=2)
                                                                     cagacgcagaTAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATCCAGACAGCGATTCTGACGCAGACAGTGAttcagacccaga                 >  1:4509746/1‑100 (MQ=11)
                                                                      cgatgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTAAGACAGCAAATCAGACGCAGATAGCGCtgcagaagcagat                >  1:4583025/6‑99 (MQ=0)
                                                                                    aCTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTc  <  1:2252624/100‑1 (MQ=255)
                                                                                                 |                                                                                      
ACTCAGATAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTC  >  NC_007793/613398‑613581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑