Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,941,650 A→G intergenic (‑323/‑158) yieP ← / → rrsC putative transcriptional regulator YieP/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,941,6500AG100.0% 51.0 / NA 17intergenic (‑323/‑158)yieP/rrsCputative transcriptional regulator YieP/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/7);  total (10/7)

GGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTC  >  NC_000913/3941575‑3941725
                                                                           |                                                                           
ggCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaa                                                                        <  2:730892/81‑1 (MQ=37)
  cAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaaag                                                                      <  2:730893/81‑1 (MQ=255)
   agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaaagc                                                                     <  2:730894/81‑1 (MQ=255)
     agaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCgg                                                                   >  1:730917/1‑81 (MQ=255)
     agaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCgg                                                                   <  1:730895/81‑1 (MQ=255)
         aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACt                                                               <  1:730896/81‑1 (MQ=255)
          gCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTg                                                              >  2:730914/1‑81 (MQ=255)
             aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATTCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGctc                                                           <  1:730897/81‑1 (MQ=255)
                  aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTaa                                                      >  1:730921/1‑81 (MQ=255)
                         tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTa                                               >  2:730922/1‑81 (MQ=255)
                         tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTa                                               >  1:730919/1‑81 (MQ=255)
                          gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTa                                               <  1:730898/80‑1 (MQ=255)
                                 tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                                       >  2:730920/1‑81 (MQ=255)
                                                ttatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                        <  1:730899/81‑1 (MQ=18)
                                                   tGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGa                      <  1:730900/80‑1 (MQ=21)
                                                           ccGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa             <  2:730901/81‑1 (MQ=18)
                                                             gcgcCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAtt           >  1:730923/1‑81 (MQ=18)
                                                                ccgcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCt         >  1:730928/1‑80 (MQ=21)
                                                                  gcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTaa      >  1:730924/1‑81 (MQ=18)
                                                                      agaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTc  >  2:730957/1‑81 (MQ=18)
                                                                           |                                                                           
GGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTC  >  NC_000913/3941575‑3941725

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