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.93.3% 50.7 / ‑2.3 15intergenic (‑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/1);  new base . (6/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATT  >  NC_000913/3742463‑3742638
                                                                                                                                       |                                        
ggCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaat                                       <  2:264798/139‑4 (MQ=255)
          tGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt                            <  2:1075093/139‑1 (MQ=255)
           ggACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATg                           <  2:320531/139‑1 (MQ=255)
                       tCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAAtt               <  2:74751/139‑1 (MQ=255)
                          tGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcat            <  1:1355658/139‑1 (MQ=255)
                           gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatc           >  2:1224814/1‑139 (MQ=255)
                             gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaa         >  2:799237/1‑139 (MQ=255)
                                    cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt  <  2:119438/139‑6 (MQ=255)
                                                    cATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa                 >  1:471824/1‑108 (MQ=255)
                                                    cATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa                 <  2:471824/108‑1 (MQ=255)
                                                            cttTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca          <  1:993168/107‑1 (MQ=255)
                                                            cttTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca          >  2:993168/1‑107 (MQ=255)
                                                                                                aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACC‑ATTGATTTTTATGTCTTTTGaa                  <  1:1094505/63‑1 (MQ=255)
                                                                                                aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaa                  >  2:1094505/1‑63 (MQ=255)
                                                                                                                     tcGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatcctg    <  1:508367/56‑4 (MQ=255)
                                                                                                                     tcGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatcctg    >  2:508367/1‑53 (MQ=255)
                                                                                                                                       |                                        
GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATT  >  NC_000913/3742463‑3742638

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