Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,742,599 Δ1 bp intergenic (‑67/‑134) yiaJ ← / → yiaK transcriptional repressor for the yiaKLMNO‑lyxK‑sgbHUE operon/2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,742,5980C.100.0% 66.2 / NA 17intergenic (‑66/‑135)yiaJ/yiaKtranscriptional repressor for the yiaKLMNO‑lyxK‑sgbHUE operon/2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/7);  total (10/7)

CCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTT  >  NC_000913/3742468‑3742646
                                                                                                                                  |                                                
ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATACTTTATAAcaatt                                              >  2:842225/1‑131 (MQ=255)
ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt                                              <  1:842225/135‑5 (MQ=255)
ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAtt                                         >  2:85451/1‑139 (MQ=255)
ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAtt                                         >  1:311402/1‑139 (MQ=255)
  ggCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAtttt                                       >  1:16921/1‑139 (MQ=255)
                tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                         <  2:311402/139‑1 (MQ=255)
                          cATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTAt                                    >  1:1321776/1‑118 (MQ=255)
                          cATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTAt                                    <  2:1321776/118‑1 (MQ=255)
                             ctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCagcaatcctg            >  1:830962/1‑136 (MQ=255)
                                 ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttga        >  1:628418/1‑132 (MQ=255)
                                       ttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAACGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgacaatta  <  2:906428/139‑14 (MQ=255)
                                                aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                         >  1:820560/1‑107 (MQ=255)
                                                aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                         <  2:820560/107‑1 (MQ=255)
                                                                     tCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAAtt                       <  1:71768/88‑1 (MQ=255)
                                                                     tCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAAtt                       >  2:71768/1‑88 (MQ=255)
                                                                                      cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               >  2:122758/1‑45 (MQ=255)
                                                                                      cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               <  1:122758/48‑4 (MQ=255)
                                                                                      cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               <  2:1108017/79‑1 (MQ=255)
                                                                                      cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               >  1:1108017/1‑79 (MQ=255)
                                                                                               tgCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt                                              <  1:661438/40‑5 (MQ=255)
                                                                                               tgCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt                                              >  2:661438/1‑36 (MQ=255)
                                                                                                         cACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               <  2:599215/60‑1 (MQ=255)
                                                                                                         cACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               >  1:599215/1‑60 (MQ=255)
                                                                                                                                  |                                                
CCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTT  >  NC_000913/3742468‑3742646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: