Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,991,548	Δ200 bp	coding (396‑595/2547 nt)	cyaA →	adenylate cyclase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3991548	3991747	200	61 [0]	[0] 62	cyaA	adenylate cyclase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3991547	0 (0.000)	61 (0.540)	43/598	1.7	100%	coding (395/2547 nt)	cyaA	adenylate cyclase
	?	NC_000913	3991748 =	0 (0.000)					coding (596/2547 nt)	cyaA	adenylate cyclase

AGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGC  >  NC_000913/3991250‑3991840                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                      agttctaCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtg                                                                                                                                                                                                                                                                                                     <  2:10431‑M1/297‑3 (MQ=255)          gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCNGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtgnaatatg                                                                                                                                                                                                                                                                                            >  2:752881‑M1/1‑289 (MQ=255)          gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTACGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatg                                                                                                                                                                                                                                                                                            >  2:545640‑M1/1‑289 (MQ=255)            cTACTGCCTACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggt                                                                                                                                                                                                                                                                                          <  2:428227‑M1/300‑14 (MQ=255)            cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggt                                                                                                                                                                                                                                                                                          <  1:829523‑M1/300‑14 (MQ=255)              aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACTTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgc                                                                                                                                                                                                                                                                                        <  1:143139‑M1/300‑16 (MQ=255)              aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgc                                                                                                                                                                                                                                                                                        <  1:143156‑M1/300‑16 (MQ=255)              aCTGCCGACATTGTTGCACTACCACCATCCGCTAATGCCGGTTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCGGATGAAACTCAACGCCACTACCTGAACGAACTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAATGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCGCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCGGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTggtgggatatggtgc                                                                                                                                                                                                                                                                                        <  2:199147‑M1/300‑16 (MQ=255)                             gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaag                                                                                                                                                                                                                                                                         <  1:684458‑M1/300‑31 (MQ=255)                             gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaag                                                                                                                                                                                                                                                                         <  2:830138‑M1/300‑31 (MQ=255)                                     aCCATCCGCGAATGCCTGGTTACCTTGATGGTAACGTTCCCAAATGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTTCAACTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgctacgaagaagagcattac                                                                                                                                                                                                                                                                 <  2:49039‑M1/300‑39 (MQ=255)                                       cATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGTTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacga                                                                                                                                                                                                                                                               <  1:186607‑M1/300‑41 (MQ=255)                                       cATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacga                                                                                                                                                                                                                                                               <  2:579378‑M1/300‑41 (MQ=255)                                                         tACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgct                                                                                                                                                                                                                                             <  2:509121‑M1/300‑59 (MQ=255)                                                                 tGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgc                                                                                                                                                                                                                                     <  1:65647‑M1/300‑67 (MQ=255)                                                                           cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgct                                                                                                                                                                                                                           <  2:486584‑M1/300‑77 (MQ=255)                                                                           cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgct                                                                                                                                                                                                                           <  1:886329‑M1/300‑77 (MQ=255)                                                                                        gCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTAGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattaagacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaa                                                                                                                                                                                                              >  2:95283‑M1/1‑211 (MQ=255)                                                                                           tttACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatgg                                                                                                                                                                                                           >  2:377925‑M1/1‑208 (MQ=255)                                                                                                       atgaAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAANCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggt                                                                                                                                                                                               >  2:224907‑M1/1‑196 (MQ=255)                                                                                                            aCTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggctt                                                                                                                                                                                          >  2:716315‑M1/1‑191 (MQ=255)                                                                                                                         aCCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGANCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatntggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttacgctcgctttct                                                                                                                                                                             >  2:725281‑M1/1‑178 (MQ=255)                                                                                                                                  aGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagag                                                                                                                                                                    >  2:577061‑M1/1‑169 (MQ=255)                                                                                                                                     ttGAACTGTATCGTTGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtac                                                                                                                                                                 <  1:915216‑M1/300‑135 (MQ=255)                                                                                                                                     ttGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatttgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtac                                                                                                                                                                 <  1:915207‑M1/300‑135 (MQ=255)                                                                                                                                       gAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactt                                                                                                                                                               <  1:464805‑M1/300‑137 (MQ=255)                                                                                                                                            gTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtg                                                                                                                                                          >  2:519698‑M1/1‑159 (MQ=255)                                                                                                                                              aTCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttacgcttgctttctgctgaagagtactttggtgcc                                                                                                                                                        >  2:321504‑M1/1‑157 (MQ=255)                                                                                                                                                  tGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatttgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcc                                                                                                                                                    <  1:408490‑M1/300‑148 (MQ=255)                                                                                                                                                  tGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcc                                                                                                                                                    <  1:408162‑M1/300‑148 (MQ=255)                                                                                                                                                       tGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttgg                                                                                                                                               >  2:347641‑M1/1‑148 (MQ=255)                                                                                                                                                            gTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagttcgcttcctgctgaagagtacgttggtgccagcatttggcagct                                                                                                                                          >  2:324299‑M1/1‑143 (MQ=255)                                                                                                                                                            gTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagct                                                                                                                                          >  1:263995‑M1/1‑143 (MQ=255)                                                                                                                                                            gTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgatttctgctgaagggtactttggtgccagcctttggcagct                                                                                                                                          >  2:324285‑M1/1‑143 (MQ=255)                                                                                                                                                                   aTCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaag                                                                                                                                   <  1:27735‑M1/300‑165 (MQ=255)                                                                                                                                                                       gccgAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtagtttggtgcaagcctttggcagctctacaagagta                                                                                                                               >  2:256525‑M1/1‑132 (MQ=255)                                                                                                                                                                           aaaGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGNTAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcga                                                                                                                           <  1:673305‑M1/300‑173 (MQ=255)                                                                                                                                                                              ggTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattc                                                                                                                        <  1:4456‑M1/300‑176 (MQ=255)                                                                                                                                                                                      tCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtggtgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaanagtactttggtgccagcctttggcagctctacaagagtatcgattccccataca                                                                                                                >  2:673407‑M1/1‑117 (MQ=255)                                                                                                                                                                                                ggTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctaaaagagtatcgattccccatacaaagcggtact                                                                                                      >  1:521108‑M1/1‑107 (MQ=255)                                                                                                                                                                                                ggTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcctaagctcgctttctgctgaagagtaatttggtgccagcctttggcagctctacaagagtatcgattccacatacaaagcggtact                                                                                                      >  2:779060‑M1/1‑107 (MQ=255)                                                                                                                                                                                                   gtATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGANATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgncagcatttggcagctctacaagagtatcgattccccatacaaagcggtactgaa                                                                                                   >  2:841958‑M1/1‑104 (MQ=255)                                                                                                                                                                                                        cacCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacac                                                                                              <  1:514762‑M1/300‑202 (MQ=255)                                                                                                                                                                                                         acCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtaatttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacact                                                                                             >  2:867479‑M1/1‑98 (MQ=255)                                                                                                                                                                                                          ccATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcacttcctgctgaagagtaatttggtgccagcctttggcagctctaaaagagtatcgcttccccatacaaagcggtactgaaaacaatg                                                                                            >  2:273549‑M1/1‑97 (MQ=255)                                                                                                                                                                                                                       cTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtagtttggtgccagcctttggcagctctacnagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcct                                                                               >  2:693998‑M1/1‑84 (MQ=255)                                                                                                                                                                                                                                aGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctggg                                                                      <  1:272971‑M1/300‑226 (MQ=255)                                                                                                                                                                                                                                 gggCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttcagctcgctttctgctgaagagtactgtggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctggga                                                                     >  2:886393‑M1/1‑74 (MQ=255)                                                                                                                                                                                                                                                        ggATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagtttgctttctgctgaggggtactttggtgttagcatttggcagctctacaagggtatcgatttccaatatggagtggtactgaaaaaactgctgctggaagcctgttcctgggaatgcccgaaaacaagtctgctgg                                              >  2:144607‑M1/1‑51 (MQ=255)                                                                                                                                                                                                                                                           tatCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagttcgctttctgctgaagagtactttggtgccagcatttggcagctctacaagagtatcgattcccaatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcga                                           >  2:189665‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                                           tatCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcga                                           >  2:874276‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                                                      gTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctngaagcctattcctgggaatacccgaacccacgtctgctggcgaaaggtatcaaa                                >  2:688747‑M1/1‑37 (MQ=255)                                                                                                                                                                                                                                                                      gTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaa                                >  1:854579‑M1/1‑37 (MQ=255)                                                                                                                                                                                                                                                                             aTCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagttcgctttctgctgaagagtactttngtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaacagcgtt                         >  2:233273‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                                                             aTCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatccaacagcgtt                         >  1:440241‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                                                             aTCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagcccgctttctgctgaagagtactttggtgccagcctttggcagctctaaaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaacagcgtt                         >  2:432408‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                                                             aTCCTGGCTCGATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagcccgctttctgctgaagagtactttggtgacagcatttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcatgggaatacccgaaccaacgtctgctggcgaaagatatcaaccagcgtt                         >  2:290734‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                                                                       gATAGCGAAGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaacagcgtttgcacgacgg               <  1:770459‑M1/300‑281 (MQ=255)                                                                                                                                                                                                                                                                                              aaGAGCGCCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaacagcgtttgcacgacggcgagatt        <  1:513053‑M1/300‑288 (MQ=255)                                                                                                                                                                                                                                                                                                  gcgcCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaacagcgtttgcacgacggcgagattggat    <  1:234246‑M1/300‑292 (MQ=255)                                                                                                                                                                                                                                                                                                    gcCAATTtgtggaatatggtgccgtgcgacgaagaagagcattacgacgactatgtgatgacgctttacgcgcagggcgtgctgacgccaaatgaatggctggatctcggtggcttaagctcgctttctgctgaagagtactttggtgccagcctttggcagctctacaagagtatcgattccccatacaaagcggtactgaaaacactgctgctggaagcctattcctgggaatacccgaacccacgtctgctggcgaaagatatcaaacagcgtttgcacgacggcgagattgtatcg  <  1:576022‑M1/300‑294 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                      AGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGC  >  NC_000913/3991250‑3991840

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