Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 389,898 Δ9,908 bp yaiT[yaiZ] yaiT, insF‑2, insE‑2, yaiT, iprA, ampH, sbmA, yaiW, yaiY, [yaiZ]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 389898 399805 9908 48 [0] [0] 51 yaiT–[yaiZ] yaiT,insF‑2,insE‑2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,[yaiZ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3898970 (0.000)48 (1.580) 42/98 0.0 100% intergenic (‑170/‑354) hemB/yaiT porphobilinogen synthase/putative autotransporter YaiT
?NC_000913 399806 = 0 (0.000)intergenic (+1/+23) yaiZ/ddlA DUF2754 domain‑containing protein YaiZ/D‑alanine‑‑D‑alanine ligase A

GGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  NC_000913/399758‑399854
                                                |                                                
aatagaaacagcaataattttatattcactgaaaatatttttaatcttg                                                  <  1:499936‑M2/1‑1 (MQ=255)
 atagaaacagcaataattttatattcactgaaaatatttttaatcttGaa                                                >  1:2106533‑M2/48‑50 (MQ=255)
 atagaaacagcaataattttatattcactgaaaatatttttaatcttGa                                                 <  1:2607889‑M2/2‑1 (MQ=255)
 atagaaacagcaataattttatattcactgaaaatatttttaatcttGa                                                 <  1:2255480‑M2/2‑1 (MQ=255)
 atagaaacagcaataattttatattcactgaaaatatttttaatcttGa                                                 <  1:2063767‑M2/2‑1 (MQ=255)
  tagaaacagcaataattttatattcactgaaaatatttttaatcttGaa                                                >  1:305744‑M2/47‑49 (MQ=255)
   agaaacagcaataattttatattcactgaaaatatttttaatcttGAata                                              >  1:1184467‑M2/46‑50 (MQ=255)
    gaaacagcaataattttatattcactgaaaatatttttaatcttGAata                                              <  1:2644294‑M2/5‑1 (MQ=255)
     aaacagcaataattttatattcactgaaaatatttttaatcttGAAtatt                                            <  1:746279‑M2/7‑1 (MQ=255)
     aaacagcaataattttatattcactgaaaatatttttaatcttGAAtatt                                            <  1:850569‑M2/7‑1 (MQ=255)
      aacagcaataattttatattcactgaaaatatttttaatcttGAAtatta                                           >  1:427265‑M2/43‑50 (MQ=255)
         agcaataattttatattcactgaaaatatttttaatcttGAATATTATCg                                        <  1:366025‑M2/11‑1 (MQ=255)
         agcaataattttatattcactgaaaatatttttaatcttGAATATTATCg                                        <  1:2439819‑M2/11‑1 (MQ=255)
         agcaataattttatattcactgaaaatatttttaatcttGAATATTATCg                                        >  1:2211946‑M2/40‑50 (MQ=255)
          gcaataattttatattcactgaaaatatttttaatcttGAATATTATCgg                                       >  1:2036241‑M2/39‑50 (MQ=255)
            aataattttatattcactgaaaatatttttaatcttGAATATTATCGGGc                                     <  1:538114‑M2/14‑1 (MQ=255)
             ataattttatattcactgaaaatatttttaatcttGAATATTATCGGGCg                                    <  1:942825‑M2/15‑1 (MQ=255)
             ataattttatattcactgaaaatatttttaatcttGAATATTATCGGGCg                                    >  1:543475‑M2/36‑50 (MQ=255)
             ataattttatattcactgaaaatatttttaatcttGAATATTATCGGGCg                                    >  1:270801‑M2/36‑50 (MQ=255)
                attttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTa                                 <  1:379925‑M2/18‑1 (MQ=255)
                 ttttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTa                                 >  1:2539647‑M2/32‑49 (MQ=255)
                  tttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTaa                                >  1:888088‑M2/31‑49 (MQ=255)
                   ttatattcactgaaaatatttttaatcttGAATATTGTCGGGCGTTaaaa                              >  1:1044471‑M2/30‑50 (MQ=255)
                   ttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTaaa                               <  1:684352‑M2/20‑1 (MQ=255)
                     atattcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAata                            >  1:1279329‑M2/28‑50 (MQ=255)
                     atattcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAat                             <  1:582494‑M2/22‑1 (MQ=255)
                      tattcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAatat                           >  1:520838‑M2/27‑50 (MQ=255)
                        ttcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAATAtt                          <  1:1538132‑M2/25‑1 (MQ=255)
                         tcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTa                         >  1:2558072‑M2/24‑49 (MQ=255)
                           actgaaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTACAt                      >  1:1597887‑M2/22‑50 (MQ=255)
                              gaaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTACATtgt                   >  1:1353864‑M2/19‑50 (MQ=255)
                               aaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTACATtgtg                  >  1:1179766‑M2/18‑50 (MQ=255)
                                  atatttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGtt               <  1:1064033‑M2/36‑1 (MQ=255)
                                  atatttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGtt               >  1:854411‑M2/15‑50 (MQ=255)
                                   tatttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGttt              <  1:2303259‑M2/37‑1 (MQ=255)
                                    atttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt             >  1:1101872‑M2/13‑50 (MQ=255)
                                     tttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt             >  1:964065‑M2/12‑49 (MQ=255)
                                     tttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt             <  1:1650184‑M2/38‑1 (MQ=255)
                                     tttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTc            <  1:1647042‑M2/39‑1 (MQ=255)
                                       tttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa          >  1:734475‑M2/10‑50 (MQ=255)
                                       tttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCa           <  1:654456‑M2/40‑1 (MQ=255)
                                        ttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa          <  1:2456554‑M2/41‑1 (MQ=255)
                                         taatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt         <  1:1493844‑M2/42‑1 (MQ=255)
                                         taatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt         <  1:1061396‑M2/42‑1 (MQ=255)
                                          aatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGc       >  1:2622717‑M2/7‑50 (MQ=255)
                                           atcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCg      <  1:2137610‑M2/45‑1 (MQ=255)
                                            tcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGt     >  1:2050696‑M2/5‑50 (MQ=255)
                                             cttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGt     >  1:2268651‑M2/4‑49 (MQ=255)
                                                gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGtt    <  1:1973028/47‑1 (MQ=255)
                                                gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGt     <  1:1499404/46‑1 (MQ=255)
                                                gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAt  <  1:1601107/49‑1 (MQ=255)
                                                |                                                
GGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  NC_000913/399758‑399854

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: