Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,955,094	Δ1 bp :: IS186 (+) +6 bp :: Δ1 bp	coding (905‑910/3543 nt)	recB ←	exodeoxyribonuclease V subunit RecB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2514274 =	NA (NA)	14 (0.920)	13/302	0.1	100%	noncoding (2/1345 nt)	IS186	repeat region
	?	NC_000913	= 2955099	0 (0.000)					coding (905/3543 nt)	recB	exodeoxyribonuclease V subunit RecB
*	?	NC_000913	= 2515616	NA (NA)	12 (0.790)	10/302	0.3	100%	noncoding (1344/1345 nt)	IS186	repeat region
	?	NC_000913	2955094 =	0 (0.000)					coding (910/3543 nt)	recB	exodeoxyribonuclease V subunit RecB

ATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACC  >  NC_000913/2954859‑2955311                                                                                                                                                                                                                                            |                                                                                                                                                                                                                          acataatccagccatcgctggatttcccccccagaagtgccggatccgaaaagaaggactaactcgttgtggagaataacaaaaatggtcatctggagcttacaggtggccattcgtgggacagtatccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCtt                                                                                                                                                                                                            <  1:123745‑M2/16‑1 (MQ=255)                 gcgTACGGATCGCCGCTGCCAACACCTAACCGCTCTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGAGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgt                                                                                                                                                                                            <  1:299040‑M2/251‑27 (MQ=255)                            gccgcTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTTTGTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGGTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcc                                                                                                                                                                                 >  1:356229‑M2/1‑214 (MQ=255)                                   ncaACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGCCAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCccccataagcgctaacttaagggttgtggtattacgcctgatat                                                                                                                                                                          <  1:36694‑M2/250‑45 (MQ=255)                                                                       ggagaataacaaaaatggtcatctggagcttacaggtggccattcgtgggacagtatccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCGATCTTCTa                                                                                                                                                                                              >  1:296764‑M2/166‑195 (MQ=255)                                                                            cGACTTAACATGTCATCCCAACCCACTTCGGCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGAGCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtc                                                                                                                                   <  1:280324‑M2/249‑84 (MQ=255)                                                                                                              ccattcgtgggaaagtatacctgacagcctacaaaacgccattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGcc                                                                                               >  1:301417‑M2/127‑251 (MQ=255)                                                                                                                  gcgTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGAAGTTGATCGATCGCCTCAAACAGTTGATGTCGCGGGGTTTCCCcccataagcgctaacgtaagggttgtggtattacgcctgatattatttaacgtgccgattaattactctcacgataactggtcagcaattctggcccatattgttatgcccgaagaactggat                                                                                           >  1:116426‑M2/1‑128 (MQ=255)                                                                                                                        ttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGGGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCTCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgccagatatgatttaacgggccgatgaattactctcacgataactggtcagcaagtctggcccatattggtaagcccgaagaactggatacttc                                                                                      >  1:249866‑M2/1‑122 (MQ=255)                                                                                                                               gcTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGAGGATCGATCATGGTTCTGCCAGCAGTTAATCGATAGACTCAAACAGTGGATGTCGCGGGGTTTAACcccatcagagctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgacttactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgt                                                                               <  1:230860‑M2/250‑136 (MQ=255)                                                                                                                                        gcctacaaaacgcaattgaagaacgcgaggcatcggcttaacgaggcaccgaggcgtcgcagtcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGtgn                                                                                                                                                                                                       <  1:41890‑M2/21‑2 (MQ=255)                                                                                                                                                                 cgaggcatctgcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCGTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTc                                                                                                                                                                        <  1:199245‑M2/52‑1 (MQ=255)                                                                                                                                                                  gATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgatagctgttcagcaattctggcccatattggtaagcccgaagaactgtatacttcggcacgtaatgccgggtctctaacccgccgccgcgaaattcg                                            >  1:287613‑M2/1‑80 (MQ=255)                                                                                                                                                                      catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggatcaccccttaagttagcgcttatggcTTTCCCCACCGGCCTTCCTGCGATCATCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCTCTTGATAACTGTCTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCCTTTAGCCTGATTGCTACGGTTCAACTTGCGTCGATCACTCCCAGAAGATTCGATCAGCGCATc                                       >  1:351730‑M2/71‑251 (MQ=255)                                                                                                                                                                      aCCAGATCGCGGCTCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTAGCGGGGTTTCCCaccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaatgccggggctctaacccgccgccgcgaaattcgtgatg                                       <  1:297296‑M2/251‑176 (MQ=255)                                                                                                                                                                         cgtcttaacgaggcaccgatgcgtctcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCTATCTTCTAAGAATCTCTGTTAGAATTTTTCCATCGACTCCGGCTACTGATTTCTTTTTTTCTCTTCTTCTTCCCTGGCTCTGATCTTGTCGATCCCTTTATCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATcc                                      <  1:192711‑M2/182‑1 (MQ=255)                                                                                                                                                                          gATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccggtgaattactctcacgataactggtgagcaatgctggcccatagtggtaagcccgaggaactgtatacttcggcacgtaatgcgggggctctaacccgccgccgcgaaattcgtgatgctgc                                   >  1:334480‑M2/1‑72 (MQ=255)                                                                                                                                                                                    ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCACCCCGGACTTAGTGCGATCTTCTAAGa                                                                                                                                                                                           >  1:331539‑M2/57‑89 (MQ=255)                                                                                                                                                                                           ggTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctgtcccatattggtaagcccgaagaactggatacttcggcacgtaatgccggggctctaacccgccgccgcgaaattcgtgatgctgcaactctgctacgtctgg                  >  1:253882‑M2/1‑55 (MQ=255)                                                                                                                                                                                             gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATcc                                                                                                       <  1:141061‑M2/117‑1 (MQ=255)                                                                                                                                                                                              cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCAcc                              >  1:325554‑M2/47‑236 (MQ=255)                                                                                                                                                                                                  gcattcttcagatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCNgccnnnnngcnnnnnn                                                                                                                >  2:210224‑M2/43‑144 (MQ=255)                                                                                                                                                                                                  cAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaatgccggggctctaacccgccgccgcgaaattcgtgatgctgcaactctg                            <  1:311698‑M2/234‑187 (MQ=255)                                                                                                                                                                                                           agatggttcaacccttaagttagcgcttatgggTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTacn  <  1:210224‑M2/218‑2 (MQ=255)                                                                                                                                                                                                                     cTCAAACAGTGGATGTCGCGGGGTTTCCCcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaatgccggggctctaacccgccgccgc                                                    <  1:200486‑M2/191‑163 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                          ATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACC  >  NC_000913/2954859‑2955311

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