Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	3,888,720	Δ6 bp	intergenic (+211/‑5)	tnaC → / → tnaA	tnaAB operon leader peptide/tryptophanase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3888719	3 (0.090)	20 (0.590)	18/484	0.7	87.0%	intergenic (+210/‑11)	tnaC/tnaA	tnaAB operon leader peptide/tryptophanase
	?	NC_000913	3888726 =	3 (0.090)					intergenic (+217/‑4)	tnaC/tnaA	tnaAB operon leader peptide/tryptophanase

ATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTA  >  NC_000913/3888545‑3888968                                                                                                                                                                               |                                                                                                                                                                                                                                                          aTTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCCTCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGCTGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgc                                                                                                                                                                               >  1:151611‑M1/1‑175 (MQ=255)      aTTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgctt                                                                                                                                                                           >  1:43751‑M1/1‑170 (MQ=255)          actgtggtCTATTTGTTTGCTATGTCTGAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGAAATATTCATAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatatcactgaaccattcagcattcgtgttattgagccagtaacacgtaccactcgcgattatcgt                                                                                                                                                                      >  2:271172‑M1/3‑166 (MQ=255)                                         tttGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATCTCTCGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtactggaaaactttaaacatcaccctgaaccgttccgcattcgtgttattgagcaagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtctgaacc                                                                                                                                       >  1:39257‑M1/1‑135 (MQ=255)                                                 ttGCATAGGCCATCTCTCGTTTATTTACTTGTTGTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTAGAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctcgctgaacctttccgcatccgtgtgatttagccagtgaaacggacctctcgcgcttatcgtgaagacgcattgattgaatccggtgtgaacccgttcct                                                                                                                                <  2:340718‑M1/250‑124 (MQ=255)                                                       aGGTCATCTCTCGTTTATTTACTTGTTTTAATAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggat                                                                                                                         >  2:103467‑M1/1‑121 (MQ=255)                                                           cATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGCAATTAAAATAAATGAAGGATTATGTAATGGAAAACCTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAAGCCGGTATGAACCCGTTCctgctg                                                                                                                      <  2:125694/250‑1 (MQ=255)                                                                   gTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAAGAAATGAAtttaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttattgtgaagaggcaattattaaatccggtatgaacccgtgcctgctggatagcgaagatgct                                                                                                             <  2:230477‑M1/251‑143 (MQ=255)                                                                      tatttaCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgttttt                                                                                                          >  2:182478‑M1/1‑106 (MQ=255)                                                                        tttaCTTGTTTTAGTAAATGATGGTGCTTGAATATATCTCTGGCGAAGTCATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaacctttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaacacgtacaactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgcagatgtttttat                                                                                                        >  1:325093‑M1/1‑104 (MQ=255)                                                                             ttGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCa                                                                                                                                                                              >  2:97587/1‑176 (MQ=255)                                                                             ttGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCa                                                                                                                                                                              <  1:97587/176‑1 (MQ=255)                                                                                     ataAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttgctgacc                                                                                           <  1:103467‑M1/250‑161 (MQ=255)                                                                                            atgGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGCTGTAATATTCACAGTGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactaaccgacagcg                                                                                    >  1:122815‑M1/1‑84 (MQ=255)                                                                                                 gCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcacc                                                                               >  2:52320‑M1/1‑79 (MQ=255)                                                                                                               cTGGCGACTTAATCGGGATAGCAGATGTAATATTCACAGGGATCACTGTACTTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattcttaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgttttactgaccgacagcggcaccggggcggtgacgcc                                                                 >  2:141647‑M1/1‑65 (MQ=255)                                                                                                                 ggCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttagcgatttactgaccgacagcggcaccggggcggtgacgcaga                                                               <  2:43751‑M1/251‑189 (MQ=255)                                                                                                                             ggTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttactgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcaccggggcggtgacgcagagcatgcaggctg                                                   <  1:246520‑M1/251‑201 (MQ=255)                                                                                                                                   gCAGAGGTAATCTTCACAGGGATCCCTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcacttattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcaccggggcggtgacgcagagcatgcaggctgcgatg                                              >  2:19779‑M1/1‑45 (MQ=255)                                                                                                                                                   caGGGATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcaccggggcggtgacgcagagcatgcaggctgcgatgatgcgcggcgacgaagc                             >  1:123094‑M1/1‑29 (MQ=255)                                                                                                                                                     gggATCACTGTAATTAAAATAAATGAAtgtaatggaaaactttaaacatctccctgaacagatccgcattagtgttattgagccagtaaaacgtaccactcgcgcttatagtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcaccggggcggtgacgaaaagcatgcaggctgcgatgatgcgcggagacgaggcct                           >  2:307186‑M1/1‑27 (MQ=255)                                                                                                                                                      ggATCAATGTAATTAAAATAAATGAAtgtactggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccgggatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcaccggggcggtgacgcagagcatgcaggccgcgaggatgcgcggcgacgaagcct                           >  2:113208‑M1/1‑26 (MQ=255)                                                                                                                                                                              aatgtaatggaaaactttaaacatctccctgaaccgttccgcattcgtgttattgagccagtaaaacgtaccactcgcgcttatcgtgaagaggcaattattaaatccggtatgaacccgttcctgctggatagcgaagatgtttttatcgatttactgaccgacagcggcaccggggcggtgacgcagagcatgcaggctgcgatgatgcgcggcgacgaagcctacagcggcagtcgtagctactatgc  >  1:155113‑M1/1‑2 (MQ=255)                                                                                                                                                                               |                                                                                                                                                                                                                                                          ATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTA  >  NC_000913/3888545‑3888968

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.