Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 616,783 G→A Q809Q (CAG→CAA)  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,7830GA86.4% 28.7 / 4.3 44Q809Q (CAG→CAA) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base A (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.36e-01

GACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTAAAACAGATAAAGATGGTAA  >  NC_007793/616699‑616878
                                                                                       |                                                                                               
gACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGAATTCAACTGaaaaa                                                                                  >  1:4827596/1‑100 (MQ=255)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:4724715/1‑100 (MQ=9)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:1379578/1‑100 (MQ=9)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:2174082/1‑100 (MQ=9)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:458013/1‑100 (MQ=9)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:6087349/1‑100 (MQ=9)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:598353/1‑100 (MQ=9)
         tAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtt                                                                            >  1:5867380/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:5595585/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:5669408/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:5767901/1‑99 (MQ=2)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:6302212/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:5089173/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:6352395/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:6413507/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:790008/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3146026/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:2394953/1‑99 (MQ=2)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2576655/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:280835/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:2811448/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3098494/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:4761131/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3385712/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3691568/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:3889604/1‑99 (MQ=9)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:4218874/1‑99 (MQ=2)
          aGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTAtta                                                                           >  1:4616051/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)
                                                                                       |                                                                                               
GACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAATTTAGGTAATTATGTATGGGAAGATACAAATAAAGATGGTAAGCAGGATTCAACTGAAAAAGGTATTTCAGGCGTAACAGTTACATTGAAAAATGAAAACGGTGAAGTTTTACAAACAACTAAAACAGATAAAGATGGTAA  >  NC_007793/616699‑616878

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