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

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGC  >  NC_000913/3090236‑3090324
                                              |                                          
aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                          >  1:815413/1‑49 (MQ=255)
 gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                         >  1:356159/1‑49 (MQ=255)
 gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                         <  1:2557650/49‑1 (MQ=255)
   tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                       >  1:195632/1‑49 (MQ=255)
   tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                       >  1:373387/1‑49 (MQ=255)
    ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                       >  1:2391384/1‑48 (MQ=255)
     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                       <  1:2361157/47‑1 (MQ=255)
      ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                    <  1:1751986/49‑1 (MQ=255)
        gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                   <  1:2674115/48‑1 (MQ=255)
        gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                  >  1:731064/1‑49 (MQ=255)
         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                  <  1:2127179/48‑1 (MQ=255)
         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                 <  1:1694135/49‑1 (MQ=255)
         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                 <  1:1701148/49‑1 (MQ=255)
          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                <  1:345894/49‑1 (MQ=255)
          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                <  1:582913/49‑1 (MQ=255)
          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                <  1:1672620/49‑1 (MQ=255)
           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                >  1:976904/1‑48 (MQ=255)
           aaaTAACCATCAGAACTATCAGGAACTTTCCTGATTTGGCTGAATGCa                                >  1:1424370/1‑48 (MQ=38)
            aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                              <  1:2003323/49‑1 (MQ=255)
            aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                              <  1:1019735/49‑1 (MQ=255)
             aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                             >  1:1769887/1‑49 (MQ=255)
             aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                             <  1:134635/49‑1 (MQ=255)
              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                             >  1:2038182/1‑48 (MQ=255)
                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                          <  1:1151803/49‑1 (MQ=255)
                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                          <  1:2331554/49‑1 (MQ=255)
                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                           <  1:2209147/48‑1 (MQ=255)
                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                           >  1:654127/1‑48 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                         <  1:2197496/49‑1 (MQ=255)
                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                          >  1:782486/1‑48 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        >  1:1842393/1‑49 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        >  1:2732364/1‑49 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        >  1:1935256/1‑49 (MQ=255)
                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        <  1:1088850/49‑1 (MQ=255)
                   aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        >  1:2426969/1‑48 (MQ=255)
                   aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        >  1:2042481/1‑48 (MQ=255)
                    tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                        >  1:1406514/1‑47 (MQ=255)
                     cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                     >  1:87028/1‑49 (MQ=255)
                     cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                     >  1:1820656/1‑49 (MQ=255)
                      tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                    >  1:821307/1‑49 (MQ=255)
                       gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                   >  1:1545789/1‑49 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                      <  1:1958971/45‑1 (MQ=255)
                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                      <  1:568351/45‑1 (MQ=255)
                         aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                 <  1:1528676/49‑1 (MQ=255)
                         aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                 <  1:196714/49‑1 (MQ=255)
                         aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                 <  1:1855359/49‑1 (MQ=255)
                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                >  1:1695729/1‑49 (MQ=255)
                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                <  1:1264977/49‑1 (MQ=255)
                           tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg               >  1:1121769/1‑49 (MQ=255)
                           tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg               >  1:665480/1‑49 (MQ=255)
                           tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg               >  1:1490381/1‑49 (MQ=255)
                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg               <  1:923371/47‑1 (MQ=255)
                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc              <  1:2375994/48‑1 (MQ=255)
                              cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa            >  1:2813538/1‑49 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa            <  1:955523/48‑1 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa            <  1:543559/48‑1 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc           <  1:1760449/49‑1 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc           >  1:1340369/1‑49 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc           <  1:1525252/49‑1 (MQ=255)
                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc           <  1:1789916/49‑1 (MQ=255)
                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt        >  1:1155801/1‑49 (MQ=255)
                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt        <  1:1462159/49‑1 (MQ=255)
                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt        >  1:1024276/1‑49 (MQ=255)
                                   aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg       >  1:719131/1‑49 (MQ=255)
                                     tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt     >  1:1332117/1‑49 (MQ=255)
                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt     <  1:2365562/47‑1 (MQ=255)
                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg    <  1:1432741/48‑1 (MQ=255)
                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg    <  1:967587/48‑1 (MQ=255)
                                        ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg   <  1:2791995/48‑1 (MQ=255)
                                        ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc  >  1:2830660/1‑49 (MQ=255)
                                        ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc  >  1:1601649/1‑49 (MQ=255)
                                              |                                          
AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGC  >  NC_000913/3090236‑3090324

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: