Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,049,903	IS2 (–) +5 bp	noncoding (5‑9/109 nt)	spf →	Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4049903 =	0 (0.000)	91 (1.870)	68/282	0.0	100%	noncoding (5/109 nt)	spf	Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 4049907	0 (0.000)	105 (2.150)	76/282	0.0	100%	noncoding (9/109 nt)	spf	Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGG  >  NC_000913/4498042‑4498320                                                                                                                                            |                                                                                                                                            ggggggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGa                                                                                                                                          <  2:441728/143‑1 (MQ=255)  gggggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGCAgcc                                                                                                                                         >  1:144535/1‑141 (MQ=255)    gggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATGTGATGGTCCCTTGTAGACTg                                                                                                                                       >  2:441318/1‑143 (MQ=255)        ttATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGCTTTGATGGTCCCTTGTAGACTGGccc                                                                                                                                   >  1:240829/1‑143 (MQ=255)           tCAGAACCCTCGGTTATTAATTATTTATATGCATTCAATAAAAAAATAACCTAATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCt                                                                                                                                <  2:483558/143‑1 (MQ=255)              taaTCCGTGGTTATTCATTATTTATATATATTCAATAAAAAAGTAACCTCATAAAAACTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCTCTGaa                                                                                                                             <  2:33445/142‑1 (MQ=255)              gAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATCTGATGGTCCCTTGTAGACTGGCCCCCTGaa                                                                                                                             >  2:326246/1‑143 (MQ=255)                  ccGCGGTTATTCATTATTTGTATACATTCAACAAAGAAGTAACCTCATGAGAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAAtctc                                                                                                                         <  2:200807/143‑1 (MQ=255)                   cgcgGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAACAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTcc                                                                                                                        <  1:192763/143‑1 (MQ=255)                    gcgGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCa                                                                                                                       >  2:200742/1‑143 (MQ=255)                      ggTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGa                                                                                                                     <  1:744737/143‑1 (MQ=255)                          atgctttagttatAGACATTCGATAAATAAGTATCCTCATCAACAAGCCTGTACGCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAAc                                                                                                                 <  2:246551/138‑1 (MQ=255)                          aTTCATTATTTATATCCATTCAATAAAAGAGTAACCTCATAAAAAAGCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAAc                                                                                                                 <  2:245278/143‑1 (MQ=255)                            tCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                                                               >  1:279340/1‑143 (MQ=255)                             cATTCTTTATACACACGCAATAAAAAAGTACCCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCaa                                                                                                              <  1:296989/143‑1 (MQ=255)                               ttatttatATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTTATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                                                            <  2:226964/143‑1 (MQ=255)                                tatttatATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                                                           <  2:602578/143‑1 (MQ=255)                                    tacataCCTTCCATAAAAAAGTAACCTCATAACAACTTCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTTGCCCCCTGCCTCTCCAGACAACCAATATCACt                                                                                                       <  2:444041/140‑1 (MQ=255)                                       ataCATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTaa                                                                                                    <  1:549889/143‑1 (MQ=255)                                       agaCAGTCAATACAAAAGTAACCTCATAAACGAGCCTGTACTGACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTaa                                                                                                    <  2:449565/141‑1 (MQ=38)                                           aTTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                                                                                <  1:117858/143‑1 (MQ=255)                                           aTTCAATAAAAAAGTAACCTCATAAAAAATCCGGTACTCACAGGGTGATGCGGGATAACTAACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                                                                                >  1:224953/1‑143 (MQ=255)                                               aaTAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGAGTGGCCCCCTGAATCTCCAGACAACCAATATCACGTAAATAAGTGa                                                                                            >  1:217153/1‑143 (MQ=255)                                               aaTAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                                                                            >  1:126232/1‑143 (MQ=255)                                                       aGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                                                                    >  1:256567/1‑143 (MQ=255)                                                        gTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                   <  1:525186/143‑1 (MQ=255)                                                        gTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                   <  1:290608/143‑1 (MQ=255)                                                        gTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCACTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                   <  2:36732/143‑1 (MQ=255)                                                                   aaaacTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                                                                        <  1:587957/143‑1 (MQ=255)                                                                   aaaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                                                                        <  2:514736/143‑1 (MQ=255)                                                                    taaaGCCTGTACTGACGGGGTGCTGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACGTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                       <  2:15858/142‑1 (MQ=25)                                                                    aaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                       <  1:620/143‑1 (MQ=255)                                                                    aaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                       <  2:173594/143‑1 (MQ=255)                                                                          cTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAg                                                                 <  1:244922/143‑1 (MQ=255)                                                                                tCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                                                           >  2:127247/1‑143 (MQ=255)                                                                                     gggTGATGCGGGATAACTGACGTAAGGACTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa                                                      <  1:673175/143‑1 (MQ=255)                                                                                      ggTGATGCGGTATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc                                                     <  2:614748/143‑1 (MQ=255)                                                                                           tGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCCCCAGACCACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTTTTTAGAGAACAgatg                                                >  2:136152/1‑143 (MQ=21)                                                                                             cGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAGAgatgat                                              >  2:114069/1‑143 (MQ=255)                                                                                              gggATAACTGACGTAAGGAGAGCTGCAGGGATTTGATGGTCCCTTGAAGACTGGCCCCCTAAAACACCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGACATAGGAGAACAgaaaaca                                             >  2:178455/1‑138 (MQ=18)                                                                                               ggATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattg                                            <  1:14214/143‑1 (MQ=255)                                                                                                gATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGATTTTAGACTAGTCAGTGGAGAACAGAtggttga                                           >  1:756297/1‑143 (MQ=255)                                                                                                     cTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCt                                      >  1:488470/1‑143 (MQ=255)                                                                                                      tGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                     >  2:493070/1‑143 (MQ=255)                                                                                                       gACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTa                                    <  1:130157/143‑1 (MQ=255)                                                                                                          gTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAggg                                 >  1:517400/1‑143 (MQ=255)                                                                                                            aaGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGcc                               <  1:38519/143‑1 (MQ=255)                                                                                                            aaGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGcc                               >  2:637148/1‑143 (MQ=255)                                                                                                             aGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTCGGGCCg                              >  1:297902/1‑143 (MQ=255)                                                                                                                          ggATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACGTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGGGAAGCGCAGAc                 >  2:731140/1‑143 (MQ=255)                                                                                                                          ggAGTTGATGTTCCCGTGTATACTTGCCCCCTGAATCTCCAGCCAACCCATATCACTTAAATAAGTGATATTCTTAATACTAGAGTTTAGACGAGTCACTGGAGAACACCTGATTGATTACATACGGCCGAAGAGAAGCAGAc                 >  2:848962/1‑143 (MQ=2)                                                                                                                                 aTGGTCGCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCTCTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAc          <  1:200742/143‑1 (MQ=25)                                                                                                                                 aTGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTCGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGtca          >  2:651594/1‑141 (MQ=255)                                                                                                                                       ccgtGTAGACTGGACCCTTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGGCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcacaca    <  2:240829/140‑1 (MQ=9)                                                                                                                                       ccTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGGTATAGTCTTAATACTCGTTTTTAGACAAGTCATTGGAGAACACATGAGTGATGTATTAGGGCCACGGAGACGCAGAAGGCATGCcacaca    >  2:613690/1‑143 (MQ=11)                                                                                                                                         ttGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAgg  >  1:677590/1‑143 (MQ=32)                                                                                                                                            |                                                                                                                                            GGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGG  >  NC_000913/4498042‑4498320

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