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

CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090297‑3090385
                                            |                                            
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                         <  1:5064057/49‑1 (MQ=255)
ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        >  1:3236815‑M1/1‑49 (MQ=255)
 cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       <  1:670143‑M1/50‑3 (MQ=255)
 cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:1874172‑M1/1‑48 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcct                                      <  1:3319146‑M1/50‑4 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcct                                      <  1:1598302‑M1/50‑4 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcct                                      <  1:576084‑M1/50‑4 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:4157696‑M1/1‑47 (MQ=255)
  aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:4731216‑M1/1‑47 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:4562158‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:2830043‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:5460466‑M1/1‑46 (MQ=255)
   aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:6103913‑M1/1‑46 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:2425805‑M1/50‑6 (MQ=255)
    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:8085372‑M1/50‑6 (MQ=255)
     aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                   <  1:5605494‑M1/50‑7 (MQ=255)
     aCAGCTTTCGCTACGTTGCTGGCTAGTTTTAACACGGAGAAGTGcctaca                                   <  1:15184‑M1/50‑7 (MQ=255)
       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                 <  1:6219533‑M1/50‑9 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               <  1:6006714‑M1/50‑11 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               <  1:5399937‑M1/50‑11 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               >  1:175321‑M1/1‑40 (MQ=255)
         cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               >  1:4994505‑M1/1‑40 (MQ=255)
          tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:1084698‑M1/1‑39 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            >  1:5625642‑M1/1‑37 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                             <  1:6390635‑M1/49‑13 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                             <  1:639270‑M1/49‑13 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                             <  1:7184706‑M1/49‑13 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:4821411‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:7809496‑M1/48‑12 (MQ=255)
            tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:1246469‑M1/48‑12 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            >  1:7823031‑M1/1‑36 (MQ=255)
             cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            <  1:2320585‑M1/49‑14 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          <  1:5530969‑M1/50‑16 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          <  1:2711935‑M1/50‑16 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:7727033‑M1/1‑35 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          <  1:2878478‑M1/50‑16 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          <  1:3302978‑M1/50‑16 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                           <  1:7877881‑M1/49‑15 (MQ=255)
              gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                           <  1:1753887‑M1/49‑15 (MQ=255)
               cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                         >  1:2474682‑M1/1‑34 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:4086194‑M1/1‑33 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:3114221‑M1/1‑33 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:4382924‑M1/1‑33 (MQ=255)
                tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:5227508‑M1/1‑33 (MQ=255)
                 aCGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                       >  1:5762393‑M1/1‑32 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                       <  1:4549975‑M1/49‑19 (MQ=255)
                  cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        <  1:2215277‑M1/48‑18 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:456338‑M1/1‑30 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:2604726‑M1/1‑30 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:5227348‑M1/1‑30 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:5966887‑M1/50‑21 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:6102010‑M1/50‑21 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:6287128‑M1/1‑30 (MQ=255)
                   gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:4621026‑M1/1‑30 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                    <  1:817999‑M1/50‑22 (MQ=255)
                    ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                    <  1:5251919‑M1/50‑22 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:6689548‑M1/1‑28 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:6222029‑M1/1‑28 (MQ=255)
                     tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:4253714‑M1/1‑28 (MQ=255)
                      gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   <  1:698698‑M1/49‑23 (MQ=255)
                        tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:6729132‑M1/1‑25 (MQ=255)
                         ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:2861306‑M1/1‑24 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtta              <  1:7622299‑M1/50‑28 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:6067445‑M1/49‑27 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:1614138‑M1/49‑27 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               >  1:1373090‑M1/1‑23 (MQ=255)
                          gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                <  1:3898821‑M1/48‑26 (MQ=255)
                           cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa             >  1:2453246‑M1/1‑22 (MQ=255)
                           cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa             <  1:3703153‑M1/50‑29 (MQ=255)
                            tCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac            >  1:7324119‑M1/1‑21 (MQ=255)
                            tCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac            >  1:1114458‑M1/1‑21 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:2541761‑M1/1‑19 (MQ=255)
                              gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:5293890‑M1/1‑19 (MQ=255)
                               ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         <  1:4371301‑M1/50‑33 (MQ=255)
                                tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:661104‑M1/1‑17 (MQ=255)
                                tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:4644164‑M1/1‑17 (MQ=255)
                                tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:3424646‑M1/1‑17 (MQ=255)
                                 ttAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta       >  1:3916501‑M1/1‑16 (MQ=255)
                                   aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     >  1:7579605‑M1/1‑14 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    <  1:4885671‑M1/50‑38 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    <  1:1209052‑M1/50‑38 (MQ=255)
                                    acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     <  1:664080‑M1/49‑37 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc   <  1:6938167‑M1/50‑39 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc   <  1:510879‑M1/50‑39 (MQ=255)
                                     cacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    >  1:100737‑M1/1‑12 (MQ=255)
                                      acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  <  1:598852‑M1/50‑40 (MQ=255)
                                      acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  <  1:256286‑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 < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: