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

CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090297‑3090385
                                            |                                            
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                         <  1:2425232/49‑1 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        <  1:926128‑M1/50‑2 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        <  1:5415023‑M1/50‑2 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        >  1:3747869‑M1/1‑49 (MQ=255)
 cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:5030799‑M1/1‑48 (MQ=255)
 cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:1733457‑M1/1‑48 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcct                                      <  1:796726‑M1/50‑4 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:2674567‑M1/1‑47 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:4383864‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:2783481‑M1/1‑46 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:5834318‑M1/50‑6 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:2343455‑M1/50‑6 (MQ=255)
     aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                   <  1:3595367‑M1/50‑7 (MQ=255)
      cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  <  1:4731040‑M1/50‑8 (MQ=255)
      cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  <  1:3634394‑M1/50‑8 (MQ=255)
       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                 <  1:1246910‑M1/50‑9 (MQ=255)
       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:1034694‑M1/1‑42 (MQ=255)
       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:3470431‑M1/1‑42 (MQ=255)
        gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:5915653‑M1/1‑41 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               <  1:1385105‑M1/50‑11 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               <  1:5908989‑M1/50‑11 (MQ=255)
          tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:2207317‑M1/1‑39 (MQ=255)
          tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:978053‑M1/1‑39 (MQ=255)
          tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:5254195‑M1/1‑39 (MQ=255)
          tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:2230594‑M1/1‑39 (MQ=255)
           ttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                             >  1:460027‑M1/1‑38 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            >  1:2692002‑M1/1‑37 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:1972019‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:1411023‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:4804459‑M1/48‑12 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacgctagcgg                            <  1:2451852‑M1/49‑14 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            <  1:1682409‑M1/49‑14 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            <  1:3893369‑M1/49‑14 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            <  1:3362832‑M1/49‑14 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          <  1:1974423‑M1/50‑16 (MQ=255)
               cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                         >  1:4394151‑M1/1‑34 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:3428472‑M1/1‑33 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                      >  1:4187777‑M1/1‑31 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                      >  1:1527202‑M1/1‑31 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                       <  1:5768875‑M1/49‑19 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                       <  1:4729531‑M1/49‑19 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        <  1:3983132‑M1/48‑18 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:734535‑M1/1‑30 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:5009936‑M1/50‑21 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:5235027‑M1/50‑21 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:2556025‑M1/50‑21 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                    <  1:1739413‑M1/50‑22 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                    <  1:6253758‑M1/50‑22 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:4816765‑M1/1‑28 (MQ=255)
                      gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   <  1:5363093‑M1/49‑23 (MQ=255)
                      gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   <  1:423798‑M1/49‑23 (MQ=255)
                      gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   <  1:220968‑M1/49‑23 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                  <  1:6285470‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                  <  1:1643438‑M1/49‑24 (MQ=255)
                       cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                  <  1:2871145‑M1/49‑24 (MQ=255)
                        tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:38468‑M1/1‑25 (MQ=255)
                        tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:3824494‑M1/1‑25 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:3878406‑M1/49‑27 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:5376256‑M1/49‑27 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                <  1:1461411‑M1/48‑26 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                <  1:208000‑M1/48‑26 (MQ=255)
                           cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa             >  1:5135044‑M1/1‑22 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:2680535‑M1/1‑19 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac            <  1:2298897‑M1/48‑30 (MQ=255)
                               ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:1467982‑M1/1‑18 (MQ=255)
                               ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:2450469‑M1/1‑18 (MQ=255)
                                 ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:2025152‑M1/49‑34 (MQ=255)
                                 ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:3130672‑M1/49‑34 (MQ=255)
                                  tAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta       <  1:695197‑M1/49‑35 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    <  1:2322710‑M1/50‑38 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     <  1:4537008‑M1/49‑37 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc   <  1:1552206‑M1/50‑39 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    >  1:1527019‑M1/1‑12 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    >  1:787540‑M1/1‑12 (MQ=255)
                                      acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  <  1:2910859‑M1/50‑40 (MQ=255)
                                            |                                            
CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090297‑3090385

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: