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

TGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATT  >  NC_000913/3822401‑3822894                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          tagttgaaagcctgaatcaacagattcaaacctactggccggaagaccaaatcaagcggctgcggccggcgtatctcgttgcacgtgatgctcacgaggggcaaacccgttcaagcggtgaaccctttatctcgcgcccggtagcggttgcctgcattctggacgagatgatactcgactatgaacctctgtttgcggcgctgctgcatgacgtgattgaagatactcccgccacctaccagtatcagga                                                                                                                                                                                                                                                       >  2:31074‑M1/246‑250 (MQ=255)       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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          TGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTG‑CAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATT  >  NC_000913/3822401‑3822894

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.