Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 616,786 6 bp→CAATGA coding (2430‑2435/4146 nt) 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,7860TC86.4% 28.1 / 4.8 44D810D (GAT→GAC) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base T (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 = 9.93e-01
*NC_007793616,7870TA86.4% 25.6 / 7.6 44S811T (TCA→ACA) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base T (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 = 9.37e-01
*NC_007793616,7880CA86.4% 26.6 / 7.3 44S811* (TCA→TAA) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base C (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 = 1.00e+00
*NC_007793616,7890AT86.4% 28.1 / 5.5 44S811S (TCA→TCT) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base T (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 = 1.00e+00
*NC_007793616,7900AG86.4% 28.4 / 4.9 44T812A (ACT→GCT) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base G (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 = 1.00e+00
*NC_007793616,7910CA86.4% 26.6 / 7.3 44T812N (ACT→AAT) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base C (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 = 9.08e-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: ‑