Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	3,822,651	(CAGGAT)1→2	coding (252/2109 nt)	spoT →	bifunctional (p)ppGpp synthase/hydrolase SpoT

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3822646 =	0 (0.000)	42 (1.580)	40/490	0.0	98.8%	coding (247/2109 nt)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
	?	NC_000913	= 3822651	1 (0.040)					coding (252/2109 nt)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT

TGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAG  >  NC_000913/3822407‑3822899                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          tagttgaaAGCCTGAATCAACAGATTCAAACCTACTGGCCGGAAGACCAAATCAAGCGGCTGCGGCCGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACCCGTTCAAGCGGTGAACCCTTTATCTCGCGCCCGGTAGCGGTTGCCTGCATTCTGGACGAGATGATACTCGACTATGAACCTCTGTTTGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGTATcagga                                                                                                                                                                                                                                                      >  2:31074‑M2/4‑245 (MQ=255)     cctgaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGc                                                                                                                                                                                                                                                 >  1:208885‑M1/236‑251 (MQ=255)      ctgaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCt                                                                                                                                                                                                                                                <  2:228761‑M1/17‑1 (MQ=255)          atcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCttttt                                                                                                                                                                                                                                            >  1:134639‑M1/231‑251 (MQ=255)           ccTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagc                                                                                                                                                                                                                                           >  1:208885‑M2/1‑235 (MQ=255)            cTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagct                                                                                                                                                                                                                                          <  2:228761‑M2/251‑18 (MQ=255)                aTCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagcttttt                                                                                                                                                                                                                                      >  1:134639‑M2/1‑230 (MQ=255)                     caaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGt                                                                                                                                                                                                                                 <  1:225461‑M1/32‑1 (MQ=255)                         cctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGcc                                                                                                                                                                                                                             <  1:242244‑M1/36‑1 (MQ=255)                           cAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgt                                                                                                                                                                                                                           <  1:225461‑M2/251‑33 (MQ=255)                               cctacctGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgcc                                                                                                                                                                                                                       <  1:242244‑M2/251‑37 (MQ=255)                                                    gtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACtt                                                                                                                                                                                                  <  2:36970‑M1/63‑1 (MQ=255)                                                          gTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaactt                                                                                                                                                                                            <  2:36970‑M2/251‑64 (MQ=255)                                                               gcgtatctcgttgcacgtgaggctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggctctgctgcatgacgtgattgaagatactcccgcctcctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCaa                                                                                                                                                                                       >  1:109225‑M1/178‑251 (MQ=255)                                                                 gtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGt                                                                                                                                                                                     <  1:122806‑M1/76‑1 (MQ=255)                                                                     gCGTATCTCGTTGCACGTGAGGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCTCTGCTGCATGACGTGATTGAAGATACTCCCGCCTCCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaa                                                                                                                                                                                 >  1:109225‑M2/1‑177 (MQ=255)                                                                       gTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagt                                                                                                                                                                               <  1:122806‑M2/251‑77 (MQ=255)                                                                                tgatgctcacgaggggcaaacaagttcaagccgtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgcaccgctgatggcggcgctgctgcatgacgtgcttgaagatactcccgccacctaccaggagCAGGATATGGAACAGCTTTTTGGTAACCGCGTCGCCGAGCTGGTAGAGGGGGTGTCGCACCTTGATAAACTCAAGTTCCGcgatacgac                                                                                                                                                                        <  2:41555‑M1/89‑5 (MQ=255)                                                                                 gatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaagaaaga                                                                                                                                                                     >  1:228761‑M1/160‑251 (MQ=255)                                                                                     ctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcaaccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcaagacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGCAGAGGGGGAGTCGAACCTTGATAAACTCAAGTTCCGCGATAAGAAAGAGgcg                                                                                                                                                                 <  1:225918‑M1/96‑1 (MQ=255)                                                                                      tGATGCTCACGAGGGGCAAACAAGTTCAAGCCGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGCACCGCTGATGGCGGCGCTGCTGCATGACGTGCTTGAAGATACTCCCGCCACCTACCAGGAGcaggatatggaacagctttttggtaaccgcgtcgccgagctggtagagggggtgtcgcaccttgataaactcaagttccgcgatacgac                                                                                                                                                                  <  2:41555‑M2/249‑90 (MQ=255)                                                                                       gATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaaga                                                                                                                                                               >  1:228761‑M2/1‑159 (MQ=255)                                                                                       cacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGgcg                                                                                                                                                                 >  1:148204‑M1/154‑249 (MQ=255)                                                                                           cTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCAACCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCAAGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggcagagggggagtcgaaccttgataaactcaagttccgcgataagaaagaggcg                                                                                                                                                           <  1:225918‑M2/251‑97 (MQ=255)                                                                                             gggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCg                                                                                                                                                          <  1:160175‑M1/103‑1 (MQ=255)                                                                                             cACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcg                                                                                                                                                           >  1:148204‑M2/1‑153 (MQ=255)                                                                                                   gggCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccg                                                                                                                                                    <  1:160175‑M2/250‑104 (MQ=255)                                                                                                     acgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTc                                                                                                                                                 <  2:16932‑M1/112‑1 (MQ=255)                                                                                                      cgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgacgatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccaccttccagtatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAGCTTTCg                                                                                                                                                >  1:166455‑M1/139‑251 (MQ=255)                                                                                                       gttcaaccggtgcaccctatatcacgcacccggtagcgggtgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGCGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCAATAAGAAAGAGGCGCAGGCCGAAAACTTTCGc                                                                                                                                               <  2:28320‑M1/114‑1 (MQ=255)                                                                                                          caagcggtgaaccctatatcacgcacccggtagcgggtgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTGCGCAAg                                                                                                                                            >  2:244450‑M1/135‑251 (MQ=255)                                                                                                           acGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttc                                                                                                                                           <  2:16932‑M2/251‑113 (MQ=255)                                                                                                            cGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACGATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTTCCAGTATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaagctttcg                                                                                                                                          >  1:166455‑M2/1‑138 (MQ=255)                                                                                                            agcggtgaaccctatatcacgctcccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgat                                                                                                                                          >  1:122590‑M1/133‑251 (MQ=255)                                                                                                             gTTCAACCGGTGCACCCTATATCACGCACCCGGTAGCGGGTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgcgctggtagagggggtgtcgaaacttgataaactcaagttccgcaataagaaagaggcgcaggccgaaaactttcgc                                                                                                                                         <  2:28320‑M2/251‑115 (MQ=255)                                                                                                              cggtgaaccctatatcacgcacccggtagcggttgcctgcattctggcctagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccatctaccaggatCAGGAGATGGAACAGTTTTTTCGTAAAAGCGTCGGCGAGCTGTTAGAGGTGTTTTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgatga                                                                                                                                        >  1:41555‑M1/131‑251 (MQ=255)                                                                                                                cAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGGTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaacttgcgcaag                                                                                                                                      >  2:244450‑M2/1‑134 (MQ=255)                                                                                                                  aGCGGTGAACCCTATATCACGCTCCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagat                                                                                                                                    >  1:122590‑M2/1‑132 (MQ=255)                                                                                                                    cGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCTAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCATCTACCAGGATcaggagatggaacagttttttcgtaaaagcgtcggcgagctgttagaggtgttttcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatga                                                                                                                                  >  1:41555‑M2/1‑130 (MQ=255)                                                                                                                             cacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTg                                                                                                                          >  1:254114‑M1/116‑250 (MQ=255)                                                                                                                                   cACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtg                                                                                                                    >  1:254114‑M2/1‑115 (MQ=255)                                                                                                                                    ccggtatcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAATATGATTATGGCGATGGTG‑CAGGatat                                                                                                                 >  2:27755‑M1/109‑251 (MQ=255)                                                                                                                                     cggtaggggttgcctgcagtctggccgagatgaaactcgactatgaaacgttgatgtcggcgcttctgcatgacgtgattgaagatactcccgccacctaccagtatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACGTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCTGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑TTGGATATc                                                                                                                >  2:211213‑M1/108‑251 (MQ=255)                                                                                                                                      ggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATcc                                                                                                               >  1:36970‑M1/107‑251 (MQ=255)                                                                                                                                          ccGGTATCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaatatgattatggcgatggtgcaggatat                                                                                                            >  2:27755‑M2/1‑108 (MQ=255)                                                                                                                                           cGGTAGGGGTTGCCTGCAGTCTGGCCGAGATGAAACTCGACTATGAAACGTTGATGTCGGCGCTTCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGTATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacgtgataaactcaagttccgcgataagaaagaggcgctggccgaaaactttcgcaagatgattatggcgatggtgttggatatc                                                                                                           >  2:211213‑M2/1‑107 (MQ=255)                                                                                                                                            ggTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatcc                                                                                                          >  1:36970‑M2/1‑106 (MQ=255)                                                                                                                                                            gccgagctgaaactcgactatgacacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgacacctaccaggatCAGGCTATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGCGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGc                                                                                         <  1:244734‑M1/167‑1 (MQ=255)                                                                                                                                                              cgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCg                                                                                       <  2:133985‑M1/169‑1 (MQ=255)                                                                                                                                                                 gatgaaactcgactatgaaacgctgatggcgacgctgctgcatgacgtgattgaagatactcccgccacctcccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGCGTCGAAACTTGATAAACTCAAGTTCCGCGCTAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGAcc                                                                                    <  2:33064‑M1/172‑1 (MQ=255)                                                                                                                                                                  gCCGAGCTGAAACTCGACTATGACACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGACACCTACCAGGATcaggctatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagcggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgc                                                                                    <  1:244734‑M2/251‑168 (MQ=255)                                                                                                                                                                    cGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccg                                                                                  <  2:133985‑M2/251‑170 (MQ=255)                                                                                                                                                                     aaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGGCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAATATGATTATGGCGA‑GGTGGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTAc                                                                                >  2:114704‑M1/76‑251 (MQ=255)                                                                                                                                                                       gaTGAAACTCGACTATGAAACGCTGATGGCGACGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTCCCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggcgtcgaaacttgataaactcaagttccgcgctaagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgacc                                                                               <  2:33064‑M2/251‑173 (MQ=255)                                                                                                                                                                           aaaCTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtggcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaatatgattatggcgaggtggcaggatatccgcgtcatcctcatcaaacttgccgaccgtac                                                                           >  2:114704‑M2/1‑75 (MQ=255)                                                                                                                                                                            actatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCacaac                                                                         >  2:242244‑M1/69‑251 (MQ=255)                                                                                                                                                                             ctatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCacaac                                                                         <  1:16914‑M1/183‑1 (MQ=255)                                                                                                                                                                               atgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCacaaca                                                                        >  2:50631‑M1/66‑249 (MQ=255)                                                                                                                                                                                  aCTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaac                                                                    >  2:242244‑M2/1‑68 (MQ=255)                                                                                                                                                                                   cTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaac                                                                    <  1:16914‑M2/250‑184 (MQ=255)                                                                                                                                                                                     aTGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaaca                                                                   >  2:50631‑M2/1‑65 (MQ=255)                                                                                                                                                                                         gatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTgg                                                            >  1:253889‑M1/56‑251 (MQ=255)                                                                                                                                                                                               ggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCAc                                                      >  1:234270‑M1/50‑251 (MQ=255)                                                                                                                                                                                               gATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgg                                                       >  1:253889‑M2/1‑55 (MQ=255)                                                                                                                                                                                                     ggcgCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcac                                                 >  1:234270‑M2/1‑49 (MQ=255)                                                                                                                                                                                                      ctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCAc                                                                 <  1:170757‑M1/191‑1 (MQ=255)                                                                                                                                                                                                      ctgcatgacgagattgacgatactcacgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTCAAAGCGTCGCCGAGCTGGCAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCAGCGATCAGAAAGAGGCGAAGGCAGAAAACTTTCGCAAGATGATAATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCCAAACATGCGCACGCCGGGCTCCCTTCGccc                                               <  2:128055‑M1/209‑1 (MQ=255)                                                                                                                                                                                                      cggcatgacgtgagtgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCAc                                                                 >  2:170757‑M1/43‑233 (MQ=255)                                                                                                                                                                                                            ctgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcac                                                            <  1:170757‑M2/233‑192 (MQ=255)                                                                                                                                                                                                            ctgcATGACGAGATTGACGATACTCACGCCACCTACCAGGATcaggatatggaacagctttttggtcaaagcgtcgccgagctggcagagggggtgtcgaaacttgataaactcaagttcagcgatcagaaagaggcgaaggcagaaaactttcgcaagatgataatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtaccccaaacatgcgcacgccgggctcccttcgccc                                          <  2:128055‑M2/251‑210 (MQ=255)                                                                                                                                                                                                            cggcatGACGTGAGTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcac                                                            >  2:170757‑M2/3‑42 (MQ=255)                                                                                                                                                                                                                     gaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTcgccg                                <  2:117211‑M1/224‑1 (MQ=255)                                                                                                                                                                                                                           gAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccg                           <  2:117211‑M2/251‑225 (MQ=255)                                                                                                                                                                                                                                   cacctgccaggctcAGGATATGGATCAGCTTCTTGGTAAAAGCGTCGCCGAGCTGGTATAGGGGGTGTCGAAACGTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACGTTCGCAAGATGATTAGGGCGATGGTG‑CAGGATAGCCGCGTCAGCCTCTGCAAACTTGCCGACCGTACCCGCTACATGCGCACGCAGGGCTCTCTTCGCCCGGAGAAACTTCGCCGCGTCGGCCGTGaaa                  <  2:226216/238‑1 (MQ=255)                                                                                                                                                                                                                                           aggatCCGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCACGATGATTATGGCGATGGTG‑CAGGAAATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCACCCGGACAAACGACGCCTCATCGCCCGGGAAACTCTcgt           <  2:104866‑M1/245‑2 (MQ=255)                                                                                                                                                                                                                                            ggatCAGGATATGGCACAGCTTTTTGGTAACAGCATCGCCGCGCTGGTAGAGGCGGCGTCGCAACTTGAAAACCTCACGTTCCGCGATAAGAAAGAGGCGAAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACACACGCCGCCGCATCGCCCGTGACACTCTCGaa          <  2:120813‑M1/246‑1 (MQ=255)                                                                                                                                                                                                                                            ggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGGCGCCGCATCGCCCGTGAaactctccaa          <  2:213441‑M1/246‑4 (MQ=255)                                                                                                                                                                                                                                             gatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACGCAAGTGCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAACCGTCGCCGCATCGCCCGTGAAACTCTCGAAAt        >  1:133985‑M1/4‑251 (MQ=255)                                                                                                                                                                                                                                              atCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACGTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAAtt       >  2:186796‑M1/3‑251 (MQ=255)                                                                                                                                                                                                                                                 aggatccggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcacgatgattatggcgatggtgcaggaaatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcacccggacaaacgacgcctcatcgcccgggaaactctcgt      <  2:104866‑M2/250‑246 (MQ=255)                                                                                                                                                                                                                                                  ggatcaggatatggcacagctttttggtaacagcatcgccgcgctggtagaggcggcgtcgcaacttgaaaacctcacgttccgcgataagaaagaggcgaaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacacacgccgccgcatcgcccgtgacactctcgaa     <  2:120813‑M2/250‑247 (MQ=255)                                                                                                                                                                                                                                                  ggatcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacggcgccgcatcgcccgtgaaactctccaa     <  2:213441‑M2/250‑247 (MQ=255)                                                                                                                                                                                                                                                   gatcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaacgcaagtgccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaccgtcgccgcatcgcccgtgaaactctcgaaat   >  1:133985‑M2/1‑3 (MQ=255)                                                                                                                                                                                                                                                    atcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacgtgccgaccgtacccacaacatgcgcacgctgggctcacttcgcccggacaaacgtcgccgcatcgcccgtgaaactctcgaaatt  >  2:186796‑M2/1‑2 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          TGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAG  >  NC_000913/3822407‑3822899

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