Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,514,272	Δ3,561 bp		insL3–[pdeA]	insL3, [pdeA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2514272	2517832	3561	72 [0]	[0] 72	insL3–[pdeA]	insL3,[pdeA]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2514271	0 (0.000)	72 (1.100)	65/478	0.0	100%	intergenic (+27/‑60)	nupC/insL3	nucleoside:H(+) symporter NupC/putative IS186/IS421 transposase
	?	NC_000913	2517833 =	0 (0.000)					intergenic (‑1/+208)	pdeA/alaX	putative c‑di‑GMP phosphodiesterase PdeA/tRNA‑Ala

GAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTGACAGCATTACGCACAGTGGCATCCACGCCAGAAAGATATAACTGGAATGGATAATGGCAAGTGGCGAAATAAAACGGGATAGCTGAATAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTATCAGGTTATGCTCCACAAACATGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATG  >  NC_000913/2517590‑2518081                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         gaaaagtgcttccacgctctcgccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgcaggaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTacac                                                                                                                                                                                                                                                   <  1:207951‑M2/8‑1 (MQ=255) aaaattgcttccacgctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgcagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTacac                                                                                                                                                                                                                                                   >  2:495075‑M2/244‑251 (MQ=255)  aaattgcttccacgctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataacaaagccggggataaGAGTACact                                                                                                                                                                                                                                                  >  1:448421‑M2/243‑251 (MQ=255)    attgcttccacgctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACactac                                                                                                                                                                                                                                                >  1:178730‑M2/241‑251 (MQ=255)      tgcttccacgctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTActc                                                                                                                                                                                                                                              >  1:409610‑M2/239‑251 (MQ=255)           ccacgctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaagcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCAc                                                                                                                                                                                                                                          <  1:470308‑M2/17‑1 (MQ=255)               gctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaatgtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATc                                                                                                                                                                                                                                     <  1:312674‑M2/22‑1 (MQ=255)               gctctctccgcgtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATc                                                                                                                                                                                                                                     <  1:98801‑M2/22‑1 (MQ=255)                          gtgctgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGAGTTTTc                                                                                                                                                                                                                            >  2:395315‑M2/219‑249 (MQ=255)                              tgaaggcatcatctctgtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGcac                                                                                                                                                                                                                      >  2:296987‑M2/215‑251 (MQ=255)                                       catctctgtgttcctggtttccttcgctaacttctcttcaatcgggcttatcgcaggtgcagttaaaggcctgactgacgagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcggccatcgccgcaatggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAaataat                                                                                                                                                                                                             >  2:428855‑M2/206‑251 (MQ=255)                                              gtgttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGttt                                                                                                                                                                                                       <  1:488324‑M2/52‑1 (MQ=255)                                                gttcctggtttccttcgctaacttctcttcactcgggattatcgccggtgcagataaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagcgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacatcaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTAtt                                                                                                                                                                                                    >  1:417882‑M2/197‑251 (MQ=255)                                                gttcctggtttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctacccttgtgagtgtgctgtctgcggcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTCCTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTAtt                                                                                                                                                                                                    <  2:110920‑M2/55‑1 (MQ=255)                                                        tttccttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATcgc                                                                                                                                                                                            >  1:314378‑M2/189‑251 (MQ=255)                                                             ttcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGc                                                                                                                                                                                       >  2:423514‑M2/184‑251 (MQ=255)                                                              tcgctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAg                                                                                                                                                                                        <  1:101127‑M2/67‑1 (MQ=255)                                                                gctaacttctcttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCaaa                                                                                                                                                                                    <  2:475550‑M2/71‑1 (MQ=255)                                                                          cttcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGcc                                                                                                                                                                           >  1:457125‑M2/171‑250 (MQ=255)                                                                            tcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGtt                                                                                                                                                                                                                               >  2:93147‑M2/169‑196 (MQ=255)                                                                            tcaatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGtt                                                                                                                                                                                                                               <  1:93147‑M2/28‑1 (MQ=255)                                                                              aatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGACGCAAACAGGGGGCCAGCCt                                                                                                                                                                      >  1:139614‑M2/167‑251 (MQ=255)                                                                              aatcgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtgggttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATGATAAGTGTATTCGCATCGCGAAGCAAACAGGTGGCCAGCCt                                                                                                                                                                      <  2:453966‑M2/85‑1 (MQ=255)                                                                                 cgggattatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGtt                                                                                                                                                                   >  1:181744‑M2/164‑251 (MQ=255)                                                                                    gattagcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacggggtttcgcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataataaagccggggataaGAGTACACAACTCTCACTGATCGGGGTTTTCAGGCACACAACTAATAAGTTTATTGGCATCGCGAAGAAAACAGGGGGCCAGCCTGTTGc                                                                                                                                                                 >  2:81104‑M2/161‑250 (MQ=255)                                                                                     attatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCg                                                                                                                                                               <  1:290840‑M2/92‑1 (MQ=255)                                                                                       tatcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGAt                                                                                                                                                             >  1:46195‑M2/158‑251 (MQ=255)                                                                                        atcgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGAtt                                                                                                                                                            >  2:463978‑M2/157‑251 (MQ=255)                                                                                          cgcaggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgctgcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGGTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGAttt                                                                                                                                                           <  1:50625‑M2/96‑1 (MQ=255)                                                                                            caggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGCGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTc                                                                                                                                                        >  2:19336‑M2/153‑251 (MQ=255)                                                                                              ggtgcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCAt                                                                                                                                                      >  1:342702‑M2/151‑251 (MQ=255)                                                                                               gtgcagttaaaggcctgaattaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATc                                                                                                                                                     <  1:313335‑M2/102‑1 (MQ=255)                                                                                                 gcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGtt                                                                                                                                                                                                                               >  2:388075‑M2/148‑175 (MQ=255)                                                                                                 gcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGtt                                                                                                                                                                                                                               <  1:388075‑M2/28‑1 (MQ=255)                                                                                                 gcagttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCa                                                                                                                                                   <  2:317656‑M2/104‑1 (MQ=255)                                                                                                    gttaaaggcctgaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAatg                                                                                                                                                >  2:38804‑M2/145‑251 (MQ=255)                                                                                                               gaatgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCAcc                                                                                                                                     <  2:153395‑M2/118‑1 (MQ=255)                                                                                                                 atgaagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCAcc                                                                                                                                     >  2:394752‑M2/132‑249 (MQ=255)                                                                                                                    aagagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAggcg                                                                                                                                <  1:295904‑M2/123‑1 (MQ=255)                                                                                                                     agagcaaggtaacgtggtttctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAggcgg                                                                                                                               >  1:129275‑M2/128‑251 (MQ=255)                                                                                                                                        tctcgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGAc                                                                                                            >  1:342617‑M2/109‑251 (MQ=255)                                                                                                                                           cgcttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGt                                                                                                         >  1:200656‑M2/106‑251 (MQ=255)                                                                                                                                             cttcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTAt                                                                                                       >  2:382696‑M2/104‑251 (MQ=255)                                                                                                                                               tcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTc                                                                                                                                        <  1:223127‑M2/115‑1 (MQ=255)                                                                                                                                               tcggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTc                                                                                                                                        >  2:223127‑M2/102‑216 (MQ=255)                                                                                                                                                 ggtctgaagctggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGCCTGAATTAc                                                                                                    >  1:302671‑M2/100‑250 (MQ=255)                                                                                                                                                           tggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCa                                                                                          >  2:21142‑M2/90‑250 (MQ=255)                                                                                                                                                            ggtttacggctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaa                                                                                        <  1:38804‑M2/163‑1 (MQ=255)                                                                                                                                                                    gctctaccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCagag                                                                                <  2:433298‑M2/171‑1 (MQ=255)                                                                                                                                                                        taccctggtgagtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACa                                                                             >  1:65290‑M2/77‑250 (MQ=255)                                                                                                                                                                              ggtgagtgtgctgtctgcgtaaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATAGGGGTTATCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAACCAGGGGGCCAGCCTGTTGCCGCTTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGATGTAGACTGGATTACAAAAGCGGCAAAAAGCAGAGACAACAAAcc                                                                      >  2:28664‑M2/71‑251 (MQ=255)                                                                                                                                                                                  agtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCg                                                                  >  1:415987‑M2/67‑251 (MQ=255)                                                                                                                                                                                                      cgcagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGa                                               >  1:118636‑M2/47‑250 (MQ=255)                                                                                                                                                                                                        cagcactggtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCCCTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGCGGCCAGCCTGTTGCCGATTTTGCATCCAATGATGCTTTCACCAGGCGGCACTAACCCTGACTTAGACTGTATTACAAAATCGGCAAAAAGCagagacaaaaatcccccgcttgtcagctagtgtttgaaat                                            <  2:417882‑M2/207‑20 (MQ=255)                                                                                                                                                                                                                gtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGa                                    <  2:314378‑M2/215‑1 (MQ=255)                                                                                                                                                                                                                gtgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGa                                    <  1:51538‑M2/215‑1 (MQ=255)                                                                                                                                                                                                                 tgctgtaagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTtggtgg                                     >  1:84948‑M2/36‑249 (MQ=255)                                                                                                                                                                                                                       aagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGTCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTc                                                                                                                                        <  1:385883‑M2/115‑1 (MQ=255)                                                                                                                                                                                                                       aagaccatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTc                                                                                                                                        >  2:385883‑M2/30‑144 (MQ=255)                                                                                                                                                                                                                         gaccattcagaaaaaagccggggataaGAGTACACTCCTCTACCTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTTTTGCCGATTTTTCAGCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGGGAGGGTTGGAAATTTGGTGGAGCTAAGCgg                           <  2:46009‑M2/224‑1 (MQ=255)                                                                                                                                                                                                                          accatacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCggg                          <  1:136039‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                             atacataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATc                       <  2:448421‑M2/228‑1 (MQ=255)                                                                                                                                                                                                                               acataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGNTTGGAAATTTGGTGGAGCTAAGCGGGATCg                      >  2:225475‑M2/22‑250 (MQ=255)                                                                                                                                                                                                                               acataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCg                      >  2:94443‑M2/22‑250 (MQ=255)                                                                                                                                                                                                                                cataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGa                     >  1:429074‑M2/21‑250 (MQ=255)                                                                                                                                                                                                                                 ataaaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAAc                   >  1:92172‑M2/20‑251 (MQ=255)                                                                                                                                                                                                                                    aaaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGc                >  1:348079‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                                     aaaagccggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCt               <  2:269895‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                                           cggggataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctc         <  2:65290‑M2/242‑1 (MQ=255)                                                                                                                                                                                                                                           cggggataaGAGTACACTACTCTCACTGAGCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctc         <  2:409610‑M2/242‑1 (MQ=255)                                                                                                                                                                                                                                                ataaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGc     >  1:72488‑M2/5‑250 (MQ=255)                                                                                                                                                                                                                                                  aaGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATg  >  2:9872‑M2/3‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         GAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGATTGACAGCATTACGCACAGTGGCATCCACGCCAGAAAGATATAACTGGAATGGATAATGGCAAGTGGCGAAATAAAACGGGATAGCTGAATAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTATCAGGTTATGCTCCACAAACATGAGTACACTACTCTCACTGATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATG  >  NC_000913/2517590‑2518081

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