Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,737,689	IS150 (+) +3 bp	intergenic (+43/‑54)	pheL → / → pheA	phe operon leader peptide/fused chorismate mutase/prephenate dehydratase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2737689 =	0 (0.000)	14 (0.530)	14/488	0.7	100%	intergenic (+43/‑56)	pheL/pheA	phe operon leader peptide/fused chorismate mutase/prephenate dehydratase
	?	NC_000913	= 3722075	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region
*	?	NC_000913	= 2737691	0 (0.000)	23 (0.870)	23/486	0.3	100%	intergenic (+45/‑54)	pheL/pheA	phe operon leader peptide/fused chorismate mutase/prephenate dehydratase
	?	NC_000913	3720633 =	NA (NA)					noncoding (1/1443 nt)	IS150	repeat region

AACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTGATCGTGAACGCGATTTGCTGGAAAGATTAATTACGCTC  >  NC_000913/2737447‑2737930                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                agaaaaggcaattgtctggataatgctgtggtggagtgtttctttggaaccttaaagtcggagtgtttttatcttgatgagttcagtaatataagcgaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGc                                                                                                                                                                                                                                           <  2:54839‑M1/9‑1 (MQ=255)       ggTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcaccc                                                                                                                                                                                                                                     >  2:69700‑M1/1‑239 (MQ=255)                cGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttg                                                                                                                                                                                                                            <  1:69700‑M1/251‑22 (MQ=255)                       atgctgtggtggagtgtttctttggaaccttaaagtcggagtgtttttatcttgatgagttcagtaatataagcgaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCtttt                                                                                                                                                                                                                     >  1:262212‑M1/221‑251 (MQ=255)                           tcgAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTACAGTCACATAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaat                                                                                                                                                                                                                 >  2:183334‑M1/1‑219 (MQ=255)                                 ggagtgtttctttggaaccttaaagtcggagtttttttatcttgatgagttcagtaatataagcgaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtggttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATaa                                                                                                                                                                                                           >  1:202836‑M1/211‑251 (MQ=255)                                                  ccttaaagtcggagtgtttttatcttgatgagttcagtaatataagcgaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTAt                                                                                                                                                                                          <  2:53931‑M1/58‑1 (MQ=255)                                                             gAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacc                                                                                                                                                                               >  1:106388‑M1/1‑185 (MQ=255)                                                                                    cagtactataagcgcactgaaggatgctgttacggcatctattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgag                                                                                                                                                          <  2:113187‑M1/90‑1 (MQ=255)                                                                                      gcgcGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgt                                                                                                                                                      <  1:213865‑M1/251‑92 (MQ=255)                                                                                           ataagcgaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACAGCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCa                                                                                                                                                 >  2:13438‑M1/153‑251 (MQ=255)                                                                                            tttttGTTGACAGCGTGAAACCAGTACGGATACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTACACTGAATGGGAGGCGTTTCGTAGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtaaaagccaaaatacccttttgaaaagcgcgttgaagtagtgaatc                                                                                                                                                 >  2:218983‑M1/1‑154 (MQ=255)                                                                                                                       gggTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacg                                                                                                                                                                                                                        <  2:8231‑M1/152‑26 (MQ=255)                                                                                                                       gggTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacg                                                                                                                                                                                                                        >  1:8231‑M1/1‑127 (MQ=255)                                                                                                                           actgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatc                                                                                                                 >  1:6638‑M1/1‑123 (MQ=255)                                                                                                                                       tCACTTAAGGAAACAAACATGAAACCCATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAACAAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggaacatctcggaacgt                                                                                                     >  2:42483‑M1/1‑111 (MQ=255)                                                                                                                                         aCTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttt                                                                                                   >  1:37693‑M1/1‑109 (MQ=255)                                                                                                                                          cTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgacaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgtttt                                                                                                  >  2:124469‑M1/1‑108 (MQ=255)                                                                                                                                               ggAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccgtgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgt                                                                                             <  2:72761‑M1/251‑149 (MQ=255)                                                                                                                                                 aaacaaacATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcc                                                                                           >  1:243500‑M1/1‑101 (MQ=255)                                                                                                                                                                                caattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCa                                                            >  2:242469‑M1/68‑251 (MQ=255)                                                                                                                                                                                   aTTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctct                                                          >  1:126335‑M1/1‑67 (MQ=255)                                                                                                                                                                                    tgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGc                                                        <  1:29154‑M1/188‑1 (MQ=255)                                                                                                                                                                                         ttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaa                                                    <  2:237971‑M1/250‑190 (MQ=255)                                                                                                                                                                                             aCCTTCCCCTGCATGGGAGGCGTTTCGTCGTGTGAAACAGCATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgacatggaatagcccctaatatgtcaaagccaaaatacacttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaacaacat                                               >  2:122470‑M1/1‑57 (MQ=255)                                                                                                                                                                                             aCCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaatccaaaatacccttttgaaaagcgccgtgaaggcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacat                                               <  1:223464‑M1/251‑195 (MQ=255)                                                                                                                                                                                              ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatg                                              <  1:248447‑M1/251‑196 (MQ=255)                                                                                                                                                                                              aatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTg                                              >  1:113187‑M1/54‑251 (MQ=255)                                                                                                                                                                                               cTTCCACTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGCACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggggtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatgg                                             >  2:179553‑M1/1‑55 (MQ=255)                                                                                                                                                                                                                   gtttaactgtccaactttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTGATCGTGAACGCGAtt                         <  2:191813‑M1/219‑1 (MQ=255)                                                                                                                                                                                                                       taactgtccaacttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGt                                                   >  1:154571‑M1/29‑221 (MQ=255)                                                                                                                                                                                                                       taactgtccaacttttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGt                                                   <  2:154571‑M1/193‑1 (MQ=255)                                                                                                                                                                                                                        tcgtGTGAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatggagaaaaaggtttaattcccaaacct                    >  2:71891‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                              gAAACAGAATGCGAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacagtatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgctctctatgaaaaacatggagaaaaaggtttaattcccaaacctaaaggc              <  1:218983‑M1/251‑228 (MQ=255)                                                                                                                                                                                                                                     ttggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAAGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTGATCGTGAACGCGATTTGCTGGCAAGattaatt        <  2:58251‑M1/236‑1 (MQ=255)                                                                                                                                                                                                                                       ggggtcagtacaCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTGATCGTGAACGCGATTTGCTGGAAAGATTAATTACg     <  2:172185‑M1/239‑1 (MQ=255)                                                                                                                                                                                                                                          gAAGACGAACAAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatggagaaaaaggtttaattcccaaacctaaaggcgttagtgctgat  >  2:54334‑M1/1‑12 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                AACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTGATCGTGAACGCGATTTGCTGGAAAGATTAATTACGCTC  >  NC_000913/2737447‑2737930

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