Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,255,321	IS1 (–) +8 bp	intergenic (‑135/‑27)	rihB ← / → yeiL	pyrimidine‑specific ribonucleoside hydrolase RihB/DNA‑binding transcriptional activator YeiL

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1978503 =	NA (NA)	14 (0.390)	12/486	1.5	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 2255328	0 (0.000)					intergenic (‑142/‑27)	rihB/yeiL	pyrimidine‑specific ribonucleoside hydrolase RihB/DNA‑binding transcriptional activator YeiL
*	?	NC_000913	= 1979270	NA (NA)	25 (0.690)	25/486	0.5	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	2255321 =	0 (0.000)					intergenic (‑135/‑34)	rihB/yeiL	pyrimidine‑specific ribonucleoside hydrolase RihB/DNA‑binding transcriptional activator YeiL

TGCCTAATAAATCTATTGCCGGATGTTTCGCCGCCATCATTATAGCAATAGCATCATCATGACCCGGATCACAATCCAGAATAATTTTTCTCTTTTCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATGGTCGCGTGTCGCTGATCGATTTCTTTGCCGCCCCCTGTTTTA  >  NC_000913/2255087‑2255562                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           gccatatcaatgattttctggtgcgtaccgggttgagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccagcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATg                                                                                                                                                                                                                                   <  2:40372‑M2/17‑1 (MQ=255)  gCCTAATAAATCTATTGCCGGATGTTTCGCCGCCATCATTATAGCAATAGCATCATCATGACCCGGATCACAATCCAGAATAATTTTTCTCTTTTCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgact                                                                                                                                                                                                                                  >  2:294104‑M2/1‑241 (MQ=255)          atgattttctggtgcgtaccgggttgagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAagc                                                                                                                                                                                                                          <  1:337683‑M2/26‑1 (MQ=255)                         gtaccgggttgagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTg                                                                                                                                                                                                           >  2:74937‑M2/211‑251 (MQ=255)                                    agaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGttt                                                                                                                                                                                                <  1:15866‑M2/52‑1 (MQ=255)                                             gtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGc                                                                                                                                                                                       >  1:134323‑M2/191‑251 (MQ=255)                                              caatagcaTCATCATGACCCGGATCACAATCCAGAATAATTTTTCTCTTTTCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgac                                                                                                                                                                                      <  1:324692‑M2/251‑55 (MQ=255)                                                agtgaactgcagttgccatgttttacggcagtgagagcagagctagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCttt                                                                                                                                                                                    >  1:233753‑M2/188‑251 (MQ=255)                                                agtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCttt                                                                                                                                                                                    <  1:224274‑M2/64‑1 (MQ=255)                                                 gtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCtttt                                                                                                                                                                                   >  1:143768‑M2/187‑251 (MQ=255)                                                     actgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggcgattttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAAc                                                                                                                                                                               >  1:336484‑M2/183‑251 (MQ=255)                                                               ccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTGAAGGATTGCTTTTTAACGGATTTTTCa                                                                                                                                                                     <  2:47297‑M2/79‑1 (MQ=255)                                                                 cGGATCACAATCCAGAATAATTTTTCTCTTTTCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc                                                                                                                                                                   >  2:294782‑M2/1‑178 (MQ=255)                                                                 atgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGc                                                                                                                                                                   >  1:76441‑M2/171‑251 (MQ=255)                                                                         cggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATAcgc                                                                                                                                                           <  1:95152‑M2/89‑1 (MQ=255)                                                                            ccAGAATAATTTTTCTCTTTTCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAATTATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgaga                                                                                                                                                         <  2:260562‑M2/250‑84 (MQ=255)                                                                                 agagcagagatagcgctgatgtccggcggtgcttttgccgttactcaccgccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATATTTCAGCCGATACGCGGCTGttt                                                                                                                                                   <  2:158799‑M2/97‑1 (MQ=255)                                                                                              tttCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccg                                                                                                                                         >  1:147324‑M2/1‑149 (MQ=255)                                                                                               gctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGcgcg                                                                                                                                     <  2:110649‑M2/111‑1 (MQ=255)                                                                                                            ggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATtgtg                                                                                                                        <  1:217864‑M2/124‑1 (MQ=255)                                                                                                                 ctGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaactgattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacggcagcgacttccgtcccagccgtgccaggtgct                                                                                                                     >  2:231222‑M2/1‑130 (MQ=255)                                                                                                                         aGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagatt                                                                                                           <  1:53804‑M2/251‑130 (MQ=255)                                                                                                                           acgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACa                                                                                                          <  2:312997‑M2/138‑1 (MQ=255)                                                                                                                              cGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggt                                                                                                      >  1:331905‑M2/1‑117 (MQ=255)                                                                                                                                          agtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGttt                                                                                          >  2:114793‑M2/98‑251 (MQ=255)                                                                                                                                                           gACCGGGCCGATGATTGAATCTGAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgagtttgagaacgacagcgactgccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtataacgc                                                                         >  2:243297‑M2/1‑88 (MQ=255)                                                                                                                                                                                    caacaGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttgtatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgataacgacagcgacgtccgtcccagccgtgccatgtgctgcctcagattcaggttatgccgctcaattctctgcgtatatcgcttgctgattacgtgcagctttccct                                                >  1:352553‑M2/1‑63 (MQ=255)                                                                                                                                                                                      gagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACgctagcta                                               >  2:350210‑M2/54‑250 (MQ=255)                                                                                                                                                                                        gcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATg                                            <  1:206068‑M2/200‑1 (MQ=255)                                                                                                                                                                                            tacgtaTGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgg                                        <  1:151418‑M2/251‑197 (MQ=255)                                                                                                                                                                                              caaaaacaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATGGTcgcg                                      >  2:150300‑M2/46‑251 (MQ=255)                                                                                                                                                                                                   acaccatcatacactaaatcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATGGTCGCGTGTCg                                 <  1:252633‑M2/211‑1 (MQ=255)                                                                                                                                                                                                                     tcagtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATGGTCGCGTGTCGCTGATCGATTTCTTTgcc               >  2:131056‑M2/23‑251 (MQ=255)                                                                                                                                                                                                                      aTAGCGGATAAACATCGTGAAAGAAATCCggtactgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgt              <  1:158799‑M2/251‑223 (MQ=255)                                                                                                                                                                                                                        gtaagttggcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGATTTACGCCACGCTAGCTAATGGTCGCGTGTCGCTGATCGATTTCTTTgccgcc            >  2:122942‑M2/20‑251 (MQ=255)                                                                                                                                                                                                                                tacATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgttccagccgtgccaggtgctgcctcagattcagggtatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatc    <  2:233549‑M2/250‑233 (MQ=255)                                                                                                                                                                                                                                gcagcatcaccAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATGGTCGCGTGTCGCTGATCGATTTCTTTGCCGCCCCCTGttt    <  2:134323‑M2/240‑1 (MQ=255)                                                                                                                                                                                                                                  cATCGTGAAAGAAATCCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatcac  <  1:170658‑M2/251‑235 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           TGCCTAATAAATCTATTGCCGGATGTTTCGCCGCCATCATTATAGCAATAGCATCATCATGACCCGGATCACAATCCAGAATAATTTTTCTCTTTTCCATTGTTTATTTCCTCTGTTTCCAGTTGCGTTATTTTTTCTACAGCAAAGAAAATTCGACCGGGCCGATGATTGAATCTTAACAACAGCGTACGTATGCTAAATATGAGAAATCTCATAGCGGATAAACATCGTGAAAGAAATCCACAATAATGATCTTAAGCAGCAATTGATGAGTGAATCCGCGTTTAAGGATTGCTTTTTAACGGATGTTTCAGCCGATACGCGGCTGTTTCATTTTTTAGCGCGTGACTACATTGTGCAGGAAGGGCAACAGCCGTCCTGGCTGTTTTACCTGACGCGAGGCCGCGCCAGGCTTTACGCCACGCTAGCTAATGGTCGCGTGTCGCTGATCGATTTCTTTGCCGCCCCCTGTTTTA  >  NC_000913/2255087‑2255562

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