Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 46,305 T→G I406L (ATA→CTA)  SAUSA300_RS00200 ← recombinase RecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00779346,3050TG100.0% 29.6 / NA 13I406L (ATA→CTA) SAUSA300_RS00200recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (5/8);  total (5/8)

TTGAATACTTGGTCTGATTACTGATGTTAAGTCCGGATTATCCTCGATGGTTTTAATCAAGTTATTTATTTTGATTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTAC  >  NC_007793/46237‑46401
                                                                    |                                                                                                
ttGAATACTTGGTCTGATTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGc                                                                   <  1:1157303/100‑1 (MQ=255)
  gctcttccgatCTGATTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCtt                                                                 <  1:151088/90‑1 (MQ=255)
  gAATACTTGGTCTGATTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCtt                                                                 <  1:2728753/100‑1 (MQ=255)
     cttccgatCTGATTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTtata                                                              <  1:4538772/93‑1 (MQ=255)
         tGGTCTGATTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCa                                                          <  1:4569899/100‑1 (MQ=255)
               gATTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCg                                                    >  1:2346698/1‑100 (MQ=255)
                aTTACTGATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTAGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGt                                                   <  1:3334506/100‑1 (MQ=255)
                      gATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGaa                                             >  1:2996864/1‑100 (MQ=255)
                      gATGTTAAGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGaa                                             >  1:3552329/1‑100 (MQ=255)
                            aaGTCCGGATTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTATACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAg                                       >  1:5379926/1‑100 (MQ=255)
                                    aTTATCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCa                               <  1:5024940/100‑1 (MQ=255)
                                        tCCTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACta                           <  1:4140197/100‑2 (MQ=255)
                                         ccTCTATGGTTTTAATCAAGTTATTTAGTTTGGTTTGTACTTCATCATATTGTTGTTGCTTGTAAGCAATATCGTGATGAAGTGCAGCGCCATCAACttt                          >  1:4801495/1‑99 (MQ=255)
                                                               atttagttTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAAcacgtgct    <  1:4463483/100‑8 (MQ=255)
                                                                 ttagttTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCtac  <  1:215338/100‑1 (MQ=255)
                                                                 ttagttTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCtac  <  1:963282/100‑1 (MQ=255)
                                                                    |                                                                                                
TTGAATACTTGGTCTGATTACTGATGTTAAGTCCGGATTATCCTCGATGGTTTTAATCAAGTTATTTATTTTGATTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTAC  >  NC_007793/46237‑46401

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