Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 47,368 G→A D51D (GAC→GAT)  SAUSA300_RS00200 ← recombinase RecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00779347,3680GA100.0% 43.0 / NA 16D51D (GAC→GAT) SAUSA300_RS00200recombinase RecB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/4);  total (12/4)

TGTATACCATAACACAGTCGATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGGTCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGTATTCTTTGATTTGATTGATTTGTCCGTCGATACTATAACCTTCTGTGCTTTGC  >  NC_007793/47274‑47464
                                                                                              |                                                                                                
tgtAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTTCCTGAAATACCACGATCAGc                                                                                             <  1:5582081/100‑1 (MQ=255)
         tAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATc                                                                                    >  1:1992261/1‑100 (MQ=255)
         tAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATc                                                                                    >  1:869067/1‑100 (MQ=255)
         aaaCACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATc                                                                                    >  1:3732217/2‑100 (MQ=255)
                    atatTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCaa                                                                         <  1:4442668/100‑1 (MQ=255)
                    atatTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCaa                                                                         >  1:5270428/1‑100 (MQ=255)
                    atatTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTACaa                                                                         >  1:2063660/1‑100 (MQ=255)
                      atTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCacca                                                                       >  1:1681703/1‑97 (MQ=255)
                      atTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCaaaa                                                                       >  1:1474721/1‑100 (MQ=255)
                      atTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTCCtaaa                                                                       >  1:894861/1‑100 (MQ=255)
                        tGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAACGCATATATCTTTAACTTCAAAata                                                                     >  1:5422960/1‑99 (MQ=255)
                                         ttCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGTAtt                                                    >  1:5147065/1‑100 (MQ=255)
                                                                   gTTCCTAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGTATTCTTTGATTTGATTGATTTGTCCGTCg                          >  1:3151257/1‑100 (MQ=255)
                                                                     tCATAGATTTACCTGAAATACCACGATCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGAATTCTTTGATTTGATTGATTTGTCCGTCGAt                        >  1:1637633/1‑100 (MQ=255)
                                                                                         ccACGATCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGTATTCTTTGATTTGATTGATTTGTCCGTCGATACTATAACCTTCTGTGCttt    <  1:4250771/100‑1 (MQ=255)
                                                                                           aCGATCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGTATTCTTTGATTTGATTGATTTGTCCGTCGATACTATAACCTTCTGTGCTTTGc  <  1:1492124/100‑1 (MQ=255)
                                                                                              |                                                                                                
TGTATACCATAACACAGTCGATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGGTCAGCGTATATATCTTTAACTTCAAAATGATGGAAGTCACAGTATTCTTTGATTTGATTGATTTGTCCGTCGATACTATAACCTTCTGTGCTTTGC  >  NC_007793/47274‑47464

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