Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	3,993,337	Δ11 bp	100%	coding (2185‑2195/2547 nt)	cyaA →	adenylate cyclase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3993337	3993347	11	100 [0]	[0] 99	cyaA	adenylate cyclase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3993336	0 (0.000)	90 (1.060)	69/264	NT	100%	coding (2184/2547 nt)	cyaA	adenylate cyclase
	?	NC_000913	3993348 =	0 (0.000)					coding (2196/2547 nt)	cyaA	adenylate cyclase

CAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGT  >  NC_000913/3993202‑3993495                                                                                                                                                   |                                                                                                                                                    gcgtttctggtcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGt                                                                                                                                                   <  4:607430‑M2/3‑1 (MQ=255)   gtttctggtcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTgg                                                                                                                                                 <  2:241514‑M2/5‑1 (MQ=255)       ctggtcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACgg                                                                                                                                             >  3:586748‑M2/141‑149 (MQ=255)         ggtcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCt                                                                                                                                           <  1:432321‑M2/11‑1 (MQ=255)          gtcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCtt                                                                                                                                          >  2:318729‑M2/138‑149 (MQ=255)          gtcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCtt                                                                                                                                          >  3:717284‑M2/138‑149 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaTGTGGACGGCttt                                                                                                                                         >  4:517794‑M2/137‑149 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTg                                                                                                                                                  <  3:741960‑M2/4‑1 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTg                                                                                                                                                  >  4:741960‑M2/137‑140 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCttt                                                                                                                                         <  4:705565‑M2/13‑1 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCttt                                                                                                                                         <  1:165074‑M2/13‑1 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCttt                                                                                                                                         <  1:390282‑M2/13‑1 (MQ=255)           tcaaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCttt                                                                                                                                         >  4:281213‑M2/137‑149 (MQ=255)             aaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGc                                                                                                                                       >  1:202484‑M2/135‑149 (MQ=255)             aaacctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGc                                                                                                                                       <  4:499476‑M2/15‑1 (MQ=255)                 ctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAg                                                                                                                                    <  4:23235‑M2/18‑1 (MQ=255)                 ctgggggttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAg                                                                                                                                    >  3:23235‑M2/131‑148 (MQ=255)                       gttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGGCGGCTTTGCCAGCGAAggg                                                                                                                             <  3:598416‑M2/25‑1 (MQ=255)                       gttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAggg                                                                                                                             <  4:505901‑M2/25‑1 (MQ=255)                        ttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                                                            <  1:394199‑M2/26‑1 (MQ=255)                        ttgttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                                                            <  4:276975‑M2/26‑1 (MQ=255)                          gttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGatc                                                                                                                          <  1:517162‑M2/28‑1 (MQ=255)                          gttcttcgaacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGatc                                                                                                                          >  4:747366‑M2/122‑149 (MQ=255)                                  aacgcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGtt                                                                                                                  >  1:572829‑M2/114‑149 (MQ=255)                                     gcctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCtt                                                                                                               <  1:349073‑M2/39‑1 (MQ=255)                                       ctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGatc                                                                                                                          <  2:671723‑M2/28‑1 (MQ=255)                                       ctgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGatc                                                                                                                          >  1:671723‑M2/109‑136 (MQ=255)                                        tgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                                                            <  3:709518‑M2/42‑1 (MQ=255)                                        tgaatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                                                            <  4:239358‑M2/42‑1 (MQ=255)                                          aatgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgttaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCg                                                                                                          <  3:517794‑M2/44‑1 (MQ=255)                                            tgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaa                                                                                                        <  4:54419‑M2/46‑1 (MQ=255)                                            tgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaa                                                                                                        <  2:707889‑M2/46‑1 (MQ=255)                                            tgtatcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaa                                                                                                        <  4:216759‑M2/46‑1 (MQ=255)                                                tcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGaaa                                                                                                    <  4:510074‑M2/50‑1 (MQ=255)                                                tcggtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGAACATCCAGTTCTTTTTCGAAGaaa                                                                                                    <  4:355794‑M2/50‑1 (MQ=255)                                                   gtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttaaaaccaatcacgtcaaag                                                                                                                                                     <  4:596515‑M2/1‑1 (MQ=255)                                                   gtacagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttaaaaccaatcacgtcaaag                                                                                                                                                     >  3:596515‑M2/97‑97 (MQ=255)                                                     acagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCaa                                                                                               >  3:430807‑M2/95‑149 (MQ=255)                                                     acagaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCaa                                                                                               <  2:693130‑M2/55‑1 (MQ=255)                                                       agaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGa                                                                                             >  4:24305‑M2/93‑149 (MQ=255)                                                       agaaactagaaaacgcaatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgcaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGa                                                                                             <  4:558235‑M2/57‑1 (MQ=255)                                                        gaaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGAc                                                                                            <  2:572829‑M2/58‑1 (MQ=255)                                                         aaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                                                            >  4:588150‑M2/91‑116 (MQ=255)                                                         aaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                                                            <  3:588150‑M2/26‑1 (MQ=255)                                                         aaactggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACg                                                                                           >  1:591605‑M2/91‑149 (MQ=255)                                                            ctggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACgaga                                                                                        >  4:189247‑M2/88‑149 (MQ=255)                                                             tggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCt                                                                                                                                           >  2:64087‑M2/87‑97 (MQ=255)                                                             tggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCt                                                                                                                                           <  1:64087‑M2/11‑1 (MQ=255)                                                              ggaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAt                                                                                      >  1:324669‑M2/86‑149 (MQ=255)                                                               gaaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATg                                                                                     >  2:152299‑M2/85‑149 (MQ=255)                                                                aaaacgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATgg                                                                                    >  2:465728‑M2/84‑149 (MQ=255)                                                                   acgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCtt                                                                                 <  2:695361‑M2/69‑1 (MQ=255)                                                                    cgccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCttt                                                                                <  1:152299‑M2/70‑1 (MQ=255)                                                                      ccatcgagttttatggctcgatttcgcataacaagctgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTaa                                                                              <  1:72098‑M2/72‑1 (MQ=255)                                                                      ccatcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTaa                                                                              <  1:318729‑M2/72‑1 (MQ=255)                                                                       catcgagttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAat                                                                             <  1:560958‑M2/73‑1 (MQ=255)                                                                              ttttatggcgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACa                                                                      >  1:695780‑M2/70‑149 (MQ=255)                                                                                      cgcgatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgac                                                              >  2:540762‑M2/62‑149 (MQ=255)                                                                                         gatttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                                           <  3:592595‑M2/91‑1 (MQ=255)                                                                                          atttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                                          <  4:714027‑M2/92‑1 (MQ=255)                                                                                          atttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                                          <  3:747366‑M2/92‑1 (MQ=255)                                                                                           tttcgcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAg                                                         >  4:11014‑M2/57‑149 (MQ=255)                                                                                               gcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtc                                                                   <  1:179515‑M2/83‑1 (MQ=255)                                                                                               gcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtc                                                                   >  2:179515‑M2/53‑135 (MQ=255)                                                                                               gcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAc                                                     <  3:281213‑M2/97‑1 (MQ=255)                                                                                               gcataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAc                                                     >  4:525135‑M2/53‑149 (MQ=255)                                                                                                 ataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                                          >  1:131363‑M2/51‑142 (MQ=255)                                                                                                 ataacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                                          <  2:131363‑M2/92‑1 (MQ=255)                                                                                                  taacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCgg                                                  <  2:528835‑M2/100‑1 (MQ=255)                                                                                                   aacaaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCggg                                                 <  2:560256‑M2/101‑1 (MQ=255)                                                                                                     caaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtctc                                                                 >  3:663650‑M2/47‑131 (MQ=255)                                                                                                     caaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtctc                                                                 <  4:663650‑M2/85‑1 (MQ=255)                                                                                                     caaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGtt                                               >  1:757374‑M2/47‑149 (MQ=255)                                                                                                      aaactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTg                                              >  2:735905‑M2/46‑149 (MQ=255)                                                                                                       aactgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGa                                             <  1:25406‑M2/105‑1 (MQ=255)                                                                                                          tgcacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGt                                          <  4:294652‑M2/108‑1 (MQ=255)                                                                                                            cacggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGtat                                        <  2:324669‑M2/110‑1 (MQ=255)                                                                                                             acggcctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGtata                                       <  4:27025‑M2/111‑1 (MQ=255)                                                                                                                gcctgtcagtgcaggttgaaaccaatcacgtcaaaGGT‑GACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcac                                   >  3:654831‑M2/36‑149 (MQ=255)                                                                                                                  ctgtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacc                                  >  2:473625‑M2/34‑149 (MQ=255)                                                                                                                    gtcagtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcaccac                                <  3:227560‑M2/118‑1 (MQ=255)                                                                                                                        gtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGaa                                                                                       <  4:592773‑M2/63‑1 (MQ=255)                                                                                                                        gtgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGaa                                                                                       >  3:592773‑M2/28‑90 (MQ=255)                                                                                                                         tgcaggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGa                           >  1:382103‑M2/27‑149 (MQ=255)                                                                                                                            aggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAt                                                                                      <  4:720318‑M2/64‑1 (MQ=255)                                                                                                                            aggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAt                                                                                      >  3:720318‑M2/24‑87 (MQ=255)                                                                                                                            aggttgaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAgg                        <  2:382103‑M2/126‑1 (MQ=255)                                                                                                                                 gaaaccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgcagca                   >  1:768388‑M2/19‑149 (MQ=255)                                                                                                                                    accaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCtttt                                                                                                             >  2:153094‑M2/16‑56 (MQ=255)                                                                                                                                    accaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCtttt                                                                                                             <  1:153094‑M2/41‑1 (MQ=255)                                                                                                                                     ccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAaga               <  1:428781‑M2/135‑1 (MQ=255)                                                                                                                                     ccaatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAaga               <  2:143070‑M2/135‑1 (MQ=255)                                                                                                                                      caatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAagag              >  1:491100‑M2/14‑149 (MQ=255)                                                                                                                                      caatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAagag              >  3:176933‑M2/14‑149 (MQ=255)                                                                                                                                       aatcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAgagg             >  2:460754‑M2/13‑149 (MQ=255)                                                                                                                                        atcacgtcaaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAgagga            <  2:757374‑M2/138‑1 (MQ=255)                                                                                                                                                aaaGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTAGATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTAc    <  4:717284‑M2/146‑1 (MQ=255)                                                                                                                                                  aGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATGTCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGt  <  1:506333‑M2/148‑1 (MQ=255)                                                                                                                                                  aGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGt  <  3:24305‑M2/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGT  >  NC_000913/3993202‑3993495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑