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% 95.5 / NA 24intergenic (‑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 . (12/12);  total (12/12)

GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTT  >  NC_000913/3742463‑3742646
                                                                                                                                       |                                                
ggCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               <  2:747465/139‑4 (MQ=255)
  cTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTg                                            <  2:951384/139‑1 (MQ=255)
   ttCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGa                                           >  2:254213/1‑139 (MQ=255)
    tCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAt                                          >  1:744260/1‑139 (MQ=255)
      cGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAttt                                        >  2:886281/1‑139 (MQ=255)
        gCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAttttt                                      <  1:254213/139‑1 (MQ=255)
               gCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCtt                               <  1:886281/139‑1 (MQ=255)
                     tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                         <  1:1000199/139‑1 (MQ=255)
                           gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGCAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcatc                   >  2:600569/1‑139 (MQ=255)
                             gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaa                 >  1:945531/1‑139 (MQ=255)
                                 tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGTTGGTAGCGCAAAGTGTTCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcct             >  2:1050813/1‑137 (MQ=255)
                                 tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATTGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtc                      <  2:1018852/130‑1 (MQ=255)
                                 tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtc                      >  1:1018852/1‑130 (MQ=255)
                                    cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt          <  2:945531/139‑6 (MQ=255)
                                      ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGt                                  >  2:39219/1‑113 (MQ=255)
                                      ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGt                                  <  1:39219/113‑1 (MQ=255)
                                      ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttga        <  1:600569/139‑8 (MQ=255)
                                           tttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgacaatt   >  1:65953/1‑127 (MQ=255)
                                            ttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgacaatta  <  2:561679/139‑14 (MQ=255)
                                                   ccATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgt           <  2:106507/123‑5 (MQ=255)
                                                   ccATCACTTCTTTCCCCATTTTGTCGCGTCCTGAGGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgt           >  1:106507/1‑119 (MQ=255)
                                                                                         aGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               >  2:166640/1‑47 (MQ=255)
                                                                                         aGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               <  1:166640/50‑4 (MQ=255)
                                                                                                    tgCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               <  2:484376/70‑1 (MQ=255)
                                                                                                    tgCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               >  1:484376/1‑70 (MQ=255)
                                                                                                            ttCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    <  1:984574/57‑1 (MQ=255)
                                                                                                            ttCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    >  2:984574/1‑57 (MQ=255)
                                                                                                                                       |                                                
GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTT  >  NC_000913/3742463‑3742646

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