Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,016,437	Δ2 bp :: IS186 (–) +8 bp :: Δ1 bp	coding (570‑577/687 nt)	fliH →	flagellar biosynthesis protein FliH

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2016437 =	1 (0.020)	62 (0.950)	58/482	0.2	98.4%	coding (570/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	2514274 =	NA (NA)					noncoding (2/1345 nt)	IS186	repeat region
*	?	NC_000913	= 2016444	4 (0.060)	72 (1.100)	69/484	0.1	94.6%	coding (577/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	= 2515615	NA (NA)					noncoding (1343/1345 nt)	IS186	repeat region

GACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCA  >  NC_000913/2515378‑2515862                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        gacgacATAATCCAGCCAACGCTGGATTTCCCCCCCAGAAGAGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGATGCGTCGCATTCTTCAGATGGTTCAACCCGTAAGTTAGCGCGTATGGGATAATTCCCCgg                                                                                                                                                                                                                                            <  2:84731/251‑1 (MQ=255)          aTCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTAc                                                                                                                                                                                                                                   <  1:640958/251‑1 (MQ=38)          aTCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGCATCCGAGAAGAAGAACTAACTCGCTGTGGAGAATAACAAAAATGGTCATCTGGAACTTCCAGGTGGCCATTCGTCGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGCAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTAc                                                                                                                                                                                                                                   >  1:309078/1‑251 (MQ=11)            ccagccagCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGc                                                                                                                                                                                                                                 <  2:174362/251‑1 (MQ=38)                ccATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGt                                                                                                                                                                                                                             >  1:157825/1‑251 (MQ=40)                ccATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATGCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGt                                                                                                                                                                                                                             <  1:556725/251‑1 (MQ=255)                   ccGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                                                                                                                                                                                                                          <  2:57949/250‑1 (MQ=255)                    cGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                                                                                                                                                                                                         <  1:377619/251‑1 (MQ=255)                     gCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                                                                                                                                                                                                                        >  2:159736/1‑251 (MQ=255)                                   ccAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCg                                                                                                                                                                                                          <  1:566186/251‑1 (MQ=255)                                      gAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc                                                                                                                                                                                                       <  1:545168/251‑1 (MQ=255)                                                     aaaaGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc                                                                                                                                                                                        >  1:51900/1‑251 (MQ=255)                                                       aagaagAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTAGCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCATCGAGGCGTCGCATTCGTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTGAATCAGCTAATGGCTACGGGGGACCAATCAAAGACCCt                                                                                                                                                                                      <  2:612722/251‑1 (MQ=17)                                                       aagaagAACTAACTAGTTGTGGAGACTAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCt                                                                                                                                                                                      >  1:201691/1‑251 (MQ=255)                                                        agaagaACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                                                                                                                                                                                                         >  1:137649/1‑215 (MQ=255)                                                        agaagaACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                                                                                                                                                                                                         <  2:137649/215‑1 (MQ=255)                                                         gaagaaCTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTg                                                                                                                                                                                    >  1:135955/1‑251 (MQ=255)                                                                     cGTTGTGGAGAATAACGAACATGGTCATCTGGAGCTTACAGGCGCCCATTCGTGGGACAGTATCCCTGACAGCCT‑CAACAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTCCGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTACTCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCAc                                                                                                                                                                        >  2:444895/1‑251 (MQ=18)                                                                     cGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACCGGTGGCCATTCGTGGGCCAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCACATACCCTTGCAGATATTGCAc                                                                                                                                                                        >  2:581778/1‑251 (MQ=255)                                                                           ggAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCg                                                                                                                                                                    >  1:678421/1‑249 (MQ=255)                                                                                aTAACAAAAATGGTCATCGGGAGCTTACAGGTGGCCATTCGTGGGACAGTATACCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACTAGGTACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCTGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTAt                                                                                                                                                             <  1:661635/251‑1 (MQ=255)                                                                                                 cTGGAGCTTACAGGTGGCCATTCGTGGGGCAGTATCCCTGACAGCCTACAA‑AAAGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACGAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCcttctgc                                                                                                                                            >  2:406186/1‑251 (MQ=255)                                                                                                         tACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGt                                                                                                                                    >  2:432281/1‑251 (MQ=255)                                                                                                                       cGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGGTGCGTCTCGACAAACTcc                                                                                                                       <  1:179965/250‑1 (MQ=255)                                                                                                                               aGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGAGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAACTACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACcacc                                                                                                              >  1:58864/1‑250 (MQ=255)                                                                                                                                gTATCCCTGACAGCCTACAA‑ACCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGCGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCacac                                                                                                             >  1:573320/1‑251 (MQ=255)                                                                                                                                          cAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCAGCGTCTTAACGAGGCACCGAATCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACTCCGGTTAATCATCTTATGGCTGCGGGCGACCGATCAAATACCCTTGCAGATATTGCACCCCGATCTTATGCAATAGCGCCTGTTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACtt                                                                                                   <  1:435797/251‑1 (MQ=37)                                                                                                                                           aGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTc                                                                                                  <  1:87599/251‑1 (MQ=255)                                                                                                                                                aCAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCttt                                                                                             >  1:578227/1‑251 (MQ=255)                                                                                                                                                aCAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCt                                                                                               <  1:179415/249‑1 (MQ=255)                                                                                                                                                 cAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTc                                                                                            >  1:195250/1‑251 (MQ=255)                                                                                                                                                     aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGcc                                                                                         >  1:634564/1‑251 (MQ=255)                                                                                                                                                           aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCg                                                                                   <  1:307806/251‑1 (MQ=255)                                                                                                                                                             tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGt                                                                                 >  1:235335/1‑251 (MQ=255)                                                                                                                                                               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgct                                                                                                                                              <  2:564604/188‑1 (MQ=255)                                                                                                                                                               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgct                                                                                                                                              >  1:564604/1‑188 (MQ=255)                                                                                                                                                                      cgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGAtct                                                                        >  1:291032/1‑251 (MQ=255)                                                                                                                                                                       gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGAtctc                                                                       >  1:546091/1‑251 (MQ=255)                                                                                                                                                                            aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTTCACCCCGAGCTTATGCAATAGCGCCTGCTTCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTc                                                                                                  <  1:604047/219‑1 (MQ=255)                                                                                                                                                                            aTCGTCTTAACGAGGAAACGAGGCGTAGCATTCTTCAGATGGCTCAACCCTCACGTTCGCGCTTATGGGATAATCCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCATTGCAGATATTGCACCCCGAGCTTCTGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTc                                                                                                  >  2:604047/1‑219 (MQ=255)                                                                                                                                                                                   tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGc                                                            >  1:176958/1‑250 (MQ=255)                                                                                                                                                                                   tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGCTTGCGCTTATGGGATCATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCgg                                                                                                                                                                                                         >  1:649950/1‑109 (MQ=21)                                                                                                                                                                                               gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTc                                               >  1:95034/1‑251 (MQ=255)                                                                                                                                                                                                  gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTg                                            <  2:135955/251‑1 (MQ=255)                                                                                                                                                                                                       gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTCATCAGCTCATGGCTGCGGGCGAACAATCAAATACCCTTGCAGCTATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACACTCATCGCATCCAGCGTGTTCGTGACAat                                       >  1:657732/1‑251 (MQ=255)                                                                                                                                                                                                              tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTc                                >  1:271466/1‑251 (MQ=255)                                                                                                                                                                                                                 gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTcac                              >  1:328973/1‑250 (MQ=255)                                                                                                                                                                                                                        cAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTGTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGc                      <  2:95034/251‑1 (MQ=255)                                                                                                                                                                                                                         aaCCCTTAAGTTAGCGCTTAGGGGATAATTCACCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCATCCCGAGCGTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACGCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGcc                     <  2:258923/251‑1 (MQ=255)                                                                                                                                                                                                                               tAAGTGAGCGCTCATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCGAATTGCTTCGGGCGACCAATCAAATACCCTTGCAGAGATTGCACCCCGAGCTTATGCAATAGCGCCTGCGGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAAGGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACTCCATCaa                <  2:93283/250‑1 (MQ=255)                                                                                                                                                                                                                                 aGTTAGCGCTTATGGGATAATTCCCAGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGCCGCCATCAAttt             >  2:612500/1‑251 (MQ=255)                                                                                                                                                                                                                                    tAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCAGCGCATCCAGCGTGTTCGTGACAATATCTTTCGCAAAGACGCCATCAATTTtg           <  2:51900/250‑1 (MQ=255)                                                                                                                                                                                                                                            aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCa  <  1:576540/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        GACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAA‑AACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCA  >  NC_000913/2515378‑2515862

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