Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793613,0610CA75.4% ‑1.5 / 30.6 57Q639K (CAA→AAA) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base C (13/1);  new base A (43/0);  total (56/1)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CAAAATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTG  >  NC_007793/612963‑613139
                                                                                                  |                                                                                
ctaaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTGCa                                                                                 >  1:5326663/3‑100 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTGaaag                                                                               >  1:5554385/1‑99 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTGaaag                                                                               >  1:4018944/1‑99 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTACaag                                                                               >  1:3469556/1‑99 (MQ=9)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTACaag                                                                               >  1:3011002/1‑99 (MQ=9)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTACaag                                                                               >  1:2109447/1‑99 (MQ=9)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTAGCaag                                                                               >  1:1620141/1‑99 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTCTGcgca                                                                               >  1:3312767/1‑97 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTCTGaaag                                                                               >  1:2276454/1‑99 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTATGaaag                                                                               >  1:3643047/1‑99 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTATGaaag                                                                               >  1:2422664/1‑99 (MQ=255)
  aaaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACATTGaaag                                                                               >  1:3402479/1‑99 (MQ=255)
   aaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTAC‑TTGCaaaga                                                                             >  1:6423925/1‑98 (MQ=255)
    gtacagcttAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTGaaagat                                                                             >  1:1469103/8‑97 (MQ=9)
    gtacagcttAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTATGaaagat                                                                             >  1:2616836/8‑97 (MQ=255)
    ctacagcttAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTACaagag                                                                             >  1:1953243/8‑97 (MQ=9)
     tataGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACg                                                                            >  1:3663601/1‑100 (MQ=255)
     tacagcttAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTGCaagatc                                                                            >  1:4993251/7‑96 (MQ=255)
      acagcttAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACATTACaagatcg                                                                           >  1:5481942/6‑95 (MQ=255)
            tAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTaaaaaatagatcgga                                                                     >  1:5151356/1‑91 (MQ=255)
                                  cAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGaa                                             >  1:2905473/1‑100 (MQ=255)
                                               ggTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATg                                >  1:3533824/1‑100 (MQ=255)
                                                            gtgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCtt                   >  1:227761/5‑100 (MQ=255)
                                                            gctatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCtt                   >  1:27196/5‑100 (MQ=255)
                                                            atgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCtt                   >  1:3446005/5‑100 (MQ=255)
                                                            atgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCtt                   >  1:681845/5‑100 (MQ=255)
                                                             tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:129650/4‑100 (MQ=2)
                                                             tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:3787447/4‑100 (MQ=2)
                                                             tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:1267738/4‑100 (MQ=2)
                                                             tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGGGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:5113410/4‑100 (MQ=2)
                                                             ctatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttc                  >  1:2022006/4‑99 (MQ=2)
                                                             cccTGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:5096094/3‑100 (MQ=2)
                                                             agatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:2961382/4‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:4728751/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:1727864/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:521347/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:2683047/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:3285052/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:4460517/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:441559/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:3982069/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:2980992/3‑100 (MQ=2)
                                                              gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:3633759/3‑100 (MQ=2)
                                                               cTGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGa                >  1:4824275/1‑100 (MQ=255)
                                                                tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt               >  1:289533/1‑100 (MQ=2)
                                                                tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt               >  1:1889856/1‑100 (MQ=2)
                                                                tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt               >  1:3555046/1‑100 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGCTGAAAATGGTAAGTATCGCTTTGAtg              >  1:3662364/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:3048495/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:4225794/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:6389444/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:5696475/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:1367093/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:5450884/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:2895497/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:498502/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:4656005/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:2829124/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTAACGCTTTGAtg              >  1:5546717/1‑99 (MQ=2)
                                                                 gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGG‑TAACACATT‑GAAAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:6053757/1‑99 (MQ=2)
                                                                 ccaaaaGGCATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:3029585/3‑99 (MQ=255)
                                                                  aaaaaGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtga             >  1:1298781/1‑98 (MQ=255)
                                                                   aaaaGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATATTGAAGG‑TAACAAATTAA‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtgac            >  1:4295906/1‑97 (MQ=2)
                                                                   aaaaGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtgac            >  1:5858651/1‑97 (MQ=2)
                                                                   aaaaGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtgac            >  1:4329480/1‑97 (MQ=2)
                                                                     aaGGAATTAAAGGTGTTAAAGTTTCTTTGCAAAACGAAAAAGGCGATGTAATT‑G‑GGACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAAttt          <  1:4687864/100‑1 (MQ=255)
                                                                             aaaGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTg  >  1:5408555/1‑100 (MQ=255)
                                                                                                  |                                                                                
CAAAATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTG  >  NC_007793/612963‑613139

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