Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,406,612	IS2 (–) +5 bp	coding (26‑30/939 nt)	lrhA ←	DNA‑binding transcriptional dual regulator LrhA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2068941 =	NA (NA)	34 (0.760)	30/480	0.4	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	= 2406616	0 (0.000)					coding (26/939 nt)	lrhA	DNA‑binding transcriptional dual regulator LrhA
*	?	NC_000913	2406612 =	0 (0.000)	31 (0.690)	31/482	0.4	100%	coding (30/939 nt)	lrhA	DNA‑binding transcriptional dual regulator LrhA
	?	NC_000913	= 2997691	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

ATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCC  >  NC_000913/2068698‑2069177                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             aTAGTTTCTACTGCAACCGCCACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCATTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATCATTTGCAGGACATCCCCGGAGCAGAGATAACAGACCACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGACTGGATTTg                                                                                                                                                                                                                                       >  2:122131/1‑251 (MQ=255)      ttCTAATGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCt                                                                                                                                                                                                                                  >  2:333267/1‑251 (MQ=255)                  cGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGAc                                                                                                                                                                                                                      >  1:170519/1‑251 (MQ=255)                    cAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACAt                                                                                                                                                                                                                    >  1:233644/1‑251 (MQ=255)                         tcttcAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTg                                                                                                                                                                                                                 <  1:108768/249‑1 (MQ=255)                                      aaaTCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCACCTATATTTCCAGACATCTGTTATCACTTAAcca                                                                                                                                                                                                   >  1:292993/1‑249 (MQ=255)                                           aGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATGGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGATATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGTATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACa                                                                                                                                                                                             <  1:297792/251‑1 (MQ=255)                                                       aCGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                                                                                                                                 >  2:157017/1‑251 (MQ=255)                                                            cgtcgtGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAg                                                                                                                                                                              <  1:157017/249‑1 (MQ=255)                                                                        gAATCGTACGGCCAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                                                                                                                                >  1:333629/1‑251 (MQ=255)                                                                                cGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGCCAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATAGTTCCAGACATCTGTTATCACTTAACCTATTACAAGCCTGCTGCCGCAGAGATTCCCGTGGCGAGCGATa                                                                                                                                                         <  1:122131/250‑1 (MQ=255)                                                                                                    aGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                                                                                                                                    <  2:248110/251‑1 (MQ=255)                                                                                                                  aCAGGTTATTGACCAAAAAACATTGCAGAACTCCACCAAAGCCTCCATAGCCGATAAATTGCAGGACATCCCCGGAGGAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGACGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAAt                                                                                                                      <  2:364229/251‑1 (MQ=255)                                                                                                                         aTTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGa                                                                                                                >  1:445441/1‑250 (MQ=255)                                                                                                                              ccAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCct                                                                                                          >  1:157721/1‑251 (MQ=255)                                                                                                                                  caaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGc                                                                                                      <  1:203569/251‑1 (MQ=255)                                                                                                                                      caTTGCAGAACTCCCCCAATGCCTCCATAGCAGATAATTTGCAGGACATCCCCGGCGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCAAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAgg                                                                                                  >  2:62299/1‑251 (MQ=255)                                                                                                                                                                     cGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTgg                                                                   <  2:233644/251‑1 (MQ=255)                                                                                                                                                                     cGATAATATGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATCTTCCCGTGGCGAGCGATAACCCCGCGCACTATGCGGATGCCATTCGTTATAATGCTCGATCGCCTCTGCAAGGTTCTTTGCTGCCGTTACCCCGTCTGGTTTgg                                                                   >  2:272841/1‑251 (MQ=255)                                                                                                                                                                                        cccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTc                                                >  1:79421/1‑251 (MQ=255)                                                                                                                                                                                                         aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTc                                >  1:219791/1‑250 (MQ=255)                                                                                                                                                                                                                           ccGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTc             >  2:387844/1‑251 (MQ=40)                                                                                                                                                                                                                                      tCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTcc  >  2:92721/1‑251 (MQ=35)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             ATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCC  >  NC_000913/2068698‑2069177

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