Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	1,301,457	(CTGGGATAAT)1→2	coding (276/1632 nt)	oppA →	oligopeptide ABC transporter periplasmic binding protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1301448 =	0 (0.000)	18 (0.400)	18/482	0.9	100%	coding (267/1632 nt)	oppA	oligopeptide ABC transporter periplasmic binding protein
	?	NC_000913	= 1301457	0 (0.000)					coding (276/1632 nt)	oppA	oligopeptide ABC transporter periplasmic binding protein

TGCGCTAATGGCAGGGAATGTCGCGCTGGCAGCTGATGTACCCGCAGGCGTCACACTGGCGGAAAAACAAACACTGGTACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATT  >  NC_000913/1301229‑1301690                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   gcagggaatgtcgcgctggcagctgatgtacccgcaggcgtcacactggcggaaaaacaaacactgggacgtaacaatggttcagaagttcagtcattagatccgcacaaaattgcaggtgtttcggagtataatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGgacat                                                                                                                                                                                                                     <  2:61097‑M1/32‑3 (MQ=255)           gCAGGGAATGTCGCGCTGGCAGCTGATGTACCCGCAGGCGTCACACTGGCGGAAAAACAAACACTGGGACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGCAGGTGTTTCGGAGTATAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggacat                                                                                                                                                                                                           <  2:61097‑M2/251‑33 (MQ=255)               gctggcagctgatgtacccgcaggcgtcacactggcggaaaaacaaacactggtacgtaacaatggttcagaagttcagtcattagatccgcacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTaa                                                                                                                                                                                                        <  1:167225‑M1/45‑1 (MQ=255)                        cgatgtacccgcaggcgtcacactggcggaaaaacaaacactggtacgtaacaatgggtcagaagttcagtcattagatccgcacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaaggcttacgggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGaaa                                                                                                                                                                                              >  2:14134‑M1/197‑251 (MQ=255)                         gcTGGCAGCTGATGTACCCGCAGGCGTCACACTGGCGGAAAAACAAACACTGGTACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaa                                                                                                                                                                                              <  1:167225‑M2/250‑46 (MQ=255)                                  cgaTGTACCCGCAGGCGTCACACTGGCGGAAAAACAAACACTGGTACGTAACAATGGGTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACGGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaa                                                                                                                                                                                    >  2:14134‑M2/2‑196 (MQ=255)                                   caggcgtcacactggcggaaaaacaaacactggtacgaaacaatggttaagaagttcagtcattagatccgcacaaaattgaaggtgttccggagactaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGCTAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATgg                                                                                                                                                                                   <  1:172088‑M1/66‑1 (MQ=255)                                             cAGGCGTCACACTGGCGGAAAAACAAACACTGGTACGAAACAATGGTTAAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGACTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggctaataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatgg                                                                                                                                                                         <  1:172088‑M2/251‑67 (MQ=255)                                             actggcggaaaaacaaacaatggtacgtaaaaatggttcagaaggtcagtcattagatccgcacaaaattgaaggggttccggagtctaatatcagacgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctggggtaatCTGGGATAATAAAGACGCGACAGTCTGGCCCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTc                                                                                                                                                                         <  2:34798‑M1/76‑1 (MQ=255)                                                 gcggaaaaacaaacactggtacgtaacaatggttcagaagttcagtcattagatccgcacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaagtcttactggtcagcgatcttgacgttcatccagcacctggcgtcgctgactcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAg                                                                                                                                                                     >  1:243798‑M1/172‑251 (MQ=255)                                                      aaaacaaacactggtacgtaacaatggttcagaagttcagtcattagatccgcacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTTCATTTGCGTTAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACaa                                                                                                                                                                >  1:164098‑M1/167‑251 (MQ=255)                                                       acTGGCGGAAAAACAAACAATGGTACGTAAAAATGGTTCAGAAGGTCAGTCATTAGATCCGCACAAAATTGAAGGGGTTCCGGAGTCTAATATCAGACGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGGTAATctgggataataaagacgcgacagtctggcccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtc                                                                                                                                                               <  2:34798‑M2/251‑77 (MQ=255)                                                           gCGGAAAAACAAACACTGGTACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGTCTTACTGGTCAGCGATCTTGACGTTCATCCAGCACCTGGCGTCGCTGACTCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacag                                                                                                                                                           >  1:243798‑M2/1‑171 (MQ=255)                                                           aaacactggtacgtaacaatggttcagaagttcagtcattagatccgcacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACtt                                                                                                                                                           <  2:102177‑M1/90‑1 (MQ=255)                                                                aaaaCAAACACTGGTACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggacctttcatttgcgttaagatgcgaaatggtctgatggcacgccagtcacagcacaa                                                                                                                                                      >  1:164098‑M2/1‑166 (MQ=255)                                                                    tacgtaacaatggttcagaagttcagtcattagatccgcacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGCCCTTCCATTTGCGTAACGATGAGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGtat                                                                                                                                                    <  2:447240‑M1/97‑1 (MQ=255)                                                                     aaacaCTGGTACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactt                                                                                                                                                 <  2:102177‑M2/251‑91 (MQ=255)                                                                             atggttcagaagttcagtcatttgatccgcccaaaattgaaggtgatccggagtataatatcagccgagacctgtttgaatgcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTGACATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTGTGTGTATAGCTGGCAACg                                                                                                                                         <  1:350413‑M1/108‑1 (MQ=255)                                                                              tACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggcccttccatttgcgtaacgatgagaaatggtctgatggcacgccagtcacagcacaagactttgtgtat                                                                                                                                          <  2:447240‑M2/249‑98 (MQ=255)                                                                                       aTGGTTCAGAAGTTCAGTCATTTGATCCGCCCAAAATTGAAGGTGATCCGGAGTATAATATCAGCCGAGACCTGTTTGAATGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggacctgacatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactgtgtgtatagctggcaacg                                                                                                                               <  1:350413‑M2/251‑109 (MQ=255)                                                                                                       ccgtacaaaattgaaggtgttccggagtctaatatcagccgagccctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctgtcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTtct                                                                                                                >  2:12477‑M1/118‑250 (MQ=255)                                                                                                          cacaaaattgaaggtgttccggagtctaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGt                                                                                                            <  1:393701‑M1/137‑1 (MQ=255)                                                                                                                 ccGTACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGCCCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGTCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaacactgcttct                                                                                                      >  2:12477‑M2/1‑117 (MQ=255)                                                                                                                    cacaAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaacactgcttctccgt                                                                                                  <  1:393701‑M2/251‑138 (MQ=255)                                                                                                                               gagtctaatatcagccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATg                                                                                       <  1:209705‑M1/158‑1 (MQ=255)                                                                                                                                         gAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaacactgcttctccgtatgccagttatctgcaatatg                                                                             <  1:209705‑M2/251‑159 (MQ=255)                                                                                                                                             ccgagacctgtttgaaggcttactggtcagcgatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGt                                                                         >  2:209705‑M1/80‑251 (MQ=255)                                                                                                                                                       ccGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaacactgcttctccgtatgccagttatctgcaatatgggcatatcgccggt                                                               >  2:209705‑M2/1‑79 (MQ=255)                                                                                                                                                                            gatcttgacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAACTGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCg                                           <  1:422544‑M1/202‑1 (MQ=255)                                                                                                                                                                                  gacggtcatccagcacctggcgtcgctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTAcc                                     >  1:129448‑M1/43‑250 (MQ=255)                                                                                                                                                                                     gggtatccatcccccggcgtggctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGGCTGGCCCGTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACTCTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATGGATGAAATTCTTGAAGGGAAAAAACCGATTACCGAtc                                 >  2:259882‑M1/40‑251 (MQ=255)                                                                                                                                                                                      gATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaactggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaacactgcttctccgtatgccagttatctgcaatatgggcatatcgccggtattgatgaaattcttgaagggaaaaaaccg                                 <  1:422544‑M2/250‑203 (MQ=255)                                                                                                                                                                                            gACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaacactgcttctccgtatgccagttatctgcaatatgggcatatcgccggtattgatgaaattcttgaagggaaaaaaccgattacc                           >  1:129448‑M2/1‑42 (MQ=255)                                                                                                                                                                                               gggtatccatcccccggcgtggctgaatcctgGGATAATctgggataataaagacgcgaaaggctggcccgtccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttctgttgatccgaactctgcttctccgtatgccagttatctgcaatatgggcatatcgccggtatggatgaaattcttgaagggaaaaaaccgattaccgatc                       >  2:259882‑M2/16‑39 (MQ=255)                                                                                                                                                                                                          gctgaatcctgggataatCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCAGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATtgatg            >  2:162064‑M1/19‑251 (MQ=255)                                                                                                                                                                                                                    gCTGAATCCTGGGATAATctgggataataaagacgcgaaagtctggaccttccatttgcgtaaagatgcgaaatggtctgatggcacgccagtcacagcacaagactttgtgtatagctggcaacgttcagttgatccgaacactgcttctccgtatgccagttatctgcaatatgggcatatcgccggtattgatgaaattcttgaagggaaaaaaccgattaccgatctcggcgtgaaagctattgatg  >  2:162064‑M2/1‑18 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   TGCGCTAATGGCAGGGAATGTCGCGCTGGCAGCTGATGTACCCGCAGGCGTCACACTGGCGGAAAAACAAACACTGGTACGTAACAATGGTTCAGAAGTTCAGTCATTAGATCCGCACAAAATTGAAGGTGTTCCGGAGTCTAATATCAGCCGAGACCTGTTTGAAGGCTTACTGGTCAGCGATCTTGACGGTCATCCAGCACCTGGCGTCGCTGAATCCTGGGATAATAAAGACGCGAAAGTCTGGACCTTCCATTTGCGTAAAGATGCGAAATGGTCTGATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATT  >  NC_000913/1301229‑1301690

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