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% 64.5 / 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 . (8/9);  total (8/9)

CTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTG  >  NC_000913/3742465‑3742644
                                                                                                                                     |                                              
cTTCCGGCTGGACGCTCTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGa                                         >  2:1624247/1‑139 (MQ=255)
              ctctTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCttt                            <  2:1045950/139‑1 (MQ=255)
                          cgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcatca                <  2:1253721/139‑1 (MQ=255)
                          cgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcatca                <  1:418167/139‑1 (MQ=255)
                           gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaa               >  1:1722535/1‑139 (MQ=255)
                               tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcct           >  2:941714/1‑137 (MQ=255)
                               tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcct           >  2:1648259/1‑137 (MQ=255)
                                   gTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttg       <  1:1343172/139‑7 (MQ=255)
                                        ctTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgacaat  <  1:941714/139‑12 (MQ=255)
                                                                 cATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                             >  1:798324/1‑69 (MQ=255)
                                                                 cATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                             <  2:798324/72‑4 (MQ=255)
                                                                             tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaa                        <  1:2375/80‑1 (MQ=255)
                                                                             tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaa                        >  2:2375/1‑80 (MQ=255)
                                                                                        gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt        <  2:1464761/85‑6 (MQ=255)
                                                                                        gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt        >  1:1464761/1‑80 (MQ=255)
                                                                                         cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt        <  1:2851/84‑6 (MQ=255)
                                                                                         cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt        >  2:2851/1‑79 (MQ=255)
                                                                                          gcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                       >  2:1225117/1‑68 (MQ=255)
                                                                                          gcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                       <  1:1225117/68‑1 (MQ=255)
                                                                                                                                     |                                              
CTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTG  >  NC_000913/3742465‑3742644

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