Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,955,443 C→A W287L (TGG→TTG)  SAUSA300_RS09680 ← leucotoxin LukDv

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,955,4430CA90.0% 20.7 / ‑4.9 10W287L (TGG→TTG) SAUSA300_RS09680leucotoxin LukDv
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (0/9);  total (0/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGTACCGATTAATTTAACAGTATGGTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCCAGTGTAGTCGATTCCATTGATTAGTATATCTATCCATTTCTCTTTGGTAAGTTACTTTGATTTT  >  NC_007793/1955349‑1955507
                                                                                              |                                                                
cGTACCGATTAATTTAACAGTATGGTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTTCCAACCCAgtgt                                                             <  1:3787908/100‑1 (MQ=255)
   aCCGATTAATTTAACAGTATGGTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGt                                                          <  1:2260029/100‑1 (MQ=255)
   aCCGATTAATTTAACAGTATGGTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGt                                                          <  1:4545737/100‑1 (MQ=255)
     cGATTAATTTAACAGTATGGTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGTca                                                        <  1:4362869/100‑2 (MQ=255)
                        gTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGTCGATTCCATTGAtgtgtgcta                                     <  1:4891743/100‑9 (MQ=255)
                                            aaGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGTCGATTCCATTGATTAGTATATCTATCCATTTCTCTTTGGTa                 <  1:4566538/100‑1 (MQ=255)
                                              gTAGATGTAATCGTTACTGTATTTTGATTTTTGTAGTTATGACCAACCAAGTTTAGTCGATTCCATTGATTAGTATATCTATCCATTtctcgggcatact               <  1:1440341/100‑10 (MQ=255)
                                               tAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGTCGATTCCATTGATTAGTATATCTATCCATTTCTCTTTGGTAAgg              <  1:1831869/100‑2 (MQ=255)
                                                           ttACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGTCGATTCCATTGATTAGTATATCTATCCATTTCTCTTTGGTAAGTTACTTTGAtttt  <  1:2153261/100‑1 (MQ=255)
                                                           ttACTGTATTTTGATTTTTGTAGTTATTACCAACCAAGTGTAGTCGATTCCATTGATTAGTATATCTATCCATTTCTCTTTGGTAAGTTACTTTGAtttt  <  1:5163548/100‑1 (MQ=255)
                                                                                              |                                                                
CGTACCGATTAATTTAACAGTATGGTTTTGCCAGTCAACTTCATAAGTAGATGTAAACGTTACTGTATTTTGATTTTTGTAGTTATTACCAACCCAGTGTAGTCGATTCCATTGATTAGTATATCTATCCATTTCTCTTTGGTAAGTTACTTTGATTTT  >  NC_007793/1955349‑1955507

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