Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 616,780 G→C K808Y (AAG→TAC)  SAUSA300_RS02920 → serine‑aspartate repeat‑containing protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793616,7800GC86.4% 24.7 / 9.2 44K808N (AAG→AAC) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base C (38/0);  total (43/1)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AATGGTGCTGATAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTAAAACAGATAAAGATGGTAA  >  NC_007793/616682‑616878
                                                                                                     |                                                                                                  
caTGGTGCTGATAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGGTGGTAAga                                                                                                   >  1:2394225/2‑99 (MQ=255)
aaTGGTGCTGATAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATAC‑AATAAAGATCGGAAgag                                                                                                  >  1:2209997/1‑98 (MQ=21)
         gATAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAGATAAAGATCGGAAgagcacacgtc                                                                                          >  1:201334/1‑90 (MQ=18)
         gATAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATTAAGATCGGAAgagcacacgtc                                                                                          >  1:5120659/1‑90 (MQ=18)
                 gACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGAATTCAACTGaaaaa                                                                                  >  1:4827596/1‑100 (MQ=255)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5595585/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:458013/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:4724715/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:3017529/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5451988/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5568816/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5867380/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:598353/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:2174082/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:6087349/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:1379578/1‑100 (MQ=9)
                          tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGAGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5428455/1‑100 (MQ=9)
                          gtgTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5520636/3‑100 (MQ=9)
                          gaGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:6236295/2‑100 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAGATACAGCTTAGGTGATTATGTATGGGACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:1982653/1‑99 (MQ=2)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTATGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3317493/1‑99 (MQ=2)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2811448/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:790008/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:6413507/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:6352395/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:6302212/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2160550/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2165247/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:230942/1‑99 (MQ=2)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:5767901/1‑99 (MQ=2)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:5669408/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2394953/1‑99 (MQ=2)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2576655/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:280835/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3098494/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3146026/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:5089173/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:4761131/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3385712/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:4616051/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:4218874/1‑99 (MQ=2)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3889604/1‑99 (MQ=9)
                           aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3691568/1‑99 (MQ=9)
                                     aCAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGt                                                              >  1:3409248/1‑100 (MQ=255)
                                                                         tatgGGAAGATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTaca                             >  1:4858338/1‑100 (MQ=255)
                                                                                 gATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTa                     <  1:3410092/100‑1 (MQ=255)
                                                                                 gATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTa                     >  1:5147447/1‑100 (MQ=255)
                                                                                                    aGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTAAAACAGATAAAGAtcggaa  >  1:2181843/1‑95 (MQ=255)
                                                                                                     |                                                                                                  
AATGGTGCTGATAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTAAAACAGATAAAGATGGTAA  >  NC_007793/616682‑616878

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