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% 114.9 / NA 29intergenic (‑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/17);  total (12/17)

GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTT  >  NC_000913/3742463‑3742646
                                                                                                                                       |                                                
ggCTTCCGGCTGGACGCTCTTTTTCCTGCGCGATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               >  2:1201747/1‑136 (MQ=255)
ggCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                               >  2:1069341/1‑136 (MQ=255)
       ggCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGAtttt                                       >  2:481195/1‑139 (MQ=255)
         cTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTa                                     >  1:692318/1‑139 (MQ=255)
                ctctTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCttt                              <  1:1265644/139‑1 (MQ=255)
                     tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGaaa                         <  2:1634118/139‑1 (MQ=255)
                       tCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAAtt                       <  1:810065/139‑1 (MQ=255)
                          tGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    <  2:1346253/139‑1 (MQ=255)
                           gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcatc                   >  1:910803/1‑139 (MQ=255)
                            cgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcatca                  <  2:910803/139‑1 (MQ=255)
                             gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    <  1:390232/136‑1 (MQ=255)
                             gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    >  2:390232/1‑136 (MQ=255)
                             gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaa                 >  2:1429806/1‑139 (MQ=255)
                                 tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcct             >  1:1336923/1‑137 (MQ=255)
                                 tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               >  1:1334100/1‑137 (MQ=255)
                                 tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatc               <  2:1334100/137‑1 (MQ=255)
                                    cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt          <  1:1201747/139‑6 (MQ=255)
                                    cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt          <  2:231443/139‑6 (MQ=255)
                                    cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt          <  2:692318/139‑6 (MQ=255)
                                       tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgac       <  2:191026/139‑9 (MQ=255)
                                            ttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgacaatta  <  1:295347/139‑14 (MQ=255)
                                                               tCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    >  1:564911/1‑102 (MQ=255)
                                                               tCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    <  2:564911/102‑1 (MQ=255)
                                                                   cATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt                                              >  2:1549344/1‑69 (MQ=255)
                                                                   cATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt                                              <  1:1549344/73‑5 (MQ=255)
                                                                               tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    <  2:1302067/86‑1 (MQ=255)
                                                                               tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATtcat                    >  1:1302067/1‑86 (MQ=255)
                                                                                        tAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctg            >  2:11714/1‑82 (MQ=255)
                                                                                        tAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctg            <  1:11714/85‑4 (MQ=255)
                                                                                           cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt          >  1:528947/1‑79 (MQ=255)
                                                                                           cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt          <  2:528947/84‑6 (MQ=255)
                                                                                                aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTAt                                    <  2:1541241/53‑1 (MQ=255)
                                                                                                aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTAt                                    >  1:1541241/1‑53 (MQ=255)
                                                                                                                                       |                                                
GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTT  >  NC_000913/3742463‑3742646

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