Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 613,064 A→G N640D (AAC→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,0640AG89.6% 21.4 / 2.2 48N640D (AAC→GAC) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base G (43/0);  total (47/1)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

AATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTG  >  NC_007793/612966‑613139
                                                                                                  |                                                                             
aaTACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAAAGGTGTTAAAGTTAC‑TTGCaaaga                                                                             >  1:6423925/1‑98 (MQ=255)
  tataGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACg                                                                            >  1:3663601/1‑100 (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:681845/5‑100 (MQ=255)
                                                         atgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCtt                   >  1:3446005/5‑100 (MQ=255)
                                                          tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:3787447/4‑100 (MQ=2)
                                                          tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:1267738/4‑100 (MQ=2)
                                                          tgatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCttt                  >  1:129650/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:441559/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:521347/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:4728751/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:4460517/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:3982069/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:1727864/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:2683047/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:3633759/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:3285052/3‑100 (MQ=2)
                                                           gatGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTg                 >  1:2980992/3‑100 (MQ=2)
                                                            cTGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGa                >  1:4824275/1‑100 (MQ=255)
                                                             tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt               >  1:1889856/1‑100 (MQ=2)
                                                             tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt               >  1:3555046/1‑100 (MQ=2)
                                                             tGAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAt               >  1:289533/1‑100 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGCTGAAAATGGTAAGTATCGCTTTGAtg              >  1:3662364/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:2895497/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:5450884/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:2829124/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:1367093/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:498502/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:3048495/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:4656005/1‑99 (MQ=2)
                                                              gAAAAAGGTATTAAAGGTGTTAAAGTTACTTTAAAAGATAGTGAAGGTAACACATTA‑A‑AAACAACTGAAACTGATGAAAATGGTAAGTATCGCTTTGAtg              >  1:4225794/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)
                                                                                                  |                                                                             
AATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGAAAAAGGAATTAAAGGTGTTAAAGTTACTTTGCAAAACGAAAAAGGCGAAGTAATT‑G‑GTACAACTGAAACAGATGAAAATGGTAAATACCGCTTTGATAATTTAGATAGTG  >  NC_007793/612966‑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: ‑