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

CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACT  >  NC_000913/3090298‑3090381
                                           |                                        
cAAAACAGCTTTCGCTACTTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                    <  1:1743047‑M1/49‑2 (MQ=255)
   aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                 <  1:2113948‑M1/49‑5 (MQ=255)
   aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                 <  1:1985117‑M1/49‑5 (MQ=255)
   aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                 <  1:957286‑M1/49‑5 (MQ=255)
     cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                               <  1:1039729‑M1/49‑7 (MQ=255)
     cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                               <  1:699752‑M1/49‑7 (MQ=255)
      aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                              >  1:2176‑M1/1‑42 (MQ=255)
        cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                            <  1:992565‑M1/49‑10 (MQ=255)
        cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                              >  1:1570082‑M1/1‑40 (MQ=255)
          ttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                          >  1:1735878‑M1/1‑38 (MQ=255)
           tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                           <  1:493121‑M1/47‑11 (MQ=255)
             gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                        >  1:668432‑M1/1‑35 (MQ=255)
             gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                        >  1:597282‑M1/1‑35 (MQ=255)
                 ccttgcTGGCTCTTTTTAACACGGAGAAGTGcctacactagcgggatt                    <  1:1361757‑M1/46‑18 (MQ=255)
                 cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                   <  1:733179‑M1/49‑19 (MQ=255)
                 cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                   <  1:950029‑M1/49‑19 (MQ=255)
                 cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                   <  1:1076035‑M1/49‑19 (MQ=255)
                 cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                    <  1:1296617‑M1/48‑18 (MQ=255)
                 cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                     <  1:210612‑M1/47‑17 (MQ=255)
                  gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                  >  1:647540‑M1/1‑30 (MQ=255)
                  gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                  <  1:719889‑M1/49‑20 (MQ=255)
                      cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt              <  1:1409014‑M1/49‑24 (MQ=255)
                       tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg             >  1:563547‑M1/1‑25 (MQ=255)
                         gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt           <  1:928520‑M1/49‑27 (MQ=255)
                          cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtta          <  1:141120‑M1/49‑28 (MQ=255)
                          cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt           >  1:442718‑M1/1‑22 (MQ=255)
                           tCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa         >  1:59552‑M1/1‑21 (MQ=255)
                             gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac        <  1:1615100‑M1/48‑30 (MQ=255)
                                ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc     <  1:528218‑M1/48‑33 (MQ=255)
                                ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc     <  1:279183‑M1/48‑33 (MQ=255)
                                ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc     <  1:878560‑M1/48‑33 (MQ=255)
                                 tAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct    <  1:478413‑M1/48‑34 (MQ=255)
                                   acacGGAGAAGTGcctacactagcgggattctttttgttaacccctac  <  1:1400220‑M1/48‑36 (MQ=255)
                                           |                                        
CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACT  >  NC_000913/3090298‑3090381

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: