Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	366,181	Δ93 bp	coding (33‑125/3075 nt)	lacZ ←	beta‑galactosidase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	366181	366273	93	72 [0]	[0] 72	lacZ	beta‑galactosidase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 366180	0 (0.000)	72 (0.990)	67/492	0.2	100%	coding (126/3075 nt)	lacZ	beta‑galactosidase
	?	NC_000913	366274 =	0 (0.000)					coding (32/3075 nt)	lacZ	beta‑galactosidase

AACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCT  >  NC_000913/365936‑366428                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         aacaacCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTGGGTGTAGATGGGCGCATCGTAACCGGGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgac                                                                                                                                                                                                                                                     <  2:298286‑M1/250‑6 (MQ=255)       ccGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGGGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagt                                                                                                                                                                                                                                              <  2:44984‑M1/251‑13 (MQ=255)            ggATTCTCCGTGGGAACAAACGGCGGATTGACCTTAATGGGATAGGTCACGGTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTGTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatc                                                                                                                                                                                                                                         <  2:568518‑M1/251‑18 (MQ=255)                          aaCAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAACGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGGTCGGGCCacgacggccagtgaatccgtaatcatggtc                                                                                                                                                                                                                            >  2:733870‑M1/1‑220 (MQ=255)                           aCAAACGGCGGATTGAACGGAATGGGATAGCTCACGTTGGTGTAGATGGGCGCATCGTATCCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTTTCGGCCTCAGGAAGACCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcat                                                                                                                                                                                                                          <  2:91917‑M1/251‑33 (MQ=255)                                cggcggATTGCACGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCATCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctg                                                                                                                                                                                                                     <  1:140757‑M1/251‑38 (MQ=255)                                 ggcggATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgt                                                                                                                                                                                                                    >  1:272225‑M1/1‑213 (MQ=255)                                   cggATTGACCGTAATGGGATAGGTGACTTTGGTGTAGATGGGCGCATCGTAACGCTGCTTCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCTCCGCTTCTGGTGCCTGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCGACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgttt                                                                                                                                                                                                                  <  1:697691‑M1/251‑41 (MQ=255)                                    ggATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttc                                                                                                                                                                                                                 <  2:513049‑M1/251‑42 (MQ=255)                                    ggATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttc                                                                                                                                                                                                                 >  2:598520‑M1/1‑210 (MQ=255)                                      aTTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcct                                                                                                                                                                                                               >  2:176372‑M1/1‑208 (MQ=255)                                        tGACCGTAAGGGGAGAGGTCACGTTGTTGTAGAGGGGCGCATCGTAACCGTGCATCTGCCAGGTGGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCGGCCAGCTTGCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATCCGCCATTCAGGCTTCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggtcagtgaatccttaatcacggtcatagctgtttcctg                                                                                                                                                                                                              <  2:663771‑M1/250‑45 (MQ=255)                                         gACCGTAATGGTATAGGTCCCGTGGGTGTAGATGGGCGCATCATAACCGTGCATCTGCCATGTTGAGGGGACGCCGACAGTCTCGGCCTCAGGAAGATCTCACTCCAGCCAGCTTTCCGGCACCTCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCTGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgt                                                                                                                                                                                                             <  2:553761‑M1/250‑46 (MQ=255)                                         gACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtg                                                                                                                                                                                                            <  1:377759‑M1/251‑47 (MQ=255)                                               aaTGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaat                                                                                                                                                                                                      <  1:558153‑M1/251‑53 (MQ=255)                                                      aagGTCACGTTGGTGTAGATGGGCGCATCGTAACCGGTCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCTGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTTCCGGAAACCAGGCAAAGCGCCATTCGCCAGTCAGGCTGCGCAACTGTGGGGAAGGGCGATCGGTGCGGTCCacgacggccagttaatccgtaatcatggtcatagctgtttcctgtgtgaaatggttatc                                                                                                                                                                                               <  2:348566‑M1/250‑60 (MQ=255)                                                       aGGTCACGTTGGGGTAGATGGGCGCACCGGAACCGTGCAGCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatcc                                                                                                                                                                                              <  2:635623‑M1/251‑61 (MQ=255)                                                         gTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCCGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgc                                                                                                                                                                                            >  1:288725‑M1/1‑189 (MQ=255)                                                           cACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctc                                                                                                                                                                                          >  2:255291‑M1/1‑187 (MQ=255)                                                                       gATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacac                                                                                                                                                                              <  1:528077‑M1/251‑77 (MQ=255)                                                                            gcgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGAGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcttgtgtgaaattgttatccgctcacaattccacacaacat                                                                                                                                                                         >  2:616446‑M1/1‑170 (MQ=255)                                                                                aTCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgggtcctgtgggaaattgttatccgctcacaattccacacaacatacga                                                                                                                                                                     <  2:320797‑M1/251‑86 (MQ=255)                                                                                  cGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgag                                                                                                                                                                    <  1:176372‑M1/250‑87 (MQ=255)                                                                                       ccGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccgga                                                                                                                                                               >  1:669586‑M1/1‑159 (MQ=255)                                                                                        cGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaag                                                                                                                                                             <  1:152487‑M1/251‑94 (MQ=255)                                                                                              tCTGCCAGTTTGAGGGGACGACGACAGTATCGGACTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCCCCGCTTCTGGAGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGCCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataa                                                                                                                                                        >  2:27944‑M1/1‑152 (MQ=255)                                                                                                  ccAGTTGGAGGGGACGACGACAGTAGCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgt                                                                                                                                                   <  2:527629‑M1/251‑104 (MQ=255)                                                                                                     gTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaa                                                                                                                                                 >  1:565910‑M1/1‑145 (MQ=255)                                                                                                      tttGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaa                                                                                                                                                 >  1:620252‑M1/1‑144 (MQ=255)                                                                                                         gAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGCTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcct                                                                                                                                            >  1:260059‑M1/1‑141 (MQ=255)                                                                                                         gAGGGGACGACGACAGGATCGCCCTCAGGAAGATCGCAACCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcc                                                                                                                                             <  2:214702‑M1/250‑110 (MQ=255)                                                                                                           ggggACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctgg                                                                                                                                          <  1:733870‑M1/251‑113 (MQ=255)                                                                                                                    gacAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccg                                                                                                                                                                 <  1:132944‑M1/219‑90 (MQ=255)                                                                                                                    gacAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccg                                                                                                                                                                 >  2:132944‑M1/1‑130 (MQ=255)                                                                                                                     acAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCTGTGCGGGCCacgacggccagtgaatccgtaagcatggtcatagctgtttcctgtgtgaaattgttatccgctcccaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaat                                                                                                                                <  2:266434‑M1/251‑123 (MQ=255)                                                                                                                      cAGTATCGGCCTCAGGCAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacaccacatacgagccggaagcataaagtgtaaagcctggggtgcctaatg                                                                                                                               >  1:607891‑M1/1‑128 (MQ=255)                                                                                                                      cAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatt                                                                                                                               >  2:559845‑M1/1‑128 (MQ=255)                                                                                                                       aGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatg                                                                                                                               >  2:245970‑M1/1‑127 (MQ=255)                                                                                                                          aTCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCaggacggccagtgaatccgtattcatggtcatagcttgttcctgtgtgaaattgtgatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtg                                                                                                                           <  1:298208‑M1/251‑128 (MQ=255)                                                                                                                             ggCCTCAGGAAGATCGCACCCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagc                                                                                                                        <  2:369633‑M1/251‑131 (MQ=255)                                                                                                                              gCCTCAGGAAGATCGCGCTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGGGCGGGCCaggacggccagtgaatccgtaatcatggtcatagctgtttccggtgtgaaatgtttagccgctcgcaattccacacaacaagcgagcctgaagcataaagtggaatgcctgggttgcctaatgagtgagct                                                                                                                       <  2:685429‑M1/251‑132 (MQ=255)                                                                                                                                  cAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaact                                                                                                                   >  2:28381‑M1/1‑116 (MQ=255)                                                                                                                                           cGCACTCCAGCCATCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCAGTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaat                                                                                                          <  2:565910‑M1/251‑145 (MQ=255)                                                                                                                                                  cagccagcTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttg                                                                                                   <  1:44652‑M1/251‑152 (MQ=255)                                                                                                                                                   agccagcTTTCCTGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGTCGAGCTGTGCGGGCCacgacggccagtgaatccgtaatcattgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgc                                                                                                  <  1:489165‑M1/251‑153 (MQ=255)                                                                                                                                                     ccagcATTCCGGCACCGCTTCTGGTGCCGGAAACCATGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGCGAGAGACacgacggccagtgaatccgtaatcatggtcataactgtttcctgtgtgaaattgttatccgctcacaattccacacaacatccgcgccggaagcttaaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgc                                                                                                >  2:325652‑M1/1‑97 (MQ=255)                                                                                                                                                      cagcTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaag                                                                                                                                                             >  1:710337‑M1/1‑96 (MQ=255)                                                                                                                                                      cagcTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaag                                                                                                                                                             <  2:710337‑M1/189‑94 (MQ=255)                                                                                                                                                                 cACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgct                                                                                    >  1:153964‑M1/1‑85 (MQ=255)                                                                                                                                                                   ccGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgtgttgcgctcactgcccgcttt                                                                                  >  2:564812‑M1/1‑83 (MQ=255)                                                                                                                                                                            gTGCCGGAAACCAGGCAAAGCGCCATTCGCCCTTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggacagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcggg                                                                         >  2:710299‑M1/1‑74 (MQ=255)                                                                                                                                                                            gTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcggg                                                                         >  2:103565‑M1/1‑74 (MQ=255)                                                                                                                                                                                 ggAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacc                                                                    >  1:198826‑M1/1‑69 (MQ=255)                                                                                                                                                                                   aaaCCAGGCAAAGCGCCATTCGACATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctg                                                                  >  1:274629‑M1/1‑67 (MQ=255)                                                                                                                                                                                       cAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgt                                                              <  2:669586‑M1/251‑189 (MQ=255)                                                                                                                                                                                                     ttcgccattcAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcatta                                                 >  2:407392‑M1/1‑49 (MQ=255)                                                                                                                                                                                                          cattcAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaat                                           <  1:103565‑M1/251‑208 (MQ=255)                                                                                                                                                                                                            ttcAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgc                                                                                                       <  1:691693‑M1/189‑148 (MQ=255)                                                                                                                                                                                                            ttcAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagagaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgc                                                                                                       >  2:691693‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                     ggcaATTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgatctaactcacattaattacgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcg                                <  2:607891‑M1/250‑219 (MQ=255)                                                                                                                                                                                                                               gggAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagagg                      >  2:574366‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                   aGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggg                  >  2:224533‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                        gATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaatcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccatcgcgcggtgagaggcggtttgcg             <  2:638390‑M1/251‑238 (MQ=255)                                                                                                                                                                                                                                        gATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtcaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggtttg               >  1:668523‑M1/1‑14 (MQ=255)                                                                                                                                                                                                                                               gCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgctagctgcattaatgaatcggccaacgcgcggggagaggcggtttgcgtattggg      <  1:27944‑M1/251‑245 (MQ=255)                                                                                                                                                                                                                                               gCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggttttcgtctttgg      >  1:100373‑M1/1‑7 (MQ=255)                                                                                                                                                                                                                                                cGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcagaaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaagtgcgttgcgctcactgctcgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggtttgcgtattgggc     <  1:363808‑M1/251‑246 (MQ=255)                                                                                                                                                                                                                                                 gggCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggtttgcgtattgggcg    <  1:113832‑M1/251‑247 (MQ=255)                                                                                                                                                                                                                                                  ggCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggttttcggttttggggg   >  1:707293‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                  ggCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggttttcgggatcggg     >  1:31853‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                  ggCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgggaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggtttgcg             <  2:31853‑M1/241‑238 (MQ=255)                                                                                                                                                                                                                                                   gCCacgacggccagtgaatccgtaatcatggtcatagctgttccctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcggggagaggcggtttgcgtattgggcgcc  <  1:564384‑M1/251‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         AACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCT  >  NC_000913/365936‑366428

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