Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	REL606	2,322,340	IS186 (+) +6 bp	66.0%	coding (131‑136/963 nt)	menC ←	O‑succinylbenzoate synthase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2322340 =	9 (0.360)	7 (0.270)	5/272	NT	42.6%	coding (136/963 nt)	menC	O‑succinylbenzoate synthase
	?	REL606	= 2774196	NA (NA)					noncoding (1343/1343 nt)	IS186	repeat region
*	?	REL606	= 2322345	9 (0.360)	28 (1.070)	19/272	NT	74.8%	coding (131/963 nt)	menC	O‑succinylbenzoate synthase
	?	REL606	2772854 =	NA (NA)					noncoding (1/1343 nt)	IS186	repeat region

GAACACCCCTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGC  >  REL606/2772706‑2772993                                                                                                                                                     |                                                                                                                                            gAACACCCCTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACCCCTGGCCGAAATGTTGACTCGTCGAGCAAAACGGTTTATCCCCGCTGGCGCGGGGGACCCCCTGAACGGTGCGCAAAAccca                                                                                                                                            >  1:352024/1‑149 (MQ=255) gAACACCCCTCACACCGATTCGCCAAACGGCGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGATCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAacc                                                                                                                                            >  1:681466/1‑150 (MQ=255)      ccccTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTTATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                                                                                       >  1:195188/1‑150 (MQ=255)      ccccTCACACCGATTCGCCAAACGGTGGAGAAGCGGGTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGACCGGTGGGGAAAAACCACATaa                                                                                                                                       >  1:538040/1‑150 (MQ=255)       cccTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                                                                                      <  1:569387/150‑1 (MQ=255)       cccTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                                                                                      <  2:352917/150‑1 (MQ=255)         cTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcg                                                                                                                                    >  2:377894/1‑150 (MQ=255)                     cGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg                                                                                                                        <  1:392570/150‑1 (MQ=255)                      gCCAAGCGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTATGATGCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTCACCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                                                                                       <  1:632655/150‑1 (MQ=255)                      gCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                                                                                       <  2:539598/150‑1 (MQ=255)                      gCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                                                                                       <  2:229514/150‑1 (MQ=255)                      gCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                                                                                       <  2:683468/150‑1 (MQ=255)                      gCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                                                                                       <  2:312356/150‑1 (MQ=255)                            cGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                                                                                 >  2:585257/1‑150 (MQ=255)                            cGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                                                                                 >  2:300277/1‑150 (MQ=255)                              gTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGGACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTGATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGTAAAACCCCATCAGCGCTAACTTAAGGGTTGTGGGAt                                                                                                               >  1:421329/1‑150 (MQ=255)                                ggagaacgagttttaccccTCTGGCGCGGTTAACACCGAGGGGCAAGACTTCGTGTCATCGTGATGCCCGGTTTATCCAGGCTGGCGCGGGGAACACCCTGAGCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTGGTGGTATTa                                                                                                             <  1:307844/135‑1 (MQ=255)                                ggAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGTCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGTGTTGTGGTTTTa                                                                                                             >  1:652068/1‑150 (MQ=255)                                   gaAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                                          <  2:53079/150‑1 (MQ=255)                                                  cGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                                                                                           >  2:478534/1‑150 (MQ=255)                                                   ataagagacaggAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                          <  1:50501/140‑1 (MQ=255)                                                              gAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACgctgtctcttat                                                                               >  2:50654/1‑139 (MQ=255)                                                                acacAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGT‑GGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                            >  2:17667/1‑150 (MQ=255)                                                                             aagttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAc                                                                <  1:497521/147‑1 (MQ=255)                                                                             aaaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAc                                                                <  2:195751/150‑1 (MQ=255)                                                                             aaaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAc                                                                <  2:422531/150‑1 (MQ=255)                                                                                  ttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                           <  2:123395/150‑1 (MQ=255)                                                                                       tCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                      <  2:376355/150‑1 (MQ=255)                                                                                               atgcCCGGTTTAGCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAGCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGcc                                              <  1:583541/150‑1 (MQ=255)                                                                                               atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGcc                                              <  1:299428/150‑1 (MQ=255)                                                                                                   ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCata                                          >  2:230320/1‑150 (MQ=255)                                                                                                      gattATCCCCGCTGGCGCGGGGAACCCCCTGAGCGGTGGGGAAAAACCACATAAGCGCTAACTCAAGGGTTGTGGTATTACTCATGATATGGTCTAAGGTGGCGATGTATTACTCTCACGATAACTGGTCAGCAATTCTGGCGCagagtg                                       >  1:412455/3‑145 (MQ=18)                                                                                                      gTTTATGCCCGCTGGCGCGGGGAACACCCTGCACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATCAATTAACGTGCCGATTAAATACGCTCCCGAAAACTGGTCAGCCATTCTGGCCCATAttt                                       >  1:345515/1‑149 (MQ=18)                                                                                                      gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                       >  2:515136/1‑150 (MQ=255)                                                                                                      gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                       >  2:174292/1‑150 (MQ=255)                                                                                                         tATCCCCGCTGGCTCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                    <  1:477180/150‑1 (MQ=255)                                                                                                         tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                    <  2:543804/150‑1 (MQ=255)                                                                                                         gaTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAATCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                    <  1:376802/149‑1 (MQ=255)                                                                                                         gaTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                    <  1:420479/149‑1 (MQ=255)                                                                                                          aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTac                                   >  2:421678/1‑149 (MQ=255)                                                                                                            ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGc                                 <  2:653989/150‑1 (MQ=255)                                                                                                                  tgacgcgGGGAACACCCTGAACGGTGGCGAAAAACCCCCTAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTACGCCCgaag                           <  1:101170/147‑1 (MQ=255)                                                                                                                     cgcgGGGAGCACCCTGAACGGTGGGGAATGAGCGCATAAGCGATAACATAAGGGTTGTCGTAGTACGCCTGATATAATTTACCGTGCCGATGAATGACGCACACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAGGCCGGAAGAAc                        <  1:573220/150‑1 (MQ=14)                                                                                                                         gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa                    >  2:70976/1‑150 (MQ=255)                                                                                                                         gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGATTCTGGTATTACGCCTAATATGATTTAACGTGCCGATCAATTACTCTCACGATAACTGGTCCGCAACACTGGCCCATATTGGTAAGACCGAAGAACTGGa                    >  1:274208/1‑150 (MQ=17)                                                                                                                               acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc              >  2:288584/1‑150 (MQ=255)                                                                                                                               acCCTCAACGGTGGGGAAAAACCCCATAAGCGCTAACTTATGGGTTGTGGTCTTACGCCTGTTGCGATTTAACGTGCCGAGGCATCACTCCCACGATAACTGGTCCGCAAGTCCGGACCATAGTGGTAAGCCCGCTGAACTGGATACtta              >  1:195746/1‑149 (MQ=14)                                                                                                                                         gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt    >  2:55857/1‑150 (MQ=35)                                                                                                                                           ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACGGGATACTTCGGCCCGTAATGc  <  1:482447/150‑1 (MQ=21)                                                                                                                                                     |                                                                                                                                            GAACACCCCTCACACCGATTCGCCAAACGGTGGAGAAGCGGTTTATCCCCGCTGGCGCGGGGAACACAGCTGGGCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGC  >  REL606/2772706‑2772993

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