Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,000,721	IS4 (–) +12 bp	coding (325‑336/765 nt)	yigE ←	DUF2233 domain‑containing protein YigE

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4000721 =	0 (0.000)	21 (0.890)	21/496	0.3	100%	coding (336/765 nt)	yigE	DUF2233 domain‑containing protein YigE
	?	NC_000913	= 4503515	NA (NA)					noncoding (1/1426 nt)	IS4	repeat region
*	?	NC_000913	= 4000732	0 (0.000)	18 (0.770)	17/490	0.5	100%	coding (325/765 nt)	yigE	DUF2233 domain‑containing protein YigE
	?	NC_000913	4502090 =	NA (NA)					noncoding (1426/1426 nt)	IS4	repeat region

ACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCA  >  NC_000913/4000474‑4000970                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          acacGGCGTTCACATGTTTATTAATCCCAACCCCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCATGACTGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAcc                                                                                                                                                                                                                                                        <  2:213242/251‑1 (MQ=255)        gTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                                                                                                                                                                                                                                 >  1:17902/1‑251 (MQ=255)        gTTCCCATGTTTATTAATCCCAACACCGTTACGAAGTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                                                                                                                                                                                                                                 <  2:111350/251‑1 (MQ=255)                     cggcgcttgcgtagcgttacagtacctgattcggcaaggcaacgagagatgagttcggggtcgaggtaatcccccagagaagtcagtgggttacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGAt                                                                                                                                                                                                                                    >  2:118718‑M1/228‑251 (MQ=255)                                       acagttcctgattcggcaaggcaacgagagatgagttcggggtcgaggtaatcccccagagaagtcagtgggttacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaacggatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggcg                                                                                                                                                                                                                   >  1:27388‑M1/210‑250 (MQ=255)                                                     ggCGACGTTGGGATGAATACGAGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggct                                                                                                                                                                                                    >  2:18683‑M1/1‑207 (MQ=255)                                                            caacgagagatgagttcggggtcgaggtaatcccccagagaagtcagtgggttacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgcc                                                                                                                                                                                             <  1:37572‑M1/63‑1 (MQ=255)                                                              acgagagatgagttcggggccgaggtaatcccccagagaagtcagtgggttacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcggttccggattcttacaccgccactggatcggtcaactgatccttaactgatcggcagtaCACCTTCTGCTGACCGTTTTAGATGTACAAAACGAGCCGCGCATAGCTTTCTTCATAGATgccgc                                                                                                                                                                                           <  2:180357‑M1/65‑1 (MQ=255)                                                                     aTACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagg                                                                                                                                                                                    >  1:108707‑M1/1‑191 (MQ=255)                                                                      tACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGATTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggccctttttcggggctgtggggtatttccaggg                                                                                                                                                                                   >  2:137958‑M1/1‑190 (MQ=255)                                                                       aCGCGGATTAATTACACAGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggc                                                                                                                                                                                  >  1:196422‑M1/1‑189 (MQ=255)                                                                                       ccGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccggtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactc                                                                                                                                                                  <  1:180523‑M1/251‑79 (MQ=255)                                                                                           cccccagagaagtcagtgggttacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctg                                                                                                                                                               >  1:231217‑M1/158‑250 (MQ=255)                                                                                            ccccagagaagtcagtgggttacgcagagaatcgtaacgggctaccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctgac                                                                                                                                                             <  2:192468‑M1/95‑1 (MQ=255)                                                                                                       ttGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttc                                                                                                                                                  >  1:73529‑M1/1‑157 (MQ=255)                                                                                                            gggttacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcag                                                                                                                                             >  1:117567‑M1/141‑251 (MQ=255)                                                                                                               gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCAGCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactcgcgggaaggccctttcccttctt                                                                                                                                           <  2:95858‑M1/250‑102 (MQ=255)                                                                                                                tacgcagagaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagc                                                                                                                                            <  2:235573‑M1/112‑1 (MQ=255)                                                                                                                        gaatcgtaacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATATGCACCTGACACTGACTATTAATATCCGCCAGCAGAGCATGTa                                                                                                                                  <  2:128679‑M1/122‑1 (MQ=255)                                                                                                                             cTGAATATCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcac                                                                                                                             >  1:123458‑M1/1‑135 (MQ=255)                                                                                                                               aacgggataccagatcaagagcctgtccaatgtgcataaaaaaatccggaaacaagtgagcgtttccggattcttccacagccactggatcggtcaactgatccttacctgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGAAACGTTccc                                                                                                                          <  1:137919‑M1/130‑1 (MQ=255)                                                                                                                                                  gAAGGCCTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccatttcccttcttgtcggtaatttcacaagttgtcccatacttgca                                                                                                          >  2:206178‑M1/1‑114 (MQ=255)                                                                                                                                                  gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcgggtctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcaccagttgtcccatacttgca                                                                                                          >  2:206174‑M1/1‑114 (MQ=255)                                                                                                                                                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagatcgcgc                                                                                                 <  1:6096‑M1/251‑144 (MQ=255)                                                                                                                                                         ccaatgtgcataaaaaaatccggaaacaagggagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCa                                                                                                >  1:18123‑M1/96‑251 (MQ=255)                                                                                                                                                                      aaaaatccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCAcc                                                                                   <  1:70929‑M1/169‑1 (MQ=255)                                                                                                                                                                            ccggaaacaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGctct                                                                             >  1:92939‑M1/77‑251 (MQ=255)                                                                                                                                                                                  acaagtgagcgtttccggattcttacacagccactggatcggtcaactgatccttaactgatcggccttaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCTGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTg                                                                        <  2:22315‑M1/180‑1 (MQ=255)                                                                                                                                                                                    cccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagatcgcgcatcagctccggtatacgtcccggtg                                                                        <  2:17902‑M1/248‑169 (MQ=255)                                                                                                                                                                                     ccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagagcgcgcatcagctccggtatacgtcccggtgaagc                                                                    <  2:154989‑M1/251‑173 (MQ=255)                                                                                                                                                                                             tttccggattcttacacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGt                                                            <  2:79326‑M1/192‑1 (MQ=255)                                                                                                                                                                                                          acacagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACg                                               <  1:116398‑M1/205‑1 (MQ=255)                                                                                                                                                                                                            acagccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGt                                             >  1:60129‑M1/45‑251 (MQ=255)                                                                                                                                                                                                                ccactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGc                                         >  1:188184‑M1/41‑251 (MQ=255)                                                                                                                                                                                                                 cactggatcggtcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCg                                        >  2:82555‑M1/40‑251 (MQ=255)                                                                                                                                                                                                                            tcaactgatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTg                               <  1:179973‑M1/221‑1 (MQ=255)                                                                                                                                                                                                                                  gatccttaactgatcggcattaCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcg                       <  1:237943‑M1/229‑1 (MQ=255)                                                                                                                                                                                                                                          gCGAGATTTAACGCCACCTTCTGCTGtaatgccagtcagttaagcaactgactggctcttttgcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagatcgcgcatcagttccggtatacgtcccggtgaagcgccctgcaatgtcatcagcattctcatcaccattccgcatgattctgagaaac               <  2:118997‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                                                 ttAACTCCTCCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttacctctctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagatcgcgcatcagctccggtatacgtccctgtgaagcgccctgcaatgtcatcagcattctcatcaccattccgcatgattctgagaaactcagtt         <  1:183158‑M1/250‑232 (MQ=255)                                                                                                                                                                                                                                                       ccACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagatcgcgcatcagctccggtatacgtcccggtgaagcgccctgcaatgtcatcagcattctcatcaccattccgcatgattctgagaaactcagttgattcgg  >  1:175436‑M1/1‑13 (MQ=255)                                                                                                                                                                                                                                                       ccACCTTCTGCTGtaatgccagtcagttaagcaactgactggctctttttcggggctgtggggtatttccagggcctctcctttaccactctcgggaaggccctttcccttcttgtcggtaatttcacaagttgtcccatacttgcaagatcgcgcatcagctccggtatacgtcccggtgaagcgccctgcaatgtcatcagcattctcatcaccattccgcatgattctgagaaactcagttgattcgg  <  1:115595‑M1/251‑239 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          ACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCA  >  NC_000913/4000474‑4000970

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