Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,579 G→A intergenic (‑30/‑97) rph ← / → yicC truncated RNase PH/UPF0701 family protein YicC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,5790GA100.0% 97.8 / NA 25intergenic (‑30/‑97)rph/yicCtruncated RNase PH/UPF0701 family protein YicC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/13);  total (12/13)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGA  >  NC_000913/3816533‑3816620
                                              |                                         
cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                         <  1:720382/49‑1 (MQ=255)
 ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                        <  1:1281121/49‑1 (MQ=255)
  gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                       >  1:1425165/1‑49 (MQ=255)
   ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                      <  1:1981485/49‑1 (MQ=255)
       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                    <  1:144268/47‑1 (MQ=255)
       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                    <  1:1947953/47‑1 (MQ=255)
       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                   >  1:1765343/1‑48 (MQ=255)
        tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                 >  1:1608402/1‑49 (MQ=255)
        tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                 >  1:1926321/1‑49 (MQ=255)
         ggACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                 >  1:151542/1‑48 (MQ=255)
          gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                <  1:1092559/48‑1 (MQ=255)
           aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                              <  1:575478/49‑1 (MQ=255)
           aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                              <  1:31956/49‑1 (MQ=255)
                tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                         >  1:1443015/1‑49 (MQ=255)
                 atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                        >  1:1942372/1‑49 (MQ=255)
                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                        <  1:77299/48‑1 (MQ=255)
                      aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                   >  1:1764810/1‑49 (MQ=255)
                          tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                >  1:1586297/1‑48 (MQ=255)
                            cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                 <  1:1709019/45‑1 (MQ=255)
                            cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt             >  1:120185/1‑49 (MQ=255)
                             ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt            <  1:450662/49‑1 (MQ=255)
                                      cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt    >  1:707690/1‑48 (MQ=255)
                                       aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa  <  1:443344/49‑1 (MQ=255)
                                       aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa  <  1:1599863/49‑1 (MQ=255)
                                        aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa  >  1:2024322/1‑48 (MQ=255)
                                              |                                         
CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGA  >  NC_000913/3816533‑3816620

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