Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,090,341 Δ1 bp intergenic (+89/‑6) 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,3410T.100.0% 276.3 / NA 56intergenic (+89/‑6)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (23/33);  total (23/33)

CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090297‑3090384
                                            |                                           
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                        <  1:91402/49‑1 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                        <  1:1870587/49‑1 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                        <  1:869015/49‑1 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                        <  1:694646/49‑1 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      <  1:1112808‑M1/49‑3 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      >  1:2168677‑M1/1‑47 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      >  1:198254‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      >  1:774566‑M1/1‑46 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                    <  1:991103‑M1/49‑5 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                    <  1:2122681‑M1/49‑5 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                    <  1:62637‑M1/49‑5 (MQ=255)
     aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                   <  1:2121627‑M1/49‑6 (MQ=255)
     aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                   <  1:1223264‑M1/49‑6 (MQ=255)
      cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  >  1:1261387‑M1/1‑43 (MQ=255)
        gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                <  1:181751‑M1/49‑9 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                             <  1:2251742‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                             <  1:879167‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                             <  1:1960360‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                              <  1:170101‑M1/47‑11 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                            <  1:2316303‑M1/48‑13 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                          <  1:1541828‑M1/49‑15 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                           >  1:1253617‑M1/1‑35 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                        >  1:1429019‑M1/1‑33 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                        >  1:2173992‑M1/1‑33 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                        >  1:2155308‑M1/1‑33 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                      <  1:1091982‑M1/49‑19 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                       <  1:2099764‑M1/48‑18 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                        <  1:516817‑M1/47‑17 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                     >  1:383009‑M1/1‑30 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                    <  1:1035962‑M1/49‑21 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                   >  1:1282529‑M1/1‑28 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                   >  1:849801‑M1/1‑28 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                   >  1:1444810‑M1/1‑28 (MQ=255)
                      gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  >  1:458226‑M1/1‑27 (MQ=255)
                      gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                   <  1:2261357‑M1/48‑22 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                 <  1:1583815‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                 <  1:1406068‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  <  1:1304063‑M1/48‑23 (MQ=255)
                        tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                >  1:1099943‑M1/1‑25 (MQ=255)
                        tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                >  1:2256332‑M1/1‑25 (MQ=255)
                        tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                >  1:2313919‑M1/1‑25 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt              <  1:612036‑M1/49‑27 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                <  1:1354527‑M1/47‑25 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                <  1:2073168‑M1/47‑25 (MQ=255)
                           cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt              >  1:730186‑M1/1‑22 (MQ=255)
                           cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt              >  1:2150984‑M1/1‑22 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc          <  1:1764350‑M1/49‑31 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc          <  1:903927‑M1/49‑31 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac           <  1:447368‑M1/48‑30 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattcttcttgttaac           <  1:823479‑M1/48‑30 (MQ=255)
                                tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc        >  1:1559602‑M1/1‑17 (MQ=255)
                                   aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctac     >  1:2000255‑M1/1‑14 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc   >  1:2288225‑M1/1‑12 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc   >  1:2239688‑M1/1‑12 (MQ=255)
                                      acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc  <  1:2295882‑M1/49‑39 (MQ=255)
                                      acGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc   >  1:1533532‑M1/1‑11 (MQ=255)
                                            |                                           
CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090297‑3090384

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