BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_012967 3,339,158 A→C 100% intergenic (+22/‑4) dusB → / → fis tRNA dihydrouridine synthase DusB/DNA‑binding transcriptional regulator Fis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129673,339,1580AC100.0% 44.3 / NA 16intergenic (+22/‑4)dusB/fistRNA dihydrouridine synthase DusB/DNA‑binding transcriptional regulator Fis
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (5/11);  total (5/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTC  >  NC_012967/3339128‑3339184
                              |                          
tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                       <  1:768788/36‑1 (MQ=38)
tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                       >  1:89192/1‑36 (MQ=38)
  tGCGTAAACAGAAATAAAGAGCTGACAGCACTATGt                     <  1:339432/36‑1 (MQ=39)
   acGTAAACAGAAATAAAGAGCTGACAGCACTATGtt                    <  1:919511/35‑1 (MQ=38)
    cGTAAACAGAAATAAAGAGCTGACAGCACCATGtcc                   >  1:366208/1‑34 (MQ=25)
     gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                  <  2:133926/36‑1 (MQ=39)
      tAAACAGAAATAAAGAGCTGACAGCACTATGTTCGa                 <  2:413743/36‑1 (MQ=39)
       aaaCAGAAATAAAGAGCTGACAGCACTATGTTCGaa                <  2:152870/36‑1 (MQ=39)
       aaaCAGAAATAAAGAACTCACAGCACTATGTTCGaa                <  1:466994/36‑1 (MQ=16)
          cAGAAATAAAAAACTGACAGCACTATGTTCGaacaa             <  1:551724/36‑1 (MQ=19)
           aGAAATAAAGAGCTGACAGCACTATGTTCGaacaac            >  1:194947/1‑36 (MQ=39)
               aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt        <  1:256298/36‑1 (MQ=39)
                 aaaGTGCTGACAGCACTATGTTCGAACAACGCGTaa      <  1:713281/36‑1 (MQ=37)
                 aaaGGGCTGACAGCACTATGTTCGAACAACGCGTaa      <  1:517003/36‑1 (MQ=37)
                    gagCTGACAGCACTATGTTCGAACAACGCGTAAAtt   >  2:166040/1‑36 (MQ=39)
                     agCTGACAGCACTATGTTCGAACAACGCGTAAATTc  >  1:800911/1‑36 (MQ=39)
                              |                          
TTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTC  >  NC_012967/3339128‑3339184

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: