Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,716 C→A A186E (GCA→GAA)  bglJ → DNA‑binding transcriptional regulator BglJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,7160CA100.0% 28.8 / NA 14A186E (GCA→GAA) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (6/8);  total (6/8)

TGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTA  >  NC_000913/4604592‑4604833
                                                                                                                             |                                                                                                                      
ggggaaaTTACCCAAAGCGGAATTTTAGCCCGGCGGAGAGGGAAATTTTTCGCTTTTTGGCCCTTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATa                                                                                                                   <  1:449377/123‑1 (MQ=255)
        tAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTaa                                                                                                            >  1:623308/1‑129 (MQ=255)
         aaCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCCAAACGATCCGTGAACATAAATTTAAt                                                                                                           >  1:312480/1‑129 (MQ=255)
                                 aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGAGGTCGAAACTGGGCgcaggt                                                                                   >  1:203686/1‑124 (MQ=255)
                                          aaaaTAATGCGGCGTTAGGTCGTGGGGCAACTGAAGACACAAATTGCCTGTCAGCTTAAACGCAATATCAAAACGTTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa                                                                          <  1:535825/128‑1 (MQ=255)
                                                                          aaaTGACACAAATTGCCGGGCAGCTTAAACGCAATGTCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                           <  1:287279/128‑1 (MQ=255)
                                                                            aTGACACAAATTGCCCGGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgt                                         <  1:601949/129‑1 (MQ=255)
                                                                                   aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTTAACATAAATTTAATGTGATGTCGCAACTGGGCGTCAGTTCTGACGCAGGGTTGTTCGAGGCCGCAGCAATTCTGTTATTTATGccc                                   >  1:459428/1‑126 (MQ=255)
                                                                                    aatttgcGAGCAACTTAAACGCAAAATCAAAACGGTCCCTGAAAATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCa                                 <  1:267156/123‑1 (MQ=255)
                                                                                                       cgaaaTACCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg             <  1:629961/126‑1 (MQ=255)
                                                                                                          aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCGGTTATGACCCGGGTTTGTGGGGGCCGCAGACATTCCTTTCAGATTGGGGCATTGGGCAACAAGATAAGGTGttc          >  1:152146/1‑128 (MQ=255)
                                                                                                          aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt            <  1:322765/127‑1 (MQ=255)
                                                                                                            atacAAAACGGTCCGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGc         <  1:240807/125‑1 (MQ=255)
                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAACCAATTAAGTGGTTTCATCCCTa  >  1:600356/1‑129 (MQ=255)
                                                                                                                             |                                                                                                                      
TGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTA  >  NC_000913/4604592‑4604833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: