Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 251,954 T→C intergenic (+1/‑51) dinB → / → yafN DNA polymerase IV/antitoxin YafN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913251,9540TC100.0% 75.8 / NA 22intergenic (+1/‑51)dinB/yafNDNA polymerase IV/antitoxin YafN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/14);  total (8/14)

CGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAATTCTCGCTGAAAAATC  >  NC_000913/251878‑252030
                                                                            |                                                                            
cgGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGta                                                                           >  1:55119/1‑80 (MQ=255)
   tgtgCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACt                                                                       <  1:55089/81‑1 (MQ=255)
            ggtggGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGta                                                               >  2:55120/1‑80 (MQ=255)
             gtggGGCTGAATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGt                                                                <  1:55090/78‑1 (MQ=255)
              tggGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTatat                                                            <  2:55091/81‑1 (MQ=255)
                      cATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGgtgt                                                    <  2:55092/81‑1 (MQ=255)
                      cATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGgt                                                      <  1:55093/79‑1 (MQ=255)
                        tgtgACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGtgtg                                                   >  2:55138/1‑80 (MQ=255)
                                 gCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATtattat                                          <  1:55094/80‑1 (MQ=255)
                                      aCCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGt                                    <  1:55095/81‑1 (MQ=255)
                                       cccGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTa                                   <  2:55096/81‑1 (MQ=255)
                                       cccGCAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTa                                   >  2:55128/1‑81 (MQ=255)
                                           cAAATGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCAc                               <  2:55097/81‑1 (MQ=255)
                                              aTGGAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGCGCATTATTATGAGGGTATCACTg                             >  1:55124/1‑80 (MQ=255)
                                                ggAAAGACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTAt                          >  1:55137/1‑81 (MQ=255)
                                                     gACAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATc                     >  1:55131/1‑81 (MQ=255)
                                                      aCAACTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCg                    <  2:55098/81‑1 (MQ=255)
                                                          cTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAAtt                <  2:55099/81‑1 (MQ=255)
                                                          cTGGTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAAtt                <  1:245689/81‑1 (MQ=255)
                                                             gTGCTGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAATtct              <  1:55100/80‑1 (MQ=255)
                                                                 tGGGATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAATTCTCGCt          >  1:55125/1‑80 (MQ=255)
                                                                   ggATTATGACGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAATTCTCGCTGa        <  2:55101/80‑1 (MQ=255)
                                                                        atgacgTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAATTCTCGCTGAAAAATc  >  2:55129/1‑81 (MQ=255)
                                                                            |                                                                            
CGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATACTATTATGTATATTCTGGTGTGCATTATTATGAGGGTATCACTGTATGCATCGAATTCTCGCTGAAAAATC  >  NC_000913/251878‑252030

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