Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,090,282 C→T intergenic (+30/‑65) yggI → / → endA conserved protein YggI/DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,2820CT100.0% 135.3 / NA 34intergenic (+30/‑65)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/18);  total (16/18)

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCT  >  NC_000913/3090236‑3090325
                                              |                                           
aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                           <  1:1982873/49‑1 (MQ=255)
     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                       <  1:469162/48‑1 (MQ=255)
      ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                     >  1:610475/1‑49 (MQ=255)
         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                   >  1:1903197/1‑48 (MQ=255)
          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                 >  1:1067026/1‑49 (MQ=255)
          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                 <  1:28725/49‑1 (MQ=255)
           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                 >  1:1945968/1‑48 (MQ=255)
              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                              >  1:216781/1‑48 (MQ=255)
              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                              >  1:2044524/1‑48 (MQ=255)
              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                              >  1:1792461/1‑48 (MQ=255)
                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                            <  1:1150431/48‑1 (MQ=255)
                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                            >  1:784802/1‑48 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                          <  1:1110365/49‑1 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                           >  1:581147/1‑48 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         <  1:983304/49‑1 (MQ=255)
                      tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                     >  1:1661181/1‑49 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                      <  1:1299117/46‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                      <  1:1586469/46‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                      <  1:211181/46‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                   >  1:659158/1‑49 (MQ=255)
                         aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  <  1:183332/49‑1 (MQ=255)
                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                 >  1:1635370/1‑49 (MQ=255)
                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                 <  1:1641075/49‑1 (MQ=255)
                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                 <  1:709536/49‑1 (MQ=255)
                           tATCAGGAACTTTCCTGATTTGGCTGTTTGCATACCAAAACAGCTTTCg                >  1:1225175/1‑49 (MQ=255)
                           tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                >  1:1890686/1‑49 (MQ=255)
                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt              <  1:898042/49‑1 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc            <  1:524757/49‑1 (MQ=255)
                                 gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          >  1:1377405/1‑49 (MQ=255)
                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt         <  1:146206/49‑1 (MQ=255)
                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt         <  1:1589853/49‑1 (MQ=255)
                                    cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc       <  1:1270104/49‑1 (MQ=255)
                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg    >  1:272707/1‑49 (MQ=255)
                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt  <  1:624449/49‑1 (MQ=255)
                                              |                                           
AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCT  >  NC_000913/3090236‑3090325

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