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	= 20703	1 (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

TCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCTGTGCAGCAGCTT  >  NC_000913/20542‑20948                                                                                                                                                                  |                                                                                                                                                                                                                                                      aatagagtggtcttaagtttaatacctcaatgtgtatcacaatatccacATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgc                                                                                                                                                             >  1:61710‑M1/23‑161 (MQ=255)   tagagtggtcttaagtttaatacctcactgtgtatcacaatatccACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgc                                                                                                                                                           >  2:274043‑M1/21‑159 (MQ=255)                 gtttaatacctcaatgtGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGCGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcaca                                                                                                                                              >  1:43009‑M1/7‑145 (MQ=255)                  tttaatacctcaatgTGTATCTCAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGATTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgcctatctggca                                                                                                                                                                                                                        <  2:263753‑M1/170‑32 (MQ=255)                  tttaatacctcaatgTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggca                                                                                                                                                                                                                        >  1:263753‑M1/6‑144 (MQ=255)                                        caATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgatta                                                                                                                       >  2:439508‑M1/1‑122 (MQ=255)                                          atatCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagc                                                                                                                    >  1:77048‑M1/1‑120 (MQ=255)                                              ccacattcttTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagt                                                                                                                >  1:26828‑M1/5‑116 (MQ=255)                                                                     ggAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggt                                                                                          >  1:160300‑M1/1‑94 (MQ=255)                                                                         atggagtAGATATTCTTGTTCAGAATTTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaat                                                                                                                                <  2:365103‑M1/206‑120 (MQ=255)                                                                         aTTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaat                                                                                                                                >  1:365103‑M1/1‑90 (MQ=255)                                                                           tGAGTAGATATTCTTGTTCAGAATGTATCAGCAGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgccagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaacgcagcaggcataacaaaacccgcttgcgcgggctttttcccaaagcttccgcaaattggcgattaagacagtttgttgctctgtgcagtccagttagcct                                                                                    >  1:469498‑M1/1‑88 (MQ=255)                                                                                     ttCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgc                                                                          >  1:116440‑M1/1‑78 (MQ=255)                                                                                               gAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgttt                                                                >  1:226169‑M1/1‑68 (MQ=255)                                                                                                 aTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaacttggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacatggagctttgttttt                                                              >  1:498852‑M1/1‑66 (MQ=255)                                                                                                         gCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTGGCGCCGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGGCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTagcagc                                                       <  1:309992/250‑1 (MQ=255)                                                                                                                                    cGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtccggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacg                            >  2:281322‑M1/1‑31 (MQ=255)                                                                                                                                            tCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaacggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgc                   >  1:194869‑M1/1‑23 (MQ=255)                                                                                                                                                  aGTACAACGGCAGGTCTtcaagaacctgtaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgcatttcg             <  2:290437‑M1/251‑235 (MQ=255)                                                                                                                                                             atgactgcaagaaccttcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgcatttcgttaaatgcttt  <  1:394088‑M1/247‑246 (MQ=255)                                                                                                                                                                  |                                                                                                                                                                                                                                                      TCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCTGTGCAGCAGCTT  >  NC_000913/20542‑20948

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.