Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	389,898	Δ9,908 bp		yaiT–[yaiZ]	yaiT, insF‑2, insE‑2, yaiT, iprA, ampH, sbmA, yaiW, yaiY, [yaiZ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	389898	399805	9908	40 [0]	[0] 39	yaiT–[yaiZ]	yaiT,insF‑2,insE‑2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,[yaiZ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 389897	0 (0.000)	38 (1.390)	30/96	0.0	100%	intergenic (‑170/‑354)	hemB/yaiT	porphobilinogen synthase/putative autotransporter YaiT
	?	NC_000913	399806 =	0 (0.000)					intergenic (+1/+23)	yaiZ/ddlA	DUF2754 domain‑containing protein YaiZ/D‑alanine‑‑D‑alanine ligase A

TAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/389852‑389897 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT  >  NC_000913/399806‑399851                                                                                              TAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAA                                             <  1:1524107/49‑1   GAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATA                                           <  1:962505/49‑1   GAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATA                                           <  1:348718/49‑1    AAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATAT                                          >  1:1500379/1‑49     AACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATT                                         <  1:1875243/49‑1      ACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTA                                        >  1:1706002/1‑49         GCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCG                                     >  1:826535/1‑49         GCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCG                                     >  1:355884/1‑49          CAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGG                                    >  1:1160068/1‑49            ATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGC                                  <  1:85605/49‑1            ATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGC                                  >  1:1103797/1‑49            ATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGC                                  >  1:313454/1‑49             TAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCG                                 <  1:1957875/49‑1               ATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTT                               <  1:683844/49‑1               ATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGT                                <  1:428164/48‑1                TTTTATATTCACTGAAAATA‑TTTTAATCTTGAATATTATCGGGCGTTAA                             <  1:1504530/49‑1                 TTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTA                              >  1:851873/1‑48                   TATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAA                           >  1:323772/1‑49                        TCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTA                      >  1:575480/1‑49                         CACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTA                      >  1:61463/1‑48                          ACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACA                    <  1:1572521/49‑1                              AAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGT                <  1:1329154/49‑1                              AAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATT                  <  1:65502/47‑1                                AATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGG              >  1:283581/1‑49                                AATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGG              <  1:1897991/49‑1                                 ATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGT             >  1:247948/1‑49                                 ATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGT             <  1:449965/49‑1                                   ATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTT           >  1:1050755/1‑49                                    TTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTT           <  1:243898/48‑1                                    TTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTT           <  1:1597641/48‑1                                       TTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCA        <  1:766264/48‑1                                         AATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAT      <  1:2033749/48‑1                                         AATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAT      <  1:1368071/48‑1                                         AATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAT      <  1:1017284/48‑1                                          ATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGC    <  1:83555/49‑1                                            CTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT  >  1:1352377/1‑49                                            CTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT  >  1:1119480/1‑49                                            CTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTCTTCAATGCGT  >  1:2118796/1‑49                                                                                              TAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/389852‑389897 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT  >  NC_000913/399806‑399851

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