Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,633,913	IS5 (–) +4 bp	coding (292‑295/591 nt)	pinQ →	putative recombinase PinQ

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1633913 =	0 (0.000)	17 (0.890)	14/276	0.2	100%	coding (292/591 nt)	pinQ	putative recombinase PinQ
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1633916	NA (NA)	20 (1.030)	16/278	0.1	NA	coding (295/591 nt)	pinQ	putative recombinase PinQ
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

ATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAAC  >  NC_000913/2101525‑2101961                                                                                                                                                                                                                                 |                                                                                                                                                                                                                      atgaGTATAATTTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCGCAGGTCAATATCTTCACAGTACATAAAGTGTCCTTGATCTAAGCCATGTACACGCACAAAATCTAAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCATCGACCGTATCATAATAGATACTTTCTTGAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCTGGGAAAttn                                                                                                                                                                                             <  1:99850/251‑2 (MQ=255)                 gCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTaagg                                                                                                                                                                                                                      >  1:425201/1‑210 (MQ=255)                 gCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTCTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAg                                                                                                                                                                            <  1:264052/251‑1 (MQ=255)                   aTGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACAGAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTGCTTTAGGAATTTTTGTTTGATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGtc                                                                                                                                                                          <  1:234770/251‑1 (MQ=255)                                                      ccagGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTAAATAAAGTAACCTTGATAGAAGCCATTTACACGCACAAAATATGAAAACAGTACCAGCATAAATGATCCTGCGCACAAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTg                                                                                                                                                      >  1:475591/1‑236 (MQ=38)                                                           cTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCggg                                                                                                                                                                                                    >  1:76445/1‑186 (MQ=255)                                                                           ggTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGATTAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTAATTTCATTTCTTATTGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                                                                                                  >  1:120502/1‑251 (MQ=255)                                                                                              cATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGtca                                                                                                                                                                         >  1:173285/1‑177 (MQ=255)                                                                                                        ncTTGATCGAAGCCATTTACACGCACAAAATCTGAAACACGTACCAGCATAAATGATCCTGCGCACCAATCAACACCCGTATCAGCATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGCATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcagctcgg                                                                                      >  1:63079/2‑243 (MQ=255)                                                                                                                     cATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGAGCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTGTTGGTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCAGTTCTTCATGTTGGAGCCGAGGTGTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAt                                                                                          <  1:229069/233‑1 (MQ=25)                                                                                                                              gCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACACCCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGACATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                                                                                                       >  1:239223/1‑211 (MQ=255)                                                                                                                               cacaAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtct                                                                                                                                        >  1:37813/1‑177 (MQ=255)                                                                                                                                         gAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCAGTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTccc                                                    <  1:322849/251‑1 (MQ=255)                                                                                                                                               cnganntgCANAAATGATCCTGCGCACCAATCAACNACCGTATCAGAATAGATACTTTCTTTAGGAANTNTTGGTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTC‑CGATCTNACTCAGGCATNTCANNTCTTNANGNTTGAGCCGATTTTTTCTCCCGTAAATGCCNTGAATCAGCCTATGTAGACCGTTNNTTNNNCANNTAAgg                                                                                    <  2:379108/206‑1 (MQ=18)                                                                                                                                               cGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCCTTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTa                                                                                                               >  1:37572/1‑186 (MQ=255)                                                                                                                                                gTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTGTTTAGTGAGATCTCTCCCACTGACg                                             <  1:389192/251‑1 (MQ=255)                                                                                                                                                      gCATAAATGCTCCTGCGCACCAATCAACAACCGTATCAGCATAGATACTTTCTTTAGGAATATTTGTTTGATTAAGGAAGGTGCGAATACGCGGAGAAATTCTTC‑CGATCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAACTGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg                                                                                    >  1:379108/1‑206 (MQ=18)                                                                                                                                                       naTAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGACTTTTTGCTTTATTAAGGAAGGTGCGAATAACCGGGGCAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCATATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGCCGTATCAtt                                      >  1:89171/2‑251 (MQ=255)                                                                                                                                                            aTGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTACCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTc                                 <  1:29493/251‑1 (MQ=255)                                                                                                                                                                 ccTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTGGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGtcatttca                                                                                                                                                                    <  1:233148/115‑1 (MQ=255)                                                                                                                                                                  cTGCGAACCAATCAACAACAGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                                                                                                   >  1:402945/1‑131 (MQ=255)                                                                                                                                                                               aacaacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGCGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGANNCAGtc                                                                                                                                                                          >  2:38487/1‑95 (MQ=255)                                                                                                                                                                                ncaacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGcc              >  1:129380/2‑250 (MQ=255)                                                                                                                                                                                      gTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGc               >  1:5432/1‑243 (MQ=255)                                                                                                                                                                                           aGAATAGATACTTTCTTTAGGAATTTTTGTTTTAGTAAGGAAGGTGCTAATAAGCGGGGAAAGTCTTCTCG‑GCTGACTCAGTCATGTCATTTCTTCATGTTTAAGCCTATTTTTTCTCCCGGAAATGCCTTGAACCAGCCTATTTTGTACTTTTGGGCCCCATGTGATGCTTTATCCCGAGTGTTTAGTTAGATCTCGTTCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAAc  <  1:503576/251‑1 (MQ=9)                                                                                                                                                                                              aTAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATACGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGCGCCg                                                                                                                                                 >  1:283697/1‑105 (MQ=255)                                                                                                                                                                                                             aNNAATTTTNNTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtct                                                                                                                                        <  2:37813/99‑1 (MQ=255)                                                                                                                                                                                                                       gTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa                                                                                       >  1:276423/1‑138 (MQ=35)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                      ATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG‑GCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAAC  >  NC_000913/2101525‑2101961

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