Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 613,114 A→G K656K (AAA→AAG)  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,1140AG86.0% 16.6 / 4.3 50K656K (AAA→AAG) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base G (43/0);  total (48/2)
Fisher's exact test for biased strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACCTGCTGGCTTAACTCAAACAGGTACAAATACAACTGA  >  NC_007793/613023‑613203
                                                                                             |                                                                                         
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:1267738/4‑100 (MQ=2)
 tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                                                                                  >  1:129650/4‑100 (MQ=2)
 tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                                                                                  >  1:3787447/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:3633759/3‑100 (MQ=2)
  gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                                                                                 >  1:3982069/3‑100 (MQ=2)
  gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                                                                                 >  1:4728751/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:2980992/3‑100 (MQ=2)
  gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                                                                                 >  1:2683047/3‑100 (MQ=2)
  gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                                                                                 >  1:1727864/3‑100 (MQ=2)
  gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                                                                                 >  1:521347/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:3555046/1‑100 (MQ=2)
    tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt                                                                               >  1:1889856/1‑100 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGCTGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:3662364/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:498502/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:5450884/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:4656005/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:5696475/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:6389444/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:4225794/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:3048495/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:2895497/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:2829124/1‑99 (MQ=2)
     gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg                                                                              >  1:1367093/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:4329480/1‑97 (MQ=2)
       aaaaGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtgac                                                                            >  1:5858651/1‑97 (MQ=2)
         aaGGAATTAAAGGTGTTAAAGTTTCTTTGCAAAACGAAAAAGGCGATGTAATT‑G‑GGACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAAttt                                                                          <  1:4687864/100‑1 (MQ=255)
                 aaaGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTg                                                                  >  1:5408555/1‑100 (MQ=255)
                                                 gtgCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACctgct                                  >  1:109980/3‑100 (MQ=255)
                                                                    aaaCTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACCTGCTGGCTTAACTCAAACAgg                 >  1:2190064/2‑100 (MQ=255)
                                                                         gAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACCTGCTGGCTTAACTCAAACAGGTACaa            >  1:724469/1‑100 (MQ=255)
                                                                             cAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACCTGCTGGCTTAACTCAAACAGGTACAAATAc        >  1:5203193/1‑100 (MQ=255)
                                                                                   aaaaTGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACCTGCTGGCTTAACTCAAACAGGTACAAATACAACTGa  <  1:1195482/100‑1 (MQ=255)
                                                                                             |                                                                                         
CGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTGGTAAATACAAAGTTATCTTTGAAAAACCTGCTGGCTTAACTCAAACAGGTACAAATACAACTGA  >  NC_007793/613023‑613203

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: ‑