Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793613,6870CT72.1% 32.7 / 16.6 43D847D (GAC→GAT) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base C (12/0);  new base T (31/0);  total (43/0)
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.65e-01

ACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGA  >  NC_007793/613602‑613770
                                                                                     |                                                                                   
actcagacAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcagatagcgactcagactcagatag                                                                       >  1:998647/1‑100 (MQ=255)
   cagacAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcagatagcgactcagactcagatagtga                                                                    >  1:1761334/1‑100 (MQ=255)
   cagacAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcagatagcgactcagactcagatagcga                                                                    >  1:2300113/1‑100 (MQ=255)
               cagactcagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGAtagcgactcagactcagatagcgactcagactcagatagcgactcagactcaga                                                        >  1:2704982/1‑100 (MQ=255)
                    tcagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGCGATTCAGATTGAGACAGCGA‑TCAGA‑TCGGaagagcac                                                 >  1:116260/1‑93 (MQ=14)
                    gcagacagcgactcagactcagacagcgactcagactCAGATAGCGATTCAGATTCAGATAGCGATTCAGATTCAGACAGAGACTCAGA‑TCGGaagagca                                                  >  1:166508/2‑94 (MQ=11)
                    gcagacAGTGACTCAGGCTCAGACAGCGACTCAGGTTCAGATAGCGACTCAGAATCAGATAGCGATTCGGATTCAGACAGGGATTCAGA‑TCGGaagagca                                                  >  1:4614617/2‑94 (MQ=18)
                    ccagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGTGACTCAGA‑TCGGaagagca                                                  >  1:1471719/2‑94 (MQ=11)
                    ccagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGATAGCGATTCAGACTCAGACAGGGACTCAGA‑TCGGaagagca                                                  >  1:4491107/2‑94 (MQ=14)
                    ccagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGAAAGTGATTCAGATTCAGACAGTGACTCAGA‑TCGGaaggagc                                                  >  1:3875303/2‑94 (MQ=14)
                    acagacAGCGATTCAGATTCAGACAGCGACCCAGACCCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGCGATTCAGA‑TCGGaagagca                                                  >  1:926586/2‑94 (MQ=17)
                     ccgacAGTGATTCCGACTCAGACAGCGATTCAGATTCCGACAGTGATTCCGACTCAGATAGCGATTCCGACTCAGATAGCGACTCAGA‑TCGGaagagcac                                                 >  1:1776018/3‑93 (MQ=21)
                     cagatAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGTGACTCAGATTCAGATAGCGAATCAGA‑TCGGaagagcac                                                 >  1:168234/1‑93 (MQ=14)
                     cagatAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGTGACTCAGATTCAGATAGAGACTCAGA‑TCGGaagagcac                                                 >  1:5009109/1‑93 (MQ=14)
                     cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGTGACTCAGA‑TCGGaagagcac                                                 >  1:1177121/1‑93 (MQ=11)
                     cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGCTTCAGACTCAGATAGCCACTCAGA‑TCGGaagagcac                                                 >  1:852263/1‑93 (MQ=14)
                     cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGTGAATCAGA‑TCGGaagagcac                                                 >  1:1537331/1‑93 (MQ=14)
                     cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcac                                                 >  1:564166/1‑93 (MQ=14)
                     cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcac                                                 >  1:2207885/1‑93 (MQ=14)
                     cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACGCAGA‑TCGGaagagcac                                                 >  1:715417/1‑93 (MQ=14)
                     cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACGCAGA‑TCGGaagagcac                                                 >  1:2464271/1‑93 (MQ=14)
                     cagacAGTGACTCGGACTCAGACAGTGACTCAGACTCAGATAGCGGTTCAGACTCAGATAGCGATTCAGACTCAGACAGCGACACAGA‑TCGGaagagcac                                                 >  1:4966753/1‑93 (MQ=14)
                     cagacAGCGATTCAGATTCCGACAGTGATTCCGACTCAGATAGCGATTCCGACTCAGATAGCGACTCAGATTCAGACAGCAATTCAGA‑TCGGaagagcac                                                 >  1:2108500/1‑93 (MQ=18)
                     cagacAGCGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGA‑TCGGaagagcac                                                 >  1:4017470/1‑93 (MQ=17)
                     cagacAGCGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGCGATTCAGA‑TCGGaagagcac                                                 >  1:1625662/1‑93 (MQ=17)
                     cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGTAACTCAGA‑TCGGaagagcac                                                 >  1:1449300/1‑93 (MQ=14)
                     cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGAGACTCAGA‑TCGGaagagcac                                                 >  1:3368889/1‑93 (MQ=14)
                     cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                                                 >  1:84879/1‑93 (MQ=11)
                     cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCAACTCAGA‑TCGGaagagcac                                                 >  1:1938395/1‑93 (MQ=14)
                     cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGAGACTCAGA‑TCGGaagagcac                                                 >  1:212089/1‑93 (MQ=14)
                     cagacAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                                                 >  1:1214158/1‑93 (MQ=14)
                     cagacAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                                                 >  1:4663812/1‑93 (MQ=14)
                     cagacAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACCCAGA‑TCGGaagagcac                                                 >  1:4394301/1‑93 (MQ=14)
                     cacacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcac                                                 >  1:4271003/4‑93 (MQ=14)
                     atgacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCAACTCAGA‑TCGGaagagcac                                                 >  1:1064533/3‑93 (MQ=14)
                      tgacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCGACGCAGA‑TCGGaagagcaca                                                >  1:1560157/2‑92 (MQ=14)
                      cgacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGATTCAGACGGCGACTCAGACTCAGACAGTGACCCAGA‑TCGGaagagcaca                                                >  1:2073823/2‑92 (MQ=11)
                      agatAGCGACTCAGATTCAGACAGCGACTCAGATTCAGACAGCGACTCAGACTCAGATAGCGATTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcaca                                                >  1:3260847/1‑92 (MQ=14)
                      agacAGTGACTCAGGCTCAGACAGCGACTCAGGTTCAGATAGCGACTCAGAATCAGATAGCGATTCGGATTCAGACAGTGGTTCAGA‑TCGGaagagcaca                                                >  1:1868920/1‑92 (MQ=18)
                      agacAGCGACTCAGATTCAGACAGCGACTCAGGCTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcaca                                                >  1:2383958/1‑92 (MQ=11)
                      agacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCTTCTCAGA‑TCGGaagagcaca                                                >  1:5336667/1‑92 (MQ=14)
                      agacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGCGATTCTGATTCAGACAGCGACCCAGA‑TCGGaagagcaca                                                >  1:4824815/1‑92 (MQ=14)
                       gacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCGACTTAGA‑TCGGaagagcacac                                               >  1:574449/1‑91 (MQ=14)
                       gacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGCGATTCAGATTCAGACAGCAACTCAGA‑TCGGaagagcacac                                               >  1:5537475/1‑91 (MQ=14)
                        acAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcacacg                                              >  1:391343/1‑90 (MQ=14)
                        acAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCGACGCAGA‑TCGGaagagcacacg                                              >  1:2227259/1‑90 (MQ=14)
                                       cagaCAGCGATTCAGACTCAGATAGCGACTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGa                                >  1:5197306/1‑100 (MQ=255)
                                                                     cagactcagatagcgactcagactcagatagcgactcagactcGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGa  >  1:2601664/1‑100 (MQ=255)
                                                                                     |                                                                                   
ACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGA  >  NC_007793/613602‑613770

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