Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	528,897	IS5 (+) +4 bp	42.5%	pseudogene (258‑261/491 nt)	ybbD →	pseudogene

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	273955 =	NA (NA)	9 (0.070)	7/544	NT	44.6%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	528897 =	11 (0.080)					pseudogene (258/491 nt)	ybbD	pseudogene
*	?	NC_000913	= 275149	NA (NA)	8 (0.060)	8/552	NT	39.6%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 528900	12 (0.090)					pseudogene (261/491 nt)	ybbD	pseudogene

TATATTAATCCAGAGCATAAAAAAACGATAAAGGCTTTACAAGATGCTGGCTATCTCAAGTAATCTATCAAAGATGATAATATTTATTTTTGCTATTATAATCATTGTTGTTTTATGCGTAATTACTTATCTTTATTTATACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAAAGTAAAATTCC  >  NC_000913/528597‑529197                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                             ccacccattactcatcaccaagcgacagttcgtcttcgtgaaagccagatacaaggggttgctgaaaaacgataaacaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                                                                                                                                                                              <  1:798214‑M2/1‑1 (MQ=255)  cacccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                                                                                                                                                                              <  2:941599‑M2/1‑1 (MQ=255)  cacccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                                                                                                                                                                              <  2:658509‑M2/1‑1 (MQ=255)  cacccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggtggctgaaaaacgataaccaacttgcgatgttattcacgctggccaacctgtgtcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                                                                                                                                                                              <  1:771384‑M2/1‑1 (MQ=255)    cccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaactgaaatgactgagttagccgagaagaatttcccacgttattcgcacattccCTaa                                                                                                                                                                                                                                                                                                           >  1:190230‑M2/298‑301 (MQ=255)    cccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTaa                                                                                                                                                                                                                                                                                                           >  1:580026‑M2/298‑301 (MQ=255)    cccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTaa                                                                                                                                                                                                                                                                                                           >  2:1147170‑M2/298‑301 (MQ=255)    cccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgaaaagaatttccccgcttattcgcacattccCTaa                                                                                                                                                                                                                                                                                                           >  1:619679‑M2/298‑301 (MQ=255)    cccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactcagtcagccgagaagaatttccccgcttattcgcactttccCTaa                                                                                                                                                                                                                                                                                                           >  2:497988‑M2/298‑301 (MQ=255)    cccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacaggaataaatgaaatgactgagtcagccgagaagcattaccccgcttattcgcaccttccCTaa                                                                                                                                                                                                                                                                                                           >  1:181912‑M2/298‑301 (MQ=255)     ccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggtgttgctgaaaaacgataaccaactggcgattttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          <  1:666916‑M2/5‑1 (MQ=255)     ccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          <  2:740254‑M2/5‑1 (MQ=255)     ccatttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          <  2:249266‑M2/5‑1 (MQ=255)     ccatttcgaatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggtggctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaaggaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          <  1:838982‑M2/5‑1 (MQ=255)       atttcgcatcatcaagcgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACaa                                                                                                                                                                                                                                                                                                        <  2:166872‑M2/7‑1 (MQ=255)                       cgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACa                                                                                                                                                                                                                                                                                                         <  1:24881‑M2/6‑1 (MQ=255)                       cgacagttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACa                                                                                                                                                                                                                                                                                                         >  2:24881‑M2/279‑284 (MQ=255)                            gttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          <  2:998425‑M2/5‑1 (MQ=255)                            gttcggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaacttccctaacctgtcccttatacacctctgac                                                                                                                                                                                                                                                                                    >  1:998386‑M2/274‑278 (MQ=255)                               cggcttcgtgaaagccagatacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAATGATGATGGGGGTAGGTTTAAGa                                                                                                                                                                                                                                                                                <  1:25608‑M2/31‑1 (MQ=255)                                                 atacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTa                                                                                                                                                                                                                                                                                                            >  2:1109533‑M2/253‑255 (MQ=255)                                                 atacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          <  2:1070715‑M2/5‑1 (MQ=255)                                                 atacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAc                                                                                                                                                                                                                                                                                                          >  1:1070672‑M2/253‑257 (MQ=255)                                                 atacaaggggttgctgaaaaacgataaccaactggcgatggtattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTa                                                                                                                                                                                                                                                                                                            <  1:1109483‑M2/3‑1 (MQ=255)                                                  tacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                                                                                                                                                                              >  1:48304‑M2/252‑252 (MQ=255)                                                  tacaaggggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                                                                                                                                                                              <  2:48304‑M2/1‑1 (MQ=255)                                                                                                                     cggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAata                                                                                                                                                                                           <  2:180153‑M2/116‑1 (MQ=255)                                                                                                                                           ataCAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTTTAAATATGGTAAGATAACAGGGCAAgggg                                                                                                                                                                    >  2:127897/1‑301 (MQ=255)                                                                                                                                                  gATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACCGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgtgt                                                                                                                                                             >  2:412575/1‑301 (MQ=255)                                                                                                                                                                      aCATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAACAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTg                                                                                                                                         >  1:737874/1‑301 (MQ=255)                                                                                                                                                                                             ggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTAtttt                                                                                                                  >  2:855138‑M2/113‑301 (MQ=255)                                                                                                                                                                                                             ttcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATTTTGGAg                                                                                                  >  1:1069964‑M2/97‑301 (MQ=255)                                                                                                                                                                                                              gTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgtttta                                                                                                 <  2:599395‑M2/301‑203 (MQ=255)                                                                                                                                                                                                                         gCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaagg                                                                                         <  2:186455‑M2/298‑211 (MQ=255)                                                                                                                                                                                                                                             aaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg                                                                  <  1:855107‑M2/237‑1 (MQ=255)                                                                                                                                                                                                                                                  gCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGAATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTTCTTGa                                                             >  1:462972/1‑301 (MQ=255)                                                                                                                                                                                                                                                       ggATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATc                                                                                                                                                                                                                                  <  1:840574/131‑1 (MQ=255)                                                                                                                                                                                                                                                       ggATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATc                                                                                                                                                                                                                                  >  2:840604/1‑131 (MQ=255)                                                                                                                                                                                                                                                           gaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAagag                                                     <  2:1070007‑M2/250‑1 (MQ=255)                                                                                                                                                                                                                                                                aaTATACAAAAATGTAGAAGAGGATTATTGTTGTCTTATTTTTCCCTAACAAATGATGAGGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATc                                               <  1:127897/301‑1 (MQ=255)                                                                                                                                                                                                                                                                 atatACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGCTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg                                                                  >  1:888267/1‑281 (MQ=255)                                                                                                                                                                                                                                                                 atatACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGCTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg                                                                  <  2:888301/281‑1 (MQ=255)                                                                                                                                                                                                                                                                        caAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccgtcgaccttatccgctggaacccatgctacgcattca                                        <  1:982546‑M2/300‑260 (MQ=255)                                                                                                                                                                                                                                                                          aaTGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggaaaccatgctacgcattcactg                                     >  1:599380‑M2/1‑39 (MQ=255)                                                                                                                                                                                                                                                                           aTGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATTATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATAcc                                                                                                                                                                               <  2:855691/162‑1 (MQ=255)                                                                                                                                                                                                                                                                           aTGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGACAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAACAAGTGTGGTGTAAATATGGTAAGATAcc                                                                                                                                                                               >  1:855660/1‑162 (MQ=255)                                                                                                                                                                                                                                                                                     tgatTATTTTTTTCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggcaaccatgcaacgc                                             <  1:1071668‑M2/282‑255 (MQ=255)                                                                                                                                                                                                                                                                                     tgatTATTTTTTTCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggaaaccatgctacgc                                             >  2:1071711‑M2/1‑28 (MQ=255)                                                                                                                                                                                                                                                                                              ttttCTTATTTTTCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggaaaccatgctacgcattcactgcatgcagcattggtacaacc                 >  1:186453‑M2/1‑19 (MQ=255)                                                                                                                                                                                                                                                                                                     aTTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAATAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTTAGGAAAAga          >  1:612211/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                          tCCCTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggaaaccatgctacgcattcactgcatgcagcattggtacaacctgagcgatggcg     <  2:430877‑M2/301‑295 (MQ=255)                                                                                                                                                                                                                                                                                                            ccTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttacccgctggaaaccatgctacgcattcactgcatgcagcattggtacaacctgagcgatggctcg   <  2:703374‑M2/301‑297 (MQ=255)                                                                                                                                                                                                                                                                                                             cTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggaaaccatgctacgcattcactgcatgcagcattggtacaagctgagcgatggcgcgg  <  1:34859‑M2/301‑298 (MQ=255)                                                                                                                                                                                                                                                                                                             cTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccggcgaccttatccgctggaaaccatgctacgcattcactgcatgcagcattggtacaacctgagcgatggcgcgg  <  1:246083‑M2/301‑298 (MQ=255)                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                             TATATTAATCCAGAGCATAAAAAAACGATAAAGGCTTTACAAGATGCTGGCTATCTCAAGTAATCTATCAAAGATGATAATATTTATTTTTGCTATTATAATCATTGTTGTTTTATGCGTAATTACTTATCTTTATTTATACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAAAGTAAAATTCC  >  NC_000913/528597‑529197

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