Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,888,795	IS1 (+) +9 bp	coding (66‑74/1416 nt)	tnaA →	tryptophanase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	42 (1.610)	40/492	0.0	97.6%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 3888803	1 (0.040)					coding (74/1416 nt)	tnaA	tryptophanase
*	?	NC_000913	= 3584195	NA (NA)	31 (1.200)	31/490	0.1	96.8%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	3888795 =	1 (0.040)					coding (66/1416 nt)	tnaA	tryptophanase

ACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAG  >  NC_000913/3583949‑3584438                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aCCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAc                                                                                                                                                                                                                                                     >  2:113944/1‑250 (MQ=31)   ccTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAtt                                                                                                                                                                                                                                                  >  2:214497/1‑251 (MQ=32)            gggtgATATGGATGACAGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                                                                                                                                                                                                                                         <  2:53468/250‑1 (MQ=25)                 atatGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                                                                                                                                                                                                                    >  1:214824/1‑251 (MQ=38)                      gATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                                                                                                                                                                                                                               <  2:63184/251‑1 (MQ=40)                            ggATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                                                                                                                                                                                                                         <  2:142981/251‑1 (MQ=255)                                      ccGCTGAAGGAATCCCGCCTGAAGGGAAAGCGGGACGTTATACGCAAGCGATATACGCAGCGAATTGAGCGGCATACCCTGAATCTGAGGAAGCACCTGGCACGGCTGTGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGGGCTGCACGACAAAGTCAGCGGTCGTGGTCTGAACAGAACACACTATCAATAAGGGGGATTCATTACTGACATGAGAAGTAATGTTTATGAAATTAGAGATAAT‑ATTgg                                                                                                                                                                                                              <  2:118713/251‑1 (MQ=2)                                             aTGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGa                                                                                                                                                                                                        >  1:7767/1‑250 (MQ=255)                                                          gAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTACATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGtt                                                                                                                                                                                         <  2:187608/251‑1 (MQ=255)                                                                  aGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAAGCAGCACCTGGCACGGCTGGGGCGGCAGCCGCTGTCGTTATCAAAATCGGTGGAGCTGCCTGACACAGTCATCGGGCCTTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAgg                                                                                                                                                                                 >  2:51080/1‑250 (MQ=18)                                                                                                                      gAATCTGAGGCAGCACCTGGCACGGATGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATc                                                                                                                              <  2:1700/249‑1 (MQ=255)                                                                                                                            gAGGCAGCCCCTGGCACGGCTGGGACCGAAGTCGCTGTTGTTCTCTAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAg                                                                                                                         <  2:269777/248‑1 (MQ=255)                                                                                                                                   caCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAA‑TAAATTTTTTGAAATTAGAGAAAATTATGGGCA‑AGAAAAATTTATCCA‑TTTGTTTAATGATTTTTCCTCATTCTTTAAAGAAAACATGGAACGGTTGTCCA‑AATAATCAAAAGGGAACaaaaa                                                                                                               >  1:114387/1‑250 (MQ=9)                                                                                                                                           aCGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGTTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTAACTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTa                                                                                                       <  2:71323/251‑1 (MQ=255)                                                                                                                                             ggCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTaa                                                                                                      >  1:108858/1‑250 (MQ=255)                                                                                                                                                  ggACGGAAGTCGCTGGCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATGGTAACCCATg                                                                                                <  2:178424/251‑1 (MQ=255)                                                                                                                                                                gTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTagag                                                                                  <  2:142811/251‑1 (MQ=255)                                                                                                                                                                  cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                                                                                                                                                                                                             <  2:143907/97‑1 (MQ=38)                                                                                                                                                                  cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                                                                                                                                                                                                             >  1:143907/1‑97 (MQ=38)                                                                                                                                                                               ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGaaaa                                                                   <  1:204471/251‑1 (MQ=255)                                                                                                                                                                               ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGaaa                                                                    >  2:100010/1‑250 (MQ=255)                                                                                                                                                                                                                              ccTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGTTAAT‑AATGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAAGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAAGTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGTAAAAACTTCCAAc                    <  2:7367/250‑1 (MQ=255)                                                                                                                                                                                                                                tATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGACATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACgg                  >  1:156009/1‑251 (MQ=255)                                                                                                                                                                                                                                 aTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGa                 <  1:159277/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cATTACCGACAATACAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTTGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGATACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGTTATTGAAAAACTTTCAACGGAGCTACACTCCTTAAg  <  2:151077/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAAT‑ATTGGCACAGAAAA‑TATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGA‑CAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAG  >  NC_000913/3583949‑3584438

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