Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793613,7110C.79.2% 21.5 / 40.2 48coding (2565/2844 nt)SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base . (38/0);  total (44/4)
Fisher's exact test for biased strand distribution p-value = 1.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-02

CAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTCAGATTCAGATAGCGATTCAGATTCAGA  >  NC_007793/613617‑613806
                                                                                               |                                                                                               
cagattcagacGGCGATTCAGATTCAGATAGCGATTCAGATTCAGACAG‑‑‑‑‑‑AGACTCAGATAGTGATTCAGACTCAGATAGCGAC‑CCAGA‑TCGGaagagcac                                                                                     >  1:921988/1‑93 (MQ=21)
cagactcagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGCC‑CCAGA‑TCGGaa                                                                                           >  1:1170529/1‑99 (MQ=14)
cagactcagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACCCAGATAGCGAC‑TCAGACTCGGa                                                                                            >  1:5107481/1‑100 (MQ=255)
 aacaggggacAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGAC‑ACAGA‑TCGGaag                                                                                          >  1:4820596/8‑98 (MQ=14)
  gatggggacAGCGATTCAGACTCAGATAGCGACTCAGACTCGGATAGCGACTCGGATTCAGATAGCGACTCAGACTCAGATAGTGAC‑TCAGA‑TCGGaaga                                                                                         >  1:594072/7‑97 (MQ=17)
   gtgcagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGAC‑CCAGA‑TCGGaagag                                                                                        >  1:560348/4‑96 (MQ=14)
    ttcagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCG‑C‑TCAGA‑TCGGaagagca                                                                                      >  1:1952147/2‑94 (MQ=17)
    tccagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGGTTCAGACAGTAAC‑TCAGA‑TCGGaagagc                                                                                       >  1:2923625/3‑95 (MQ=14)
     tcagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑GCAGA‑TCGGaagagca                                                                                      >  1:4531161/1‑94 (MQ=14)
     ccagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGAC‑CCAGA‑TCGGaagagca                                                                                      >  1:3677918/2‑94 (MQ=14)
      cagatAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGATAGTGATTCAGACTCAGATAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:838350/1‑93 (MQ=14)
      cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGTGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:669687/1‑93 (MQ=14)
      cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:534330/1‑93 (MQ=14)
      cagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCga                                                                                      >  1:3140241/1‑100 (MQ=255)
      cagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑ACAGA‑TCGGaagagcac                                                                                     >  1:5209610/1‑93 (MQ=14)
      cagacAGTGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:4450070/1‑93 (MQ=17)
      cagacAGCGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGGGAT‑TCAGA‑TCGGaagagcac                                                                                     >  1:4585515/1‑93 (MQ=17)
      cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGCC‑TCAGA‑TCGGaagagcac                                                                                     >  1:2056250/1‑93 (MQ=14)
      cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:40851/1‑93 (MQ=14)
      cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:730098/1‑93 (MQ=14)
      cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:1543593/1‑93 (MQ=14)
      cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGAC‑TCAGA‑TCGGAAAgagca                                                                                     >  1:892903/1‑96 (MQ=14)
      cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCAAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:4041854/1‑93 (MQ=14)
      cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGATAGTGATTCAGACTCAGATAGCGACTCAGATTCAGACAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:799571/1‑93 (MQ=17)
      cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGTGAC‑CCAGA‑TCGGaagagcac                                                                                     >  1:73747/1‑93 (MQ=11)
      cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCTGACAGCTTC‑TCAGA‑TCGGaagagcac                                                                                     >  1:517053/1‑93 (MQ=14)
      cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGTGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:855829/1‑93 (MQ=14)
      cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGTGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:3724478/1‑93 (MQ=14)
      cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGGGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:5439130/1‑93 (MQ=14)
      cagacAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCTAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:1981648/1‑93 (MQ=17)
      cagacAGCGACTCAGACCCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGTGAT‑TCAGA‑TCGGaagagcac                                                                                     >  1:1876247/1‑93 (MQ=14)
      acgacAGCGACTCAGATTCAGACAGCGATTCAGATTCAGACTGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGAC‑TCAGA‑TCGGaagagcac                                                                                     >  1:3458851/3‑93 (MQ=11)
       tgacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACGTCAGA‑TCGGaagagcac                                                                                     >  1:2254376/2‑93 (MQ=14)
       tgacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGTAAC‑TCAGA‑TCGGaagagcaca                                                                                    >  1:934528/2‑92 (MQ=14)
       cgacagcgactcagactcagacagcgactcagactcaGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCCAC‑TCAGA‑TCGGaagagcaca                                                                                    >  1:213288/2‑92 (MQ=11)
       agatAGCGACTCAGAATCAGATAGTGAGTCAGATTCAGACAGTGACTCGGACTCAGACAGTGATTCAGACTCAGATAGCGAT‑TCAGA‑TCGGaagagcaca                                                                                    >  1:4187686/1‑92 (MQ=18)
       agacagcgactcagactcagacagcgactcagactcAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGAC‑TCAGA‑TCGGaagagcaca                                                                                    >  1:5570226/1‑92 (MQ=14)
       agacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGTGAC‑CCAGA‑TCGGaagagcaca                                                                                    >  1:4687681/1‑92 (MQ=11)
        gacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCgatt                                                                                    >  1:2157616/1‑100 (MQ=255)
        gacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGAC‑TCAGA‑TCGGaagagcacac                                                                                   >  1:2725752/1‑91 (MQ=14)
         acAGCGACTCAGATTCAGACAGCGACTCAGACTCAGATAGTGATTCAGACTCAGATAGCGACTCAGATTCAGACAGCGAC‑TCAGA‑TCGGaagagcacacg                                                                                  >  1:2134166/1‑90 (MQ=14)
                          gaCAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACt                                                                  >  1:5519356/1‑100 (MQ=255)
                              cagacagcgattcagacAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcaga                                                              >  1:2459595/9‑100 (MQ=14)
                                 aCTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattc                                                           <  1:6105831/100‑1 (MQ=255)
                                 aCTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattc                                                           <  1:2659707/100‑1 (MQ=255)
                                 aCTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattc                                                           <  1:1634959/100‑1 (MQ=255)
                                   tcagattcagatAGCGACTCAGACCCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattcag                                                         <  1:4119688/100‑1 (MQ=255)
                                    cagattcagatAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattcaga                                                        >  1:2558486/1‑100 (MQ=255)
                                           agatagcgactcagactcagatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGAt                                                 >  1:2785405/1‑100 (MQ=255)
                                              caacgactcagattcagatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCg                                              >  1:4224335/4‑100 (MQ=255)
                                                   ccgcagactcagatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTc                                         >  1:4062645/4‑100 (MQ=255)
                                                      cagactcagatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGa                                      >  1:4998491/1‑100 (MQ=255)
                                                            cagatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGTGATTCGGACTCAGACAACGa                                >  1:4951958/1‑100 (MQ=255)
                                                            cagatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGa                                >  1:4319915/1‑100 (MQ=255)
                                                             agatagcgactcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGaa                               <  1:2747161/100‑2 (MQ=255)
                                                                     actcagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAGCGACtcagattc                       <  1:2279188/100‑1 (MQ=255)
                                                                        cagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACtcagatcggaa                    >  1:4891880/1‑95 (MQ=255)
                                                                        cagactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACtcagatcggaa                    >  1:1355114/1‑95 (MQ=255)
                                                                         agactcagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACtcagattcagat                   <  1:3307798/100‑1 (MQ=255)
                                                                              cagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTCAGATTCAGATAGCga              >  1:2951075/1‑100 (MQ=255)
                                                                              cagatagcgac‑tcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTCAGATTCAGATAGCga              >  1:1476427/1‑100 (MQ=255)
                                                                                           cagcctcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTCAGATTCAGATAGCGAttcagattcaga  >  1:3568219/1‑100 (MQ=255)
                                                                                               |                                                                                               
CAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGAC‑TCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTCAGATTCAGATAGCGATTCAGATTCAGA  >  NC_007793/613617‑613806

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