Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 702,352 Δ21 bp coding (380‑400/1149 nt) nagA ← N‑acetylglucosamine‑6‑phosphate deacetylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 702352 702372 21 20 [0] [0] 20 nagA N‑acetylglucosamine‑6‑phosphate deacetylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 7023510 (0.000)20 (0.770) 17/296 0.4 100% coding (401/1149 nt) nagA N‑acetylglucosamine‑6‑phosphate deacetylase
?NC_000913 702373 = 0 (0.000)coding (379/1149 nt) nagA N‑acetylglucosamine‑6‑phosphate deacetylase

CTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCACGGACCTTCCAGATGCAGACCTAACGCCTGATTCGGATGTTTTGCCAGGTACTCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTT  >  NC_000913/702214‑702499
                                                                                                                                         |                                                                                                                                                    
cTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgatt                                                                                                                                         <  1:80147‑M1/151‑14 (MQ=255)
cTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgatt                                                                                                                                         <  1:47141‑M1/151‑14 (MQ=255)
                            ggTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcat                                                                                                              <  1:251948‑M1/150‑41 (MQ=255)
                                   ttGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgca                                                                                                      >  2:93242‑M1/1‑103 (MQ=255)
                                            aCGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtt                                                                                             <  1:376606‑M1/151‑58 (MQ=255)
                                                       ttCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctca                                                                                  >  1:272810‑M1/1‑83 (MQ=255)
                                                               ggAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggt                                                                          >  2:66715‑M1/1‑75 (MQ=255)
                                                                            cgcgGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcg                                                               <  1:303623‑M1/149‑88 (MQ=255)
                                                                            cgcgGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagc                                                                <  1:362282‑M1/148‑87 (MQ=255)
                                                                                       gCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTGTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatag                                                    <  1:184057‑M1/149‑99 (MQ=255)
                                                                                            cGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtac                                             <  1:77443‑M1/151‑106 (MQ=255)
                                                                                            cGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtac                                             <  1:165647‑M1/151‑106 (MQ=255)
                                                                                                caAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcc                                         <  1:63674‑M1/151‑110 (MQ=255)
                                                                                                         gATTATGGGTGCCTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttct                                 >  1:22087‑M1/1‑33 (MQ=255)
                                                                                                                     cTTTTTTTACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttctcattggctttctg                    <  1:162055‑M1/151‑131 (MQ=255)
                                                                                                                      tttttttACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttctcattggctttctgc                   <  1:66965‑M1/151‑132 (MQ=255)
                                                                                                                      tttttttACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttctcattggctttctgc                   <  1:98892‑M1/151‑132 (MQ=255)
                                                                                                                       ttttttACCAGATTCAGCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttctcattggctttctgca                  >  2:170772‑M1/1‑19 (MQ=255)
                                                                                                                                       gCCctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttctcattggctttctgcatgatttccagcgttt   >  1:344670‑M1/1‑3 (MQ=255)
                                                                                                                                        ccctaacgcctgattcggatgttttgccaggtactcgcgcataacgcgcacgccctgtttcatcagctcatcgctggtggtgataagcgtcggcagatagttagtacagcctgatttctcattggctttctgcatgatttccagcgtttc  <  2:240995‑M1/150‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
CTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGTGCCTTTTTTTACCAGATTCAGCCACGGACCTTCCAGATGCAGACCTAACGCCTGATTCGGATGTTTTGCCAGGTACTCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTT  >  NC_000913/702214‑702499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ 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.