New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 139420812 (0.380)4 (0.130) 4/482 2.1 25.0% coding (982/1614 nt) mppA murein tripeptide ABC transporter periplasmic binding protein
?NC_000913 1396044 = NA (NA)noncoding (1/1195 nt) IS5 repeat region

CCGAAAATAATGACCAGATTAAGCGCATTATTTGTAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCT  >  NC_000913/1395870‑1396281
                                                                                                                                                                              |                                                                                                                                                                                                                                             
ccGAAAATAATGACCAGATTAAGCGCATTATTTGTAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCTTAGAAAAGGGTTCATCATGAGtcat                                                                                                                                                                   >  1:214075/1‑251 (MQ=255)
ccGAAAATAATGACCAGATTAAGCGCATTATTTGTAAACAGTTCAGCGATCATTGTTTTAACAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTTCCGCGACCATTCATTCAGATCATCAATAGTCGGGGACGGAAGTAGAAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGACCATTTTTGTCATCTGGAGCCATAGAACAAGGGTCATGATGAGtcat                                                                                                                                                                   >  2:210800/1‑251 (MQ=38)
 cGAAAATAATGACCAGATTAAGCGCATTATTTGTAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGAGCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCAGGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatc                                                                                                                                                                  <  1:39140/251‑1 (MQ=255)
                           ttattTGTAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTc                                                                                                                                        >  1:288054/1‑251 (MQ=255)
                                  tAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTACGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtc                                                                                                                                                                     >  2:101279/1‑215 (MQ=255)
                                  tAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtc                                                                                                                                                                     <  1:101279/215‑1 (MQ=255)
                                     aCAGTTCAGAGATCATTGTTTTATCAGGCTCCTCCAGATACATGTCCTCATGCCGGAAACACCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTCGATAAGGAAGGTGCGAACCAGTCCCTGAACTGAGATCATGTTTGTCATCTGGAGCCATAGAACGGGGTTCCTCATGACTCATCAACTTACCTTCGCCGCCAGTGACTTCAGCAGTCAgc                                                                                                                              >  2:43805/1‑251 (MQ=25)
                                                                                                               ttttaacccgACAATTCATTCAGATCATCCATAGTCAGGGAAGGAAGTAGCAACATTAGATAAGGAAGGTGCGAAAAAGTCCCTGATTTGAGATCATGTTTGTCATATGGAGCAATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAACCATGGTGGAAGTc                                                    >  2:215761/9‑251 (MQ=17)
                                                                                                                        gACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCg                                           >  1:176955/1‑251 (MQ=255)
                                                                                                                                    aGATCATCAATAGTCAGGGAAGCAAGTAGCAACATTAGCTAAGGAAGGAGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGCGCCGTTTTACCCCa                                 >  2:282361/1‑249 (MQ=255)
                                                                                                                                              tAGTCAGGGAAGGAAGTAGCAACATTAGCTACGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATg                     >  1:256245/1‑251 (MQ=255)
                                                                                                                                                gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGc                   >  1:131390/1‑251 (MQ=40)
                                                                                                                                                     ggaaggaaGTAGCAACATTAGCTAAGGAAGGAGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGATTTACCCCAAGGCTGGTAATGGCCGGCg              >  1:199791/1‑251 (MQ=255)
                                                                                                                                                            agTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGCGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTAt       >  1:126331/1‑251 (MQ=25)
                                                                                                                                                               aGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGGCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGGAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATcc     <  1:215761/250‑1 (MQ=21)
                                                                                                                                                                 cAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCt  <  2:122785/251‑1 (MQ=33)
                                                                                                                                                                              |                                                                                                                                                                                                                                             
CCGAAAATAATGACCAGATTAAGCGCATTATTTGTAAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCT  >  NC_000913/1395870‑1396281

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