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	0 (0.000)	26 (0.440) +GCAAGAACCTGCAA	22/462	1.1	100%	intergenic (‑195/+112)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20
	?	NC_000913	20704 =	0 (0.000)					intergenic (‑196/+111)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20

AGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCTGTGCAGCAGCTTTGTCG  >  NC_000913/20544‑20953                                                                                                                                                                |                                                                                                                                                                                                                                                           tagagtggtcttaagtttaatacctcaatgtgtatcacaatatccACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgc                                                                                                                                                                >  1:505400‑M1/21‑159 (MQ=255)          cttaagtttaatacctcaatgtgtatcACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggcttttt                                                                                                                                                       >  2:1153‑M1/12‑150 (MQ=255)             aagtttaatacctcaatgtgtATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctcgcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttca                                                                                                                                                     <  2:602418‑M1/242‑104 (MQ=255)                  taatacctcaaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaag                                                                                                                                                <  2:23199‑M1/247‑109 (MQ=255)                      acctCAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttca                                                                                                                                           <  2:40453‑M1/251‑114 (MQ=255)                       cctCAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGTAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgggttctctcacgcatcacaaaagcagcaggcagaaaaaaacccgcttgcgcgggctttttcacaaagcttcag                                                                                                                                          <  2:395091‑M1/251‑115 (MQ=255)                          cAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTGCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggcttttgcacaaagcttcagcaa                                                                                                                                       <  1:536548‑M1/251‑118 (MQ=255)                            aTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaat                                                                                                                                     >  1:125340‑M1/1‑132 (MQ=255)                             tgtgTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaatt                                                                                                                                    >  1:67‑M1/1‑131 (MQ=255)                                              acattcTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagttt                                                                                                                   >  2:486786‑M1/3‑114 (MQ=255)                                                atTCTTTGTAGGGG‑GTCTGGAGATTGAGTCGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgt                                                                                                                 >  1:191769‑M1/1‑112 (MQ=255)                                                atTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgt                                                                                                                 <  1:464459‑M1/251‑140 (MQ=255)                                                       gggagggggtctggagaTTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaatcttcagcaaattggcgattaagccagtttgttgatct                                                                                                           <  2:572054‑M1/243‑146 (MQ=255)                                                          agggggtatggagattgAGTAGATATTCTTGTTCAGCATGTATCAGCCGCTGGTTCTACGATTCTTAAGACACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagccaattggcgattaagccagtttgttgatctgtg                                                                                                        >  2:379828‑M1/9‑103 (MQ=255)                                                                              aGATATTCTTGTTCATATTGTATCAGCCTATGGTTTTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatga                                                                                    <  1:173949‑M1/251‑169 (MQ=255)                                                                                              aaTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttt                                                                    >  2:298067‑M1/1‑67 (MQ=255)                                                                                               aTGTATCAGACGATGGTTCTACGATTCTTAAGCCACGACGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggccgtgtgttctctcacgcatcacaaacgaagcaggcataaaaaaaaccgcttgcgcgggctttttcacaaagattcagcaaattggcgattaagccagtttgttgatctgtgcagtcagggtagccttatgacgtgcagctttgttttt                                                                   >  2:337244‑M1/1‑66 (MQ=255)                                                                                                 gTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgt                                                                 >  1:79437‑M1/1‑64 (MQ=255)                                                                                                      aGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatc                                                            <  2:505400‑M1/251‑193 (MQ=255)                                                                                                                        tCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctg                                           >  1:214221‑M1/1‑41 (MQ=255)                                                                                                                                cACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagcacctgcaacttttgactatctggccaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcg                                                                                                                                                                 >  2:65049‑M1/1‑33 (MQ=255)                                                                                                                                cACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcg                                                                                                                                                                 <  1:65049‑M1/124‑92 (MQ=255)                                                                                                                                        gTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagccgtacggtccacgatcggtt                          <  1:501633‑M1/251‑227 (MQ=255)                                                                                                                                                    cAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgaagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccccgatcggttgcatttcgtta               >  2:503928‑M1/1‑13 (MQ=255)                                                                                                                                                           aTGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgcgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgcatttcgttaaatgcttt       >  1:366458‑M1/1‑6 (MQ=255)                                                                                                                                                                tgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcataacaaaagcagcaggcataaaaaacaccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggatagacttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgcatttcgttaaatgctttctgtg  >  2:102395‑M1/1‑1 (MQ=255)                                                                                                                                                                |                                                                                                                                                                                                                                                           AGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCTTTCTGTGCAGCAGCTTTGTCG  >  NC_000913/20544‑20953

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