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% 168.5 / NA 43intergenic (+30/‑65)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/20);  total (23/20)

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090236‑3090327
                                              |                                             
aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                             >  1:2266552/1‑49 (MQ=255)
    ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                          >  1:2146994/1‑48 (MQ=255)
     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                        >  1:414974/1‑49 (MQ=255)
       tGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      >  1:694604/1‑49 (MQ=255)
        gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      <  1:1350508/48‑1 (MQ=255)
        gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      <  1:1838885/48‑1 (MQ=255)
         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                     <  1:262022/48‑1 (MQ=255)
         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                    <  1:1727480/49‑1 (MQ=255)
           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:369724/1‑48 (MQ=255)
             aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                >  1:1477772/1‑49 (MQ=255)
             aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                >  1:2040731/1‑49 (MQ=255)
              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                               <  1:1253613/49‑1 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:1165565/49‑1 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:806375/49‑1 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:1515225/49‑1 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                             >  1:312356/1‑48 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:1711048/1‑49 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:879752/1‑49 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:2272841/1‑49 (MQ=255)
                    tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:1841677/1‑47 (MQ=255)
                    tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:1968067/1‑47 (MQ=255)
                      tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                       >  1:1693407/1‑49 (MQ=255)
                      tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                       >  1:132984/1‑49 (MQ=255)
                       gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      >  1:1487530/1‑49 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                     >  1:1742500/1‑49 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      <  1:234693/48‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      <  1:2064429/48‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      <  1:177687/48‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      <  1:1473527/48‑1 (MQ=255)
                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   >  1:2247333/1‑49 (MQ=255)
                          cTATCAGGAACTTTCCTGATTGGGCTGATTGCATACCAAAACAGCTTTc                   >  1:787045/1‑49 (MQ=39)
                           tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                  >  1:553434/1‑49 (MQ=255)
                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                <  1:652406/49‑1 (MQ=255)
                             tCAGGAACTTTCCTCATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 <  1:1126133/48‑1 (MQ=39)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:1630005/49‑1 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:919834/49‑1 (MQ=255)
                                ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg             >  1:2258537/1‑49 (MQ=255)
                                ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg             <  1:1605510/49‑1 (MQ=255)
                                    cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc         >  1:281162/1‑49 (MQ=255)
                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt        <  1:331166/47‑1 (MQ=255)
                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       <  1:540581/48‑1 (MQ=255)
                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc     <  1:244384/48‑1 (MQ=255)
                                           gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg  >  1:1951643/1‑49 (MQ=255)
                                              |                                             
AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090236‑3090327

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