Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,262,613	IS1 (–) +9 bp	coding (1040‑1048/1209 nt)	tdcD ←	propionate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3262613 =	0 (0.000)	56 (1.320)	53/492	0.1	100%	coding (1048/1209 nt)	tdcD	propionate kinase
	?	NC_000913	3583428 =	NA (NA)					noncoding (1/768 nt)	IS1	repeat region
*	?	NC_000913	= 3262621	0 (0.000)	48 (1.130)	47/492	0.1	100%	coding (1040/1209 nt)	tdcD	propionate kinase
	?	NC_000913	= 3584195	NA (NA)					noncoding (768/768 nt)	IS1	repeat region

GTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTA  >  NC_000913/3583946‑3584436                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  gTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTAGCAGCTGTATGAATAGCGCCTGAAGGGAAAGCTGCACGTAATCAGCACGAGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGCCGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCg                                                                                                                                                                                                                                                  >  1:409980/1‑251 (MQ=11)   cACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGATGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCAGTATCTGAACATAAAACGCTATCAATAAGTTGGAGTCATTACCGa                                                                                                                                                                                                                                                 <  2:210401/250‑1 (MQ=11)    ccccTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGGCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTTCCGACa                                                                                                                                                                                                                                               <  2:187338/250‑1 (MQ=14)      ccTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGCTATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATACGTTGGAGTCATTACCGACAtt                                                                                                                                                                                                                                             >  2:403463/1‑251 (MQ=17)            aCGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGAccatgtttat                                                                                                                                                                                                                                       <  1:167989/251‑10 (MQ=9)                      atGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                                                                                                                                                                                                                              <  1:34977/250‑1 (MQ=39)                        ggATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCg                                                                                                                                                                                                                                                  <  2:324659/228‑1 (MQ=11)                        ggATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCg                                                                                                                                                                                                                                                  >  1:324659/1‑228 (MQ=11)                        ggATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                                                                                                                                                                                                                            <  1:293814/250‑1 (MQ=40)                         gATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaa                                                                                                                                                                                                                           >  1:200732/1‑250 (MQ=40)                                    gctggcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat                                                                                                                                                                                                               <  1:143172/251‑1 (MQ=255)                                        gcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTTGGACGGAAGTCGCTGGCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAGTAGAGATAATATTg                                                                                                                                                                                                           <  2:164772/251‑1 (MQ=255)                                         ccGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg                                                                                                                                                                                                          <  1:10460/251‑1 (MQ=255)                                          cGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                                                                                                                                                                                                         <  1:57381/251‑1 (MQ=255)                                             tGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcac                                                                                                                                                                                                       >  2:216155/1‑250 (MQ=255)                                                  gAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaaa                                                                                                                                                                                                  >  2:146297/1‑250 (MQ=255)                                                    aTCCCGACTGAAGGGAACGCTGCACGTAATCAGCAAGCGATATACGCCGCGAATTGAGCGGCATAACATGAATCTAAGACAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGTGGGAGCTGCATGACACAGTCCTCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAACAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaacat                                                                                                                                                                                                >  2:113918/1‑247 (MQ=17)                                                        cGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCGGTCGTACTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACTCTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt                                                                                                                                                                                           <  2:39965/251‑1 (MQ=255)                                                         gCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACTGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTTAACATAAAACACTATCAATAAGTGGGAGTCATTACCGACATTAAAAGTAACCTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTc                                                                                                                                                                                          <  1:239349/251‑1 (MQ=255)                                                                 ggAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCGCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGGGGGAGTCATTACCGACAATAAAAGTAACTTTTTGGAAATTAGAGATAATATTGGCACAGAAAAGATATTCAGTTCGtt                                                                                                                                                                                  <  2:77784/251‑1 (MQ=255)                                                                  gAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGttt                                                                                                                                                                                 >  1:300844/1‑251 (MQ=255)                                                                   aaaGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTa                                                                                                                                                                                <  1:298620/251‑1 (MQ=255)                                                                           aCGTAATCAGCAAGCGATATACGAAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATGGGCACAGAAAATATATTCAGTTCGTTTAGTGATTtc                                                                                                                                                                        <  2:410095/251‑1 (MQ=255)                                                                                   aGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCGGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCAtt                                                                                                                                                                <  2:145415/251‑1 (MQ=255)                                                                                         cGATATACGCAGCTAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGTCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGTTGGAGCTTCATGACAAAGTCATCGGGCATTTTCTGAACATAAAACACTATCAATAAGTTGGAGTTATTACCGGCATTATAAGTAACTTTTTTGAAATTAGAGTTAATAGTTGTACATAAAGTTTTTTCAGTTCGTTTAGTGATTTCTCTTCATTCTTaaa                                                                                                                                                          <  1:214120/251‑1 (MQ=17)                                                                                                 gCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCa                                                                                                                                                  <  1:296729/251‑1 (MQ=255)                                                                                                                 cATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCata                                                                                                                                  <  2:403615/251‑1 (MQ=255)                                                                                                                  aTAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATTGACAGTTGTTCataa                                                                                                                                 <  1:355023/251‑1 (MQ=255)                                                                                                                          aatttgagactgcaCTTGGCACGGCTGGGACGGAAGTCGCTGTCGATCTCAAAATCGGTGGAGCTCCCAGAAAACGTCATCGGGCTGTAACTGAACATAAAACACTATCAATAAGAGGGAGTCATTAACGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACCGAAAATATATTCAGTTCGTTTATTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCCTAATAATCat                                                                                                                          >  2:309278/12‑249 (MQ=25)                                                                                                                               gAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCaaaa                                                                                                                        <  1:101062/247‑1 (MQ=255)                                                                                                                               gAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAg                                                                                                                       <  1:308808/248‑1 (MQ=255)                                                                                                                               gAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGa                                                                                                                     <  1:291788/250‑1 (MQ=255)                                                                                                                                  gcagcaACTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACa                                                                                                                 <  2:203019/251‑1 (MQ=255)                                                                                                                                      caCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACaaaa                                                                                                              >  1:112121/1‑250 (MQ=255)                                                                                                                                                     ggACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATg                                                                                              >  1:187894/1‑251 (MQ=255)                                                                                                                                                       aCGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGcc                                                                                            >  1:342842/1‑251 (MQ=255)                                                                                                                                                          gAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATc                                                                                         <  2:342842/251‑1 (MQ=255)                                                                                                                                                                       ttctccacaTCGGAGGCGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGAAACTTTTTTGAAATTAGAGATAATATTGGCCCCGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGCCAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGAtt                                                                             >  2:138531/9‑250 (MQ=38)                                                                                                                                                                           cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAGCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTAcc                                                                          <  1:58276/249‑1 (MQ=255)                                                                                                                                                                                                  aaaGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTTTAAACCATGCCATCTTTTTAGATATTACCAGCAGAAAATTCCTTGTTCGAAAtt                                                 <  1:371114/251‑1 (MQ=255)                                                                                                                                                                                                         ccGGGCATTATCTGAACATAAATCCCTATCAATAAGCTGGAGTCATTAACGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAATGGCACAGAAAATATAGTCATTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATGGTAAACCATGCCATCTTTTTAGAGAGTACCAGCAGAAAATTCCTTGTTCGACATTTAGCAGt                                          <  2:161385/250‑1 (MQ=255)                                                                                                                                                                                                                 tATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTACCTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCt                                   >  1:370220/1‑250 (MQ=255)                                                                                                                                                                                                                                cacTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACa                    >  2:153593/1‑250 (MQ=255)                                                                                                                                                                                                                                             ttGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTc      >  2:179340/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gTCATTACCGACATTAAAAGTAAATTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTa  >  2:29601/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  GTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTA  >  NC_000913/3583946‑3584436

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