Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 92,854 (GCG)4→3 coding (1442‑1444/1767 nt) ftsI → peptidoglycan DD‑transpeptidase FtsI

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 92844 92846 3 26 [0] [0] 24 ftsI peptidoglycan DD‑transpeptidase FtsI

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 92847 =0 (0.000)18 (0.620) 16/76 0.3 100% coding (1435/1767 nt) ftsI peptidoglycan DD‑transpeptidase FtsI
?NC_000913 = 92843 0 (0.000)coding (1431/1767 nt) ftsI peptidoglycan DD‑transpeptidase FtsI

TGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAA  >  NC_000913/92801‑92895
                                              |                                                
cattgtccgcactgtggtgcatatgatggaaagcgtggcgctaccaggc                                                <  1:772836‑M1/3‑1 (MQ=255)
 attgtccgcactgtggtgcatatgatggaaagcgtggcgctaccaggcg                                               <  1:1719647‑M1/4‑1 (MQ=255)
    gtccgcactgtggtgcatatgatggaaagcgtggcgctaccaggcggcg                                            <  1:73368‑M1/7‑1 (MQ=255)
    gtccgcactgtggtgcatatgatggaaagcgtggcgctaccaggcggc                                             >  1:2280251‑M1/43‑48 (MQ=255)
     tccgcactgtggtgcatatgatggaaagcgtggcgctaccaggcggcgg                                           <  1:694071‑M1/8‑1 (MQ=255)
        gcactgtggtgcatatgatggaaagcgtggcgctaccaGGCGGCGGCGt                                        >  1:1015323‑M1/39‑49 (MQ=255)
         cactgtggtgcatatgatggaaagcgtggcgctaccaGGCGGCGGCGTg                                       >  1:1751235‑M1/38‑49 (MQ=255)
           ctgtggtgcatatgatggaaagcgtggcgctaccaGGCGGCGGCGTGaa                                     >  1:224564‑M1/36‑49 (MQ=255)
                 tgcatatgatggaaagcgtggcgctaccaGGCGGCGGCGTGAaggcgga                               >  1:477556‑M1/30‑48 (MQ=255)
                   catatgatggaaagcgtggcgctaccaGGCGGCGGCGTGAAGGCGGCGa                             <  1:312425‑M1/22‑1 (MQ=255)
                      atgatggaaagcgtggcgctaccaGGCGGCGGCGTGAAGGCGGCGAtt                           <  1:2017045‑M1/24‑1 (MQ=255)
                       tgatggaaagcgtggcgctaccaGGCGGCGGCGTGAAGGCGGCGATTa                          <  1:1969018‑M1/25‑1 (MQ=255)
                         atggaaagcgtggcgctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAg                       >  1:2251084‑M1/22‑49 (MQ=255)
                            gaaagcgtggcgctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCt                    >  1:1016227‑M1/19‑49 (MQ=255)
                               agcgtggcgctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatc                 >  1:2784301‑M1/16‑49 (MQ=255)
                                   tggcgctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgtat             >  1:2323216‑M1/12‑49 (MQ=255)
                                     gcgctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgtatcg           <  1:1692800‑M1/40‑1 (MQ=255)
                                      cgctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGc          <  1:1005464‑M1/41‑1 (MQ=255)
                                       gctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGc          <  1:1518113‑M1/41‑1 (MQ=255)
                                        ctaccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCa        >  1:2848443‑M1/7‑49 (MQ=255)
                                           ccaGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTa     >  1:2738334‑M1/4‑49 (MQ=255)
                                             aGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa   >  1:2622092‑M1/2‑49 (MQ=255)
                                             aGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa   >  1:2111649‑M1/2‑49 (MQ=255)
                                              GGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaaa  <  1:1171445‑M1/49‑1 (MQ=255)
                                              |                                                
TGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAA  >  NC_000913/92801‑92895

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.