Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 20,703 +GCAAGAACCTGCAA intergenic (‑195/+112) insA1 ← / ← rpsT IS1 family protein InsA/30S ribosomal subunit protein S20

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 207031 (0.020)19 (0.420)
+GCAAGAACCTGCAA
14/458 1.8 95.3% intergenic (‑195/+112) insA1/rpsT IS1 family protein InsA/30S ribosomal subunit protein S20
?NC_000913 20704 = 1 (0.020)intergenic (‑196/+111) insA1/rpsT IS1 family protein InsA/30S ribosomal subunit protein S20

ATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCTGTG  >  NC_000913/20529‑20939
                                                                                                                                                                               |                                                                                                                                                                                                                                           
aatagagtggtcttaagtttaatacctcaatgtgtatcacaatatccacattctttgtagggggtctggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTgc                                                                                                                                                                  >  1:61710‑M2/176‑251 (MQ=255)
  tagagtggtcttaagtttaatacctcactgtgtatcacaatatccacattctttgtagggggtctggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTgcgc                                                                                                                                                                >  2:274043‑M2/174‑251 (MQ=255)
                gtttaatacctcaatgtgtatcacaatatccacattctttgtagggggtctggagattgcgtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTcaca                                                                                                                                                   >  1:43009‑M2/160‑250 (MQ=255)
                 tttaatacctcaatgtgtatctcaatatccacattctttgtagggggtctggagattgattagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGCCTATCTGGCa                                                                                                                                                                                                                             <  2:263753‑M2/17‑1 (MQ=255)
                 tttaatacctcaatgtgtatcacaatatccacattctttgtagggggtctggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCa                                                                                                                                                                                                                             >  1:263753‑M2/159‑175 (MQ=255)
                                       caatatccacattctttgtagggggtctggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTa                                                                                                                            >  2:439508‑M2/137‑250 (MQ=255)
                                         atatccacattctttgtagggggtctggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGc                                                                                                                         >  1:77048‑M2/135‑251 (MQ=255)
                                             ccacattctttgtagggggtctggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGt                                                                                                                     >  1:26828‑M2/131‑251 (MQ=255)
                                                                   ggagattgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGt                                                                                               >  1:160300‑M2/109‑251 (MQ=255)
                                                                       attgagtagatattcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAAt                                                                                                                                     >  1:365103‑M2/105‑209 (MQ=255)
                                                                       atggagtagatattcttgttcagaatttatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAAt                                                                                                                                     <  2:365103‑M2/105‑1 (MQ=255)
                                                                         tgagtagatattcttgttcagaatgtatcagcagatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgccagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAACGCAGCAGGCATAACAAAACCCGCTTGCGCGGGCTTTTTCCCAAAGCTTCCGCAAATTGGCGATTAAGACAGTTTGTTGCTCTGTGCAGTCCAGTTAGCCt                                                                                         >  1:469498‑M2/103‑251 (MQ=255)
                                                                                   ttcttgttcagaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGc                                                                               >  1:116440‑M2/93‑251 (MQ=255)
                                                                                             gaatgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGttt                                                                     >  1:226169‑M2/83‑251 (MQ=255)
                                                                                               atgtatcagccgatggttctacgattcttaagccacgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAACTTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACATGGAGCTTTGttttt                                                                   >  1:498852‑M2/81‑251 (MQ=255)
                                                                                                                     gCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTGGCGCCGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGGCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTagcagc                                              <  1:309992/250‑1 (MQ=255)
                                                                                                                                  cgaagagttcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCCGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACg                                 >  2:281322‑M2/46‑250 (MQ=255)
                                                                                                                                          tcagatagtacaacggcatgtctgcaagaacctgcaaCTTTTGACTATCTGGCAAACGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGc                        >  1:194869‑M2/38‑251 (MQ=255)
                                                                                                                                                agtacaacggcaggtcttcaagaacctgtaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCg                  <  2:290437‑M2/220‑1 (MQ=255)
                                                                                                                                                           atgactgcaagaaccttcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCttt       <  1:394088‑M2/231‑1 (MQ=255)
                                                                                                                                                                 gcaagaacctgcaaCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCtgtg  <  2:105388‑M2/236‑1 (MQ=255)
                                                                                                                                                                               |                                                                                                                                                                                                                                           
ATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCTGTG  >  NC_000913/20529‑20939

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