Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,825,796	Δ4,947 bp		srlA–[gutQ]	srlA, srlE, srlB, srlD, gutM, srlR, [gutQ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2825796	2830742	4947	31 [0]	[0] 32	srlA–[gutQ]	srlA,srlE,srlB,srlD,gutM,srlR,[gutQ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2825795	0 (0.000)	31 (0.430)	31/492	0.8	100%	intergenic (‑219/‑37)	mltB/srlA	membrane‑bound lytic murein transglycosylase B/sorbitol‑specific PTS enzyme IIC2 component
	?	NC_000913	2830743 =	0 (0.000)					coding (931/966 nt)	gutQ	D‑arabinose 5‑phosphate isomerase GutQ

GGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGCAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTAACAACGGGCACCGCCGCCACTTCTGGCGTCGGCAACGTCACCTCAGGAAAAGCAAAATG  >  NC_000913/2830504‑2830993                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           gcataacgacgcttgaacatgaggggtccatttaacagattcaaccaggggcaagtatggtaaagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataacgcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGa                                                                                                                                                                                                                                                  <  2:145087‑M2/11‑1 (MQ=255)                  catgaggggtccatttaacagattcaaccaggggcaagtatggtaaagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataaatcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattagtttgccgtAACCTGCAGGATTTCTATCAGGCCGGGat                                                                                                                                                                                                                                >  2:345822‑M2/223‑251 (MQ=255)                      aggggtccatttaacagattcaaccaggggcaagtatggtaaagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGattatt                                                                                                                                                                                                                            >  1:295033‑M2/219‑251 (MQ=255)                        gggtccatttaacagattcaaccaggggcaagtatggtaaagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACNTGCAGGATTTCTATCAGGCCGGGATTAttt                                                                                                                                                                                                                           >  2:284180‑M2/217‑250 (MQ=255)                        gggtccatttaacagattcaaccaggggcaagtatggtaaagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTa                                                                                                                                                                                                                          <  1:318932‑M2/35‑1 (MQ=255)                                 taacagatgcaaccaggggcaagtatggtaaagcatcccgccccgcacaaggaagcggttgtcactgcccgattcggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCa                                                                                                                                                                                                                  <  1:127283‑M2/43‑1 (MQ=255)                                                            gtaaagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGACGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCgatgacg                                                                                                                                                                                      >  2:139320‑M2/181‑248 (MQ=255)                                                               aagcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGtt                                                                                                                                                                                   <  1:405663‑M2/74‑1 (MQ=255)                                                                agcatcacgccccgcacaaggaagcggaagtcactgcccgatacggactttacataactcaactcattaccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTCTTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTg                                                                                                                                                                                  >  2:575735‑M2/177‑251 (MQ=255)                                                                 gcatcacgccccgcacaaggaagcggtagtcactgcccgatacggactttacataactcaactaattcccctcgcgatccttttattcaacctttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggtagccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCTCCAGCCGATGCAGGTTgg                                                                                                                                                                                 >  2:168737‑M2/176‑251 (MQ=255)                                                                                       gcggtagtccctgcccgatacggactttacataactcaactcattcccctcgctatccttttatgcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattagtttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCcg                                                                                                                                                           <  2:67243‑M2/98‑1 (MQ=255)                                                                                                       gatacggactttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAg                                                                                                                                           <  1:111417‑M2/114‑1 (MQ=255)                                                                                                               ctttacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGtgatg                                                                                                                                   >  1:255278‑M2/130‑251 (MQ=255)                                                                                                                  tacataactcaactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCTCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGAtt                                                                                                                                >  1:427364‑M2/127‑251 (MQ=255)                                                                                                                            aactcattcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGt                                                                                                                      <  1:370266‑M2/135‑1 (MQ=255)                                                                                                                                 attcccctcgctatccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTg                                                                                                                 >  1:576938‑M2/112‑251 (MQ=255)                                                                                                                                               ccttttattcaaactttcaaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCAc                                                                                                   >  1:115832‑M2/98‑251 (MQ=255)                                                                                                                                                  tttattcaacatttcaaattaaaatatttatctttcatgttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtATTCTGCAGGATTTCTCTCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGc                                                                                                 <  2:399623‑M2/156‑1 (MQ=255)                                                                                                                                                    tattcaatctttcaaattaaaatatttatctttcattttgcgatcaaaataaaccttttaaatctttcaatcttattagattagtttgccttAACCGGCAGGATTTCTTGCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACTAATAGTTTCACGCTg                                                                                               <  2:350747‑M2/158‑1 (MQ=255)                                                                                                                                                                caaattaaaatatttatctttcattttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGt                                                                                   <  1:176121‑M2/170‑1 (MQ=255)                                                                                                                                                                     taaaatatttatctttcattttgcgatcaacataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGAGTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGCGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTc                                                                             >  1:522496‑M2/76‑251 (MQ=255)                                                                                                                                                                                        tttgcgatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTaac                                                           >  2:168281‑M2/57‑250 (MQ=255)                                                                                                                                                                                             gatcaaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTAACAACggg                                                     <  1:449576‑M2/200‑1 (MQ=255)                                                                                                                                                                                                caaaataacacttttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTAACAACGGGCAc                                                  >  1:127176‑M2/49‑251 (MQ=255)                                                                                                                                                                                                            tttaaatctttcaatctgattagattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGCTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCa                                                                                                                        >  1:356238‑M2/37‑169 (MQ=255)                                                                                                                                                                                                            ttgaaatctttcaatctgattagattaggtggccgtAACCGGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCa                                                                                                                        <  2:356238‑M2/133‑1 (MQ=255)                                                                                                                                                                                                                       caatcgaattagattagtttgccgtACCCTGCAGGTTTTCTATCAGGCCGGGTTTATTAAATCCTTCAATCCCAGACTTTCCGCCAGCCGATGCAGGTGGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGGTCCAGTGCCTGACGAATAGTTGCACGCTGGAACGCTTCTGTCGCTCCACGCAGGTTTTGCTGAACAACGGGCACCGCCGCCACTTCTGGCGTCGGCa                           <  1:511834‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                          tctgattagattcggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCAAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTAACAACGGGCACCGCCGCCACTTCTGGCGTCGGCAACg                        >  2:62676‑M2/23‑251 (MQ=255)                                                                                                                                                                                                                           gtgattagattgggttgccttACCCTGCAGGTTTTCTATCAGGCCGGGATTATTCAATCTTTCAACCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGCCGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGGGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTAACAACGGGCACCGCCGCCACTTCTGGCGTCGGCAACGt                       <  1:415889‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                 agattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCAc                                                                                                   <  1:539041‑M2/154‑1 (MQ=255)                                                                                                                                                                                                                                 agattaggttgccgtAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCAc                                                                                                   >  2:539041‑M2/16‑169 (MQ=255)                                                                                                                                                                                                                                                aaCCTGCAGGATTTCTCCCAGGCCAGGATTATTTAATCCTACAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATACGCGCGCAGGCAGCCCAGTGGTGATGATTTTGTGCCAGTGCCTGACGAATAGATTCACGACGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTTACAACGGGCACCGCCGACACTTCTGGCGCCGGCAACGACACCTCAGGACAAGCACAATg  >  2:201911/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           GGCTGGTGGCGGTATGTGACGCTCAACAACAGGTACAAGGCGTCTTTACCGACGGCGATTTACGTCGCTGGCTGGTTGGCGGCGGCGCACTCACCACGCCAGTCAATGAAGCGATGACGGTCGGCGGCACCACGTTGCAATCGCAAAGTCGCGCCATCGACGCCAAAGAGATCCTGATGAAGCGCAAAATCACTGCCGCACCGGTGGTGGATGAAAACGGCAAACTCACCGGCGCAATAAACCTGCAGGATTTCTATCAGGCCGGGATTATTTAATCCTTCAATCCCAGACGTTTCGCCAGCCGATGCAGGTTGGCGACGTCGGTTTCCAGCATCCGCGCGCAGGCAGCCCAGTTGTGATGATTTTGTGCCAGTGCCTGACGAATAGTTTCACGCTGGAACGCTTCTGTCGCTTCACGCAGGTTTTGCTTAACAACGGGCACCGCCGCCACTTCTGGCGTCGGCAACGTCACCTCAGGAAAAGCAAAATG  >  NC_000913/2830504‑2830993

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