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	2 (0.040)	26 (0.600) +GCAAGAACCTGCAA	20/450	0.9	93.2%	intergenic (‑195/+112)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20
	?	NC_000913	20704 =	2 (0.040)					intergenic (‑196/+111)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20

TCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCT  >  NC_000913/20530‑20932                                                                                                                                                                              |                                                                                                                                                                                                                                                    tcatACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGCGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCAGGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAACAACCCGCTTgcg                                                                                                                                                          >  2:454868/1‑251 (MQ=255)           aaaatagagtggtcttaagtttaatacctcaatgtgtatcacaatatccacatTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgctt                                                                                                                                                             >  1:305898‑M1/25‑163 (MQ=255)                       tcttaagtttaatacctcaatgtgtatcaCAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttt                                                                                                                                                 >  1:232909‑M1/13‑151 (MQ=255)                           aagtttaatacctcaatgtgtATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttaac                                                                                                                                             >  1:224553‑M1/9‑147 (MQ=255)                                  atacctcAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagctt                                                                                                                                      >  1:169582‑M1/2‑140 (MQ=255)                                   taccgCAATGGGTATCACAATATACACATTCTTTGTCGGGG‑GTCTGGCGTTTGAGTCGATATTCTTGTTCAGATTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACTAAGAGTTCAGATACTACCACGGCATGTCTGCAAGCACCTGCAACTTTTGACTATCTGGCCACCGGCAGTGTGTTCTCTACCGCATCACAAAAGCATCCGGCCCAAAAAATACCGCTTTTGCAGGCTTTTTCACAAAACTTc                                                                                                                                     >  2:31573/1‑251 (MQ=255)                                        cAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTCCAGAATGTATCAGCCGATGGTTCTACTATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagca                                                                                                                                 <  1:156224‑M1/250‑117 (MQ=255)                                         aaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaa                                                                                                                                >  1:348856‑M1/1‑133 (MQ=255)                                         aaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgt                                                                                                                                                                                                         <  2:147465‑M1/177‑45 (MQ=255)                                         aaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgt                                                                                                                                                                                                         >  1:147465‑M1/1‑133 (MQ=255)                                                  cacaATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGGTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgat                                                                                                                       <  2:232909‑M1/250‑127 (MQ=255)                                                          ccacattcttTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagt                                                                                                              <  1:480818‑M1/247‑136 (MQ=255)                                                                            ggtCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcagaaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagt                                                                                             <  1:315173‑M1/250‑153 (MQ=255)                                                                                  gggattTAGTAGATATACGTGGTCAGAATGTAGCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggtt                                                                                       <  2:408592‑M1/249‑159 (MQ=255)                                                                                                  tCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgagcagtctggttagccttatgacgtgc                                                                        >  2:321475‑M1/1‑77 (MQ=255)                                                                                                          aGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtt                                                               <  1:208804‑M1/251‑183 (MQ=255)                                                                                                            aaTGTATCAGCCGATGGTTCTACGTTTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttt                                                             <  1:239796‑M1/251‑185 (MQ=255)                                                                                                                tATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcaccaaagcagcaggcataaaa                                                                                                                                                                       >  1:361050‑M1/1‑63 (MQ=255)                                                                                                                tATCAGCCGATGGATCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTTCAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaa                                                                                                                                                                       <  2:361050‑M1/141‑79 (MQ=255)                                                                                                                         aTGGTTTTACGATTCTTAAGCCACGAAGAGTTCAGATAGTTCAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacc                                                <  2:270921‑M1/251‑198 (MQ=255)                                                                                                                                 aCGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcat                                                                                                                                                                           >  1:253620‑M1/1‑46 (MQ=255)                                                                                                                                 aCGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcat                                                                                                                                                                           <  2:253620‑M1/120‑75 (MQ=255)                                                                                                                                                  aaGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatc                       >  1:151795‑M1/1‑29 (MQ=255)                                                                                                                                                   agagTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaacagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgctctgtgcagtcaggtt                                                                                       >  2:403923‑M1/1‑28 (MQ=255)                                                                                                                                                   agagTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggtt                                                                                       <  1:403923‑M1/186‑159 (MQ=255)                                                                                                                                                              aGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgcatttc            <  2:151795‑M1/250‑234 (MQ=255)                                                                                                                                                               gTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcctcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattatgccagtttgttgatctgtgcagtctggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggaccacgatcggttgcatttcgt          >  2:5630‑M1/1‑16 (MQ=255)                                                                                                                                                                       gCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagtttgttgatctgtgcagtcaggttagccttatgacgtgcagctttgtttttgtggatcagacctttagcagcctgacggtccacgatcggttgcatttcgttaaatgct  >  1:194145‑M1/1‑8 (MQ=255)                                                                                                                                                                              |                                                                                                                                                                                                                                                    TCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAGCCTGACGGTCCACGATCGGTTGCATTTCGTTAAATGCT  >  NC_000913/20530‑20932

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.