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% 292.3 / NA 59intergenic (+89/‑6)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (31/28);  total (31/28)

CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGA  >  NC_000913/3090297‑3090383
                                            |                                          
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                       <  1:2120686/49‑1 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                       >  1:1539309/1‑49 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                     >  1:191226‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                     >  1:1557163‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                     >  1:684656‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                     >  1:2277728‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                     >  1:2440720‑M1/1‑46 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                   <  1:2555376‑M1/49‑5 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                   >  1:2301866‑M1/1‑45 (MQ=255)
     aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                  <  1:1896751‑M1/49‑6 (MQ=255)
      cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                 >  1:8594‑M1/1‑43 (MQ=255)
       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                >  1:2535283‑M1/1‑42 (MQ=255)
       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                >  1:984302‑M1/1‑42 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                              <  1:2679502‑M1/49‑10 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                              <  1:2019513‑M1/49‑10 (MQ=255)
          tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                             >  1:1989169‑M1/1‑39 (MQ=255)
           ttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                            >  1:1500304‑M1/1‑38 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                            <  1:2187156‑M1/48‑12 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                           <  1:261117‑M1/48‑13 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                           <  1:225230‑M1/48‑13 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                           <  1:1593103‑M1/48‑13 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                         <  1:2127402‑M1/49‑15 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                          <  1:1323765‑M1/48‑14 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                          >  1:120225‑M1/1‑35 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                          >  1:2424554‑M1/1‑35 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                          >  1:839449‑M1/1‑35 (MQ=255)
               cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                        >  1:124232‑M1/1‑34 (MQ=255)
               cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                        >  1:1516843‑M1/1‑34 (MQ=255)
               cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                        >  1:2396651‑M1/1‑34 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                       >  1:1316922‑M1/1‑33 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                     <  1:603890‑M1/49‑19 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                      <  1:770427‑M1/48‑18 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                    <  1:438713‑M1/49‑20 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                   >  1:2328265‑M1/1‑29 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                   <  1:2467669‑M1/49‑21 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                   <  1:422613‑M1/49‑21 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                >  1:2133740‑M1/1‑26 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                <  1:842085‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                <  1:1567208‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                >  1:2465482‑M1/1‑26 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                <  1:2483385‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                 <  1:958233‑M1/48‑23 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt              >  1:1257414‑M1/1‑23 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg               <  1:1007866‑M1/47‑25 (MQ=255)
                           cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtta            <  1:11610‑M1/49‑28 (MQ=255)
                            tCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa           >  1:618472‑M1/1‑21 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc         >  1:1368798‑M1/1‑19 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc         >  1:611040‑M1/1‑19 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc         >  1:2081128‑M1/1‑19 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc         <  1:1691709‑M1/49‑31 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac          <  1:2615829‑M1/48‑30 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa           <  1:1489709‑M1/47‑29 (MQ=255)
                                tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc       >  1:1505407‑M1/1‑17 (MQ=255)
                                tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc       >  1:1882840‑M1/1‑17 (MQ=255)
                                 ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc       >  1:1545342‑M1/1‑16 (MQ=255)
                                   aaCACGGAGAAGTGcctacactagcgggattctttttgttaaccccta     <  1:567858‑M1/48‑35 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc   <  1:2451209‑M1/49‑37 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctac    <  1:1502170‑M1/48‑36 (MQ=255)
                                      acGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc  >  1:514056‑M1/1‑11 (MQ=255)
                                            |                                          
CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGA  >  NC_000913/3090297‑3090383

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: