Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,888,869	IS1 (+) +9 bp	coding (140‑148/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)	24 (0.660)	23/486	0.5	82.2%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 3888877	5 (0.140)					coding (148/1416 nt)	tnaA	tryptophanase
*	?	NC_000913	= 3584195	NA (NA)	40 (1.110)	39/486	0.1	88.5%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	3888869 =	5 (0.140)					coding (140/1416 nt)	tnaA	tryptophanase

CACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGA  >  NC_000913/3583948‑3584441                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       cACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAc                                                                                                                                                                                                                                                     <  1:187496/251‑1 (MQ=31)                 gATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                                                                                                                                                                                                                                     >  1:284571/1‑251 (MQ=38)                         ggACGGAGGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATGGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTTATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                                                                                                                                                                                                             <  2:32940/250‑1 (MQ=255)                                            gTATGAATCCCGCCTGAAGGGAAAGCGGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGGCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcaca                                                                                                                                                                                                           <  2:122434/250‑1 (MQ=255)                                                                                    aaGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCtt                                                                                                                                                                  >  2:82310/1‑251 (MQ=255)                                                                                             ccGCAGCGAATTGAACGGCATAACCTGAATCTAAGGCTGGCAC‑TGGCACGTCTGGGCCGGAAGTAGCTGTCGTCAT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTAGCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGGTTAGTGATTTCTCTTCATTCTTAGACGaaac                                                                                                                                                         <  2:71428/250‑2 (MQ=21)                                                                                                 aGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATg                                                                                                                                                     <  2:164488/251‑1 (MQ=255)                                                                                                    gAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAACCGAAATCATGGAc                                                                                                                                                  >  2:286335/1‑251 (MQ=255)                                                                                                          aGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAgttgt                                                                                                                                            <  1:39533/251‑1 (MQ=255)                                                                                                                        aaTCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGGGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATGACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGCTTTCTCTTCATTCTTAACCGAAATCATGGACAGTTGTTCATAATAATCaa                                                                                                                               >  2:137281/1‑250 (MQ=255)                                                                                                                        aaTCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCaa                                                                                                                               >  2:7184/1‑250 (MQ=255)                                                                                                                        aaTCTGAGGA‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGGCATCGGGCATTATCGGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTT‑GAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCaaaa                                                                                                                             <  2:270863/251‑1 (MQ=255)                                                                                                                            tGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATACGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCCGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAgg                                                                                                                           >  2:84115/1‑250 (MQ=255)                                                                                                                                     caCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACaaaa                                                                                                                   >  1:170715/1‑250 (MQ=255)                                                                                                                                           gCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCTGTTGAGCTGAATGACAAAGTCAACGGTCATGATCTGAACATAAAACACTATCAAAAAGTTGGAGGCATTACCTACATTAAAAGTAACTTTGTTGAAATTAGAGATAATATTGTCTCATAAAATATATTCAGTTCTTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAATAGGAAACAAAAATCATTg                                                                                                            <  2:15981/251‑1 (MQ=25)                                                                                                                                                gcggcaacgaaaGTCGCAGTCGTTTTACAAAA‑CGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCACTAAGTCGGAGTTATTGCCGACATTAAAATTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAAGATATTCAGTTCGTTTAGTGATTTCTCTGCATTCTTAAACGAAATCATGGACAGTTGTTCATAATATTCAAAAGGAAACAAAAATCAtcttaaac                                                                                                       <  2:262286/241‑8 (MQ=37)                                                                                                                                                       cGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACCCTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCa                                                                                                <  2:121777/251‑1 (MQ=255)                                                                                                                                                                     gTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCAGTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGAtt                                                                                  <  2:284571/251‑1 (MQ=255)                                                                                                                                                                       tct‑cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAAGCATTGTAAACCATGCCATCTTTTTAGAGATTa                                                                                 <  1:137281/250‑1 (MQ=255)                                                                                                                                                                         t‑cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGa                                                                                                                                                            >  2:330617/1‑173 (MQ=255)                                                                                                                                                                         t‑cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGa                                                                                                                                                            <  1:330617/173‑1 (MQ=255)                                                                                                                                                                                 cGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGaa                                                                        <  1:334917/250‑1 (MQ=255)                                                                                                                                                                                             aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCg                                                           <  1:7184/251‑1 (MQ=255)                                                                                                                                                                                                            ggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCa                                            <  1:84115/251‑1 (MQ=255)                                                                                                                                                                                                                                  cTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCAAAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACg                      >  2:333886/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTaca              <  1:155849/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGt    <  2:170715/251‑1 (MQ=255)                                                                                                                                                                                                                                                       tACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTga  <  1:82310/250‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       CACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGC‑AGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCT‑CAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGA  >  NC_000913/3583948‑3584441

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