Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 613,508 T→G S788A (TCA→GCA)  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,5080TG88.0% 19.6 / 6.1 50S788A (TCA→GCA) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base T (5/1);  new base G (44/0);  total (49/1)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

CAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTC  >  NC_007793/613413‑613581
                                                                                               |                                                                         
cagatgcggaCAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGActcagacgcaga                                                                       >  1:4095994/1‑100 (MQ=11)
        gaCAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgatt                                                               >  1:1474529/1‑100 (MQ=11)
        gaCAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgatt                                                               >  1:2212082/1‑100 (MQ=11)
        cccAGCGACTCAGATTCAGACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgatt                                                               >  1:6619162/3‑100 (MQ=11)
              tactcagattcagACAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgattctgatg                                                         >  1:2471071/2‑99 (MQ=11)
                 tcagatgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCTGATGCGGACAGCGAGTCAGACGCAGATAGCGattcagattcag                                                      >  1:5869738/10‑100 (MQ=2)
                 acagattcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGAGAGCgattctgatgcgg                                                      >  1:3350997/2‑96 (MQ=2)
                  cagatgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCTGATGCGGACAGCGCTTCAGACGCAGATAGCGAttcagattcaga                                                     >  1:6053280/9‑100 (MQ=2)
                  cagatgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCTGATGCGGACAGCGCTTCAGACGCAGATAGCGAttcagattcaga                                                     >  1:3872915/9‑100 (MQ=2)
                    ggtgcggacagcgattcagacgcagacagcgattcagacgcaGATAGCGATTCAGATGCGGATAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                   >  1:3049587/7‑98 (MQ=11)
                    ggtgcggacAGCGATTCAGACGCAGATAGCGATTCAGACTCAGACAGCGATTCTGATGCGGACAGCGATTCACACGCAGATAGCGattcagatttagcca                                                   >  1:6573055/7‑97 (MQ=2)
                    gattcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGAGTCTGCTGCGGACAGCGACTCAGACGCAGATAGCgattctgatgcggaca                                                   >  1:4325731/4‑93 (MQ=2)
                    gattcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCgattctgatgcggaca                                                   >  1:5984590/4‑93 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:1848814/1‑100 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:1199722/1‑100 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:5369032/1‑100 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:2815356/1‑100 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:1650908/1‑100 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCg                                                >  1:4007986/1‑100 (MQ=11)
                       tcagacAGCGACTCAGACGCAGATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCAGACAGCg                                                >  1:2965739/1‑100 (MQ=11)
                       gcgaacAGCGATTCAGACGCGGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCg                                                >  1:970308/5‑100 (MQ=11)
                        cagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGACGCGGACAGCGACTCAGACGCAGATAGCGATTCAGATTCAGATAGCGa                                               >  1:5568457/1‑100 (MQ=2)
                        cagacAGCGATTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGa                                               >  1:823330/1‑100 (MQ=255)
                           acagcgattcagacgcagacagcgattcagacgcagATAGCGATTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTc                                            >  1:4372147/1‑100 (MQ=11)
                                gactcagacgcagatagcgactcagacgcagatagcgacTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcg                                       >  1:1957731/4‑99 (MQ=2)
                                gactcagacgcagatagcgactcagacgcagatagcgacTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcg                                       >  1:257146/4‑99 (MQ=2)
                                ctctcagacgcagatagcgactcagacgcagatagcgacTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcg                                       >  1:4254332/4‑99 (MQ=2)
                                    cAGACGCAGATAGCGATTCAGACTCAGATAGCGACTCTGATGCGGACAGCGACTCAGACTCAGATAGCGATTCTGACGCGGACAGCGATTCTGAcgcgga                                   >  1:2868377/1‑95 (MQ=2)
                                     aGACCCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGatgcggac                                  >  1:4238419/1‑93 (MQ=2)
                                      tacgcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGGTTCAGACGCAGata                                 >  1:4804173/7‑98 (MQ=2)
                                      gacgcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGCTTCTGATGCGGACAGCGATTCAGACGCAGata                                 >  1:2847570/7‑98 (MQ=2)
                                       attcagacaGCGATTCAGACGCAGATAGCGATTCTGATGCGGCCAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGATTCTAACGCGGACAg                                >  1:6291264/9‑100 (MQ=2)
                                         gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:410507/4‑98 (MQ=2)
                                         gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:4458054/4‑98 (MQ=2)
                                         gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:1914692/4‑98 (MQ=2)
                                         gcagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:1938854/4‑98 (MQ=2)
                                         ccagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTCATGCGGACAGCGTTTCAGACGCAGATAgcg                              >  1:3809547/4‑98 (MQ=2)
                                         acagATAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAgcg                              >  1:4023369/4‑98 (MQ=2)
                                          cagATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTga                             >  1:5469767/3‑100 (MQ=255)
                                             aTAGCGACTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCGGATTCAGatcggaagag                          >  1:1805449/1‑91 (MQ=11)
                                              accagatgcagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattca                         >  1:3078718/9‑100 (MQ=0)
                                                 tctggcagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagac                      >  1:5223022/6‑99 (MQ=2)
                                                   actcagacgcagATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGCgattcagattc                    >  1:5905857/1‑100 (MQ=11)
                                                    tgcagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGGTAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgca                   >  1:4572881/3‑96 (MQ=2)
                                                    gacagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACGGCgattcagacgca                   >  1:4823954/3‑96 (MQ=2)
                                                    cccagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgca                   >  1:5518651/3‑96 (MQ=2)
                                                     ccagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgcag                  >  1:5598211/2‑95 (MQ=2)
                                                      cagacgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTCTGATGCGGACAGCgattcagacgcgga                 >  1:1724519/1‑94 (MQ=2)
                                                      cagacgcagaTAGCGATTCAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATCCAGACAGCGATTCTGACGCAGACAGTGAttcagacccaga                 >  1:4509746/1‑100 (MQ=11)
                                                       cgatgcagaTAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTAAGACAGCAAATCAGACGCAGATAGCGCtgcagaagcagat                >  1:4583025/6‑99 (MQ=0)
                                                                     aCTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTc  <  1:2252624/100‑1 (MQ=255)
                                                                                               |                                                                         
CAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTC  >  NC_007793/613413‑613581

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