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

GAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGG  >  NC_000913/2514028‑2514520                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          gAAAAGTGCTTCCACGCTCTCGCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCT‑ACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGCAGGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacac                                                                                                                                                                                                                                                   <  1:207951‑M1/251‑9 (MQ=255)  aaaaTTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGCAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacac                                                                                                                                                                                                                                                   >  2:495075‑M1/1‑243 (MQ=255)   aaaTTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAACAAAGCCGGGGATAAgagtacact                                                                                                                                                                                                                                                  >  1:448421‑M1/1‑242 (MQ=255)     aTTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactac                                                                                                                                                                                                                                                >  1:178730‑M1/1‑240 (MQ=255)       tGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactc                                                                                                                                                                                                                                              >  1:409610‑M1/1‑238 (MQ=255)            ccACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAAGCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcac                                                                                                                                                                                                                                          <  1:470308‑M1/250‑18 (MQ=255)                gCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAATGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatc                                                                                                                                                                                                                                     <  1:312674‑M1/251‑23 (MQ=255)                gCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatc                                                                                                                                                                                                                                     <  1:98801‑M1/251‑23 (MQ=255)                           gTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggagttttc                                                                                                                                                                                                                            >  2:395315‑M1/1‑218 (MQ=255)                               tGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcac                                                                                                                                                                                                                      >  2:296987‑M1/1‑214 (MQ=255)                                        catcTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGCTTATCGCAGGTGCAGTTAAAGGCCTGACTGACGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGGCCATCGCCGCAATGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataat                                                                                                                                                                                                             >  2:428855‑M1/1‑205 (MQ=255)                                               gtgtTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagttt                                                                                                                                                                                                       <  1:488324‑M1/250‑53 (MQ=255)                                                 gtTCCTGGTTTCCTTCGCTAACTTCTCTTCACTCGGGATTATCGCCGGTGCAGATAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGCGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATCAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttatt                                                                                                                                                                                                    >  1:417882‑M1/1‑196 (MQ=255)                                                 gtTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTTGTGAGTGTGCTGTCTGCGGCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactcctctcactgatcggggttttcaggcacacaaataataagtttatt                                                                                                                                                                                                    <  2:110920‑M1/251‑56 (MQ=255)                                                         tttCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgc                                                                                                                                                                                            >  1:314378‑M1/1‑188 (MQ=255)                                                              ttCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagc                                                                                                                                                                                       >  2:423514‑M1/1‑183 (MQ=255)                                                               tCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaag                                                                                                                                                                                        <  1:101127‑M1/249‑68 (MQ=255)                                                                 gCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaa                                                                                                                                                                                    <  2:475550‑M1/251‑72 (MQ=255)                                                                           ctTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggcc                                                                                                                                                                           >  1:457125‑M1/1‑170 (MQ=255)                                                                             tCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggtt                                                                                                                                                                                                                               <  1:93147‑M1/196‑29 (MQ=255)                                                                             tCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggtt                                                                                                                                                                                                                               >  2:93147‑M1/1‑168 (MQ=255)                                                                               aaTCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgacgcaaacagggggccagcct                                                                                                                                                                      >  1:139614‑M1/1‑166 (MQ=255)                                                                               aaTCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGGTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaatgataagtgtattcgcatcgcgaagcaaacaggtggccagcct                                                                                                                                                                      <  2:453966‑M1/251‑86 (MQ=255)                                                                                  cGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgtt                                                                                                                                                                   >  1:181744‑M1/1‑163 (MQ=255)                                                                                     gATTAGCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGGGGTTTCGCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAATAAAGCCGGGGATAAgagtacacaactctcactgatcggggttttcaggcacacaactaataagtttattggcatcgcgaagaaaacagggggccagcctgttgc                                                                                                                                                                 >  2:81104‑M1/1‑160 (MQ=255)                                                                                      aTTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccg                                                                                                                                                               <  1:290840‑M1/251‑93 (MQ=255)                                                                                        tATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgat                                                                                                                                                             >  1:46195‑M1/1‑157 (MQ=255)                                                                                         aTCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatt                                                                                                                                                            >  2:463978‑M1/1‑156 (MQ=255)                                                                                           cGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTGCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataaggttattcgcatcgcgaagcaaacagggggccagcctgttgccgattt                                                                                                                                                           <  1:50625‑M1/250‑97 (MQ=255)                                                                                             cAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgcgtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttc                                                                                                                                                        >  2:19336‑M1/1‑152 (MQ=255)                                                                                               ggTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcat                                                                                                                                                      >  1:342702‑M1/1‑150 (MQ=255)                                                                                                gTGCAGTTAAAGGCCTGAATTAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatc                                                                                                                                                     <  1:313335‑M1/251‑103 (MQ=255)                                                                                                  gCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatcca                                                                                                                                                   <  2:317656‑M1/251‑105 (MQ=255)                                                                                                  gCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggtt                                                                                                                                                                                                                               <  1:388075‑M1/175‑29 (MQ=255)                                                                                                  gCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggtt                                                                                                                                                                                                                               >  2:388075‑M1/1‑147 (MQ=255)                                                                                                     gTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatg                                                                                                                                                >  2:38804‑M1/1‑144 (MQ=255)                                                                                                                gaatgaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcacc                                                                                                                                     <  2:153395‑M1/251‑119 (MQ=255)                                                                                                                  atgaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcacc                                                                                                                                     >  2:394752‑M1/1‑131 (MQ=255)                                                                                                                     aaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcg                                                                                                                                <  1:295904‑M1/251‑124 (MQ=255)                                                                                                                      agagCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcgg                                                                                                                               >  1:129275‑M1/1‑127 (MQ=255)                                                                                                                                         tctcGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagac                                                                                                            >  1:342617‑M1/1‑108 (MQ=255)                                                                                                                                            cGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgt                                                                                                         >  1:200656‑M1/1‑105 (MQ=255)                                                                                                                                              cTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtat                                                                                                       >  2:382696‑M1/1‑103 (MQ=255)                                                                                                                                                tCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttc                                                                                                                                        <  1:223127‑M1/216‑116 (MQ=255)                                                                                                                                                tCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttc                                                                                                                                        >  2:223127‑M1/1‑101 (MQ=255)                                                                                                                                                  ggTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagcctgaattac                                                                                                    >  1:302671‑M1/1‑99 (MQ=255)                                                                                                                                                            tGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggca                                                                                          >  2:21142‑M1/1‑89 (MQ=255)                                                                                                                                                             ggTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaa                                                                                        <  1:38804‑M1/251‑164 (MQ=255)                                                                                                                                                                     gCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagag                                                                                <  2:433298‑M1/251‑172 (MQ=255)                                                                                                                                                                         tACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagaca                                                                             >  1:65290‑M1/1‑76 (MQ=255)                                                                                                                                                                               ggTGAGTGTGCTGTCTGCGTAAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgataggggttatcaggcacacaaataataagtttattcgcatcgcgaagcaaccagggggccagcctgttgccgctttttcatccaatgatgctttcaccaggcggcactaacccggatgtagactggattacaaaagcggcaaaaagcagagacaacaaacc                                                                      >  2:28664‑M1/1‑70 (MQ=255)                                                                                                                                                                                   aGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccg                                                                  >  1:415987‑M1/1‑66 (MQ=255)                                                                                                                                                                                                       cGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttgga                                               >  1:118636‑M1/1‑46 (MQ=255)                                                                                                                                                                                                         cagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctccctgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacaggcggccagcctgttgccgattttgcatccaatgatgctttcaccaggcggcactaaccctgacttagactgtattacaaaatcggcaaaaagcagagacaaaaatcccccgcttgtcagctagtgtttgaaat                                            <  2:417882‑M1/251‑208 (MQ=255)                                                                                                                                                                                                                 gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtgga                                    <  2:314378‑M1/251‑216 (MQ=255)                                                                                                                                                                                                                 gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtgga                                    <  1:51538‑M1/251‑216 (MQ=255)                                                                                                                                                                                                                  tGCTGTAAGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtgg                                     >  1:84948‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                        aaGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacaggggtccagcctgttgccgatttttcatccaatgatgctttc                                                                                                                                        <  1:385883‑M1/144‑116 (MQ=255)                                                                                                                                                                                                                        aaGACCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttc                                                                                                                                        >  2:385883‑M1/1‑29 (MQ=255)                                                                                                                                                                                                                          gACCATTCAGAAAAAAGCCGGGGATAAgagtacactcctctacctgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagccttttgccgatttttcagccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagggagggttggaaatttggtggagctaagcgg                           <  2:46009‑M1/251‑225 (MQ=255)                                                                                                                                                                                                                           aCCATACATAAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcggg                          <  1:136039‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                              atacataAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatc                       <  2:448421‑M1/251‑229 (MQ=255)                                                                                                                                                                                                                                acataAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgaggnttggaaatttggtggagctaagcgggatcg                      >  2:225475‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                acataAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcg                      >  2:94443‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                 cataAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcga                     >  1:429074‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                  ataAAAAAGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaac                   >  1:92172‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                     aaaaaGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaaccgc                >  1:348079‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                                      aaaaGCCGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaaccgct               <  2:269895‑M1/251‑237 (MQ=255)                                                                                                                                                                                                                                            cGGGGATAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaaccgctgacctc         <  2:65290‑M1/251‑243 (MQ=255)                                                                                                                                                                                                                                            cGGGGATAAgagtacactactctcactgagcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaaccgctgacctc         <  2:409610‑M1/251‑243 (MQ=255)                                                                                                                                                                                                                                                 aTAAgagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaaccgctgacctcttgc     >  1:72488‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                   aagagtacactactctcactgatcggggttttcaggcacacaaataataagtttattcgcatcgcgaagcaaacagggggccagcctgttgccgatttttcatccaatgatgctttcaccaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaagcagagacaaaaaacccccgctttgcagcgagggttggaaatttggtggagctaagcgggatcgaaccgctgacctcttgcatg  >  2:9872‑M1/1‑2 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          GAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGG  >  NC_000913/2514028‑2514520

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: ‑