Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,588,756	Δ970 bp		[ydeS]–[ydeT]	[ydeS], [ydeT]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1588756	1589725	970	64 [0]	[0] 63	[ydeS]–[ydeT]	[ydeS],[ydeT]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1588755	0 (0.000)	63 (0.780)	56/490	0.4	100%	coding (85/531 nt)	ydeS	putative fimbrial protein YdeS
	?	NC_000913	1589726 =	0 (0.000)					coding (276/1149 nt)	ydeT	fimbrial usher domain‑containing protein YdeT

AAGAGAACTGTAACCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATGCCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTAAGCGTAAAAGCGAGTAAATGATCCCGATCGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCT  >  NC_000913/1589477‑1589958                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         ggaggccgtattgcttccttcatccagtttcaacaacgtattaatttctgcatcctcaataccgttaaatgcaacccgaaaccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggatttctgctacaggtcgttggacattgtatggcgctgtttttttgaagatctacgggaaaattaagcgactccgatgagactgtgaagccataagcgaggccgcgaccgctcatGt                                                                                                                                                                                                                                         >  2:92394‑M2/250‑251 (MQ=255)    ggccgtattgcttccctcatccagtttcaacaaagtattattttcttcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGtttt                                                                                                                                                                                                                                      <  2:432731‑M2/5‑1 (MQ=255)      ccgtattgcttccttcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGc                                                                                                                                                                                                                                    <  1:55507‑M2/7‑1 (MQ=255)         tattgcttccttcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCat                                                                                                                                                                                                                                 >  1:167621‑M2/242‑251 (MQ=255)               ttccttcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaccccgaaccacagttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcaaatccaatgagtctgtgccgccataatcgaggacgcgcccgctaatTTTTTACCATAtatt                                                                                                                                                                                                                            >  2:265968‑M2/236‑250 (MQ=255)               ttccttcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATatattt                                                                                                                                                                                                                           >  1:172791‑M2/236‑250 (MQ=255)                    tcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTgg                                                                                                                                                                                                                       <  1:488173‑M2/20‑1 (MQ=255)                      atccagtttcaacaaagtattattttctgcatcctcaattccgttaaatgcaccccgaacccccgttgtccctttgctgcattcacttaacgttatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaatattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAAt                                                                                                                                                                                                                    <  1:125118‑M2/23‑1 (MQ=255)                               caaccaagtattattttctgcatcctcaataccgatcaatgcaacccgcacccccgttgtccctttgctgcattcccttaacgtaatctgaaaagggacgactggacttgtgctaccggtcgttggcaattatctggcactgtttttttggagatctaaggtaaaattaagcgaatccgatgagactgtgcagccataatcgagcacgcgccagcgaagGTTTTGCCAtatattattggactcactatcat                                                                                                                                                                                                           >  2:838805‑M2/220‑240 (MQ=255)                               caaaacagtattattttctgcatcctcaataccgttaactgcaaccagaaccccagttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattcagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctcatGTTTTGCCATAAAatattggaatgacggatat                                                                                                                                                                                                           >  2:570929‑M2/220‑242 (MQ=255)                                aactaagtatttttttttgcatcattaatacagttaaatgcaacccgacaccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgcgaacggtcgttggaaattgtctggcactgtttttttggagatcttcggtaaaattaagcgattccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAAtt                                                                                                                                                                                                          <  1:546039‑M2/33‑1 (MQ=255)                                aacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAAtt                                                                                                                                                                                                          >  1:262928‑M2/219‑251 (MQ=255)                                       tattattttctgcatcctcaataccgttaaatgcaacccgaaccccctttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGa                                                                                                                                                                                                    >  2:241857‑M2/212‑250 (MQ=255)                                       tattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgctgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGa                                                                                                                                                                                                    >  1:120780‑M2/212‑250 (MQ=255)                                          tattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatccacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTc                                                                                                                                                                                                >  1:11774‑M2/209‑251 (MQ=255)                                                   catcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCAATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTaatata                                                                                                                                                                                        >  1:110374‑M2/200‑249 (MQ=255)                                                    atcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCAc                                                                                                                                                                                      >  2:289003‑M2/199‑251 (MQ=255)                                                         caataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAAc                                                                                                                                                                                 >  1:260422‑M2/194‑251 (MQ=255)                                                               cgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatacgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGACGTCGAAATATAACCAaacggaac                                                                                                                                                                            >  2:619390‑M2/188‑247 (MQ=255)                                                                     atgtaacccgaacccccgtggtccctttgtggctgtctcgtaacgttatttgaaaagggacgggtggccttgttctaccggtctttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccttaatcgaggacgcgcccgttaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGcacattaacc                                                                                                                                                                     <  2:63595‑M2/70‑5 (MQ=255)                                                                     atgcaacccgaaaccacgctgtcccgttgctgcattctcttaacgtaagctgaaaagttacggatggacgtgtgctaccggtcgttggaaatggtctggcactgtttttttggagatctacgggaaaattacgcgaatccgatgagactttgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAAt                                                                                                                                                                     <  2:692026‑M2/70‑1 (MQ=255)                                                                                         tgccctttgcttcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaagtgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGt                                                                                                                                                  <  2:458086‑M2/89‑1 (MQ=255)                                                                                              tttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGc                                                                                                                                            >  1:172397‑M2/157‑251 (MQ=255)                                                                                              tttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGc                                                                                                                                            >  1:249578‑M2/157‑251 (MQ=255)                                                                                                  ctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgcggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTACTATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGtgc                                                                                                                                        >  2:637796‑M2/153‑250 (MQ=255)                                                                                                    gcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCgctgcgaggcg                                                                                                                                      >  2:184234‑M2/151‑246 (MQ=255)                                                                                                         cacttaacgaaatctgaaaagggccggctggacttgtgctaccggtcgttggaaattgtctggaactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGACTAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTc                                                                                                                                 >  2:508353‑M2/146‑251 (MQ=255)                                                                                                              aacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGt                                                                                                                             >  1:43275‑M2/141‑250 (MQ=255)                                                                                                                 gtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACt                                                                                                                         >  1:857907‑M2/138‑251 (MQ=255)                                                                                                                       tgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggaagcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAAATTTGGCg                                                                                                                   >  2:215089‑M2/132‑251 (MQ=255)                                                                                                                                 cggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTg                                                                                                         <  1:200881‑M2/130‑1 (MQ=255)                                                                                                                                 cggctggacttgtgctaccggtcgctggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgaagagactgtgcagccataatcgaggaagagcccgctaatGGTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTACCTTTGGCGACTGGCCCTg                                                                                                         >  2:517811‑M2/122‑251 (MQ=255)                                                                                                                                              gctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTAc                                                                                            <  1:172039‑M2/143‑1 (MQ=255)                                                                                                                                                        gttggaaattgtatggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgccagctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTAAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGCTATTTGCTGGTCGCTGCGTGACGTGTTCCAGTACCTTTGGAGACTGGCACTGAAAAATGTCGTACAGTAATTTcc                                                                                  >  2:160732‑M2/99‑251 (MQ=255)                                                                                                                                                        gttggaaattgacgtgcactgttttgttggagatctacggtaaacttaagcgaatccgttgagactctgcagccataatcgaggacgcgccctcgaatGTTTTGTCAAATATGAGTGGAATAACTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACAGTGAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTGTGGCGACTGGCACTGCAAAATGTCGTACCGTAATTtca                                                                                  <  2:81463‑M2/153‑2 (MQ=255)                                                                                                                                                             aaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAgct                                                                             <  1:413233‑M2/158‑1 (MQ=255)                                                                                                                                                               attgtcgggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAgctgc                                                                           <  1:272862‑M2/160‑1 (MQ=255)                                                                                                                                                                 tgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGctgctg                                                                         >  2:168852‑M2/90‑251 (MQ=255)                                                                                                                                                                 ggtttggcactgtttttgtggagaactacggtaaaagtaagcgaatccgataagactgtgcatccataatcgaggacgcgcccgctaatGTTTTTCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAgctgct                                                                          <  2:249578‑M2/161‑1 (MQ=255)                                                                                                                                                                  gtctggcactggtttttgggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGctgctg                                                                         <  2:167621‑M2/162‑1 (MQ=255)                                                                                                                                                                            gtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCtt                                                               <  1:363257‑M2/172‑1 (MQ=255)                                                                                                                                                                            gtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCt                                                                <  1:276827‑M2/171‑1 (MQ=255)                                                                                                                                                                                 tttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTGCCAAGCTTCTGAGAGATGCTTATTTg                                                          <  1:275101‑M2/177‑1 (MQ=255)                                                                                                                                                                                     gagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgtggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTAGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCACGCTGCTGAGAGATGCTTATTTGTTCCt                                                     >  1:376012‑M2/70‑251 (MQ=255)                                                                                                                                                                                       gatttacgttaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGa                                                   <  2:260422‑M2/184‑1 (MQ=255)                                                                                                                                                                                         tctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGCGATGCTTATTTGTTCCTGAcc                                                 >  2:118728‑M2/66‑251 (MQ=255)                                                                                                                                                                                             cggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGc                                             >  1:61203‑M2/62‑251 (MQ=255)                                                                                                                                                                                                  aaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACt                                        <  1:637796‑M2/195‑1 (MQ=255)                                                                                                                                                                                                      taagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTa                                     <  1:495892‑M2/198‑1 (MQ=255)                                                                                                                                                                                                          cgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTAcc                                                                                           >  2:291726‑M2/49‑192 (MQ=255)                                                                                                                                                                                                          cgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTAcc                                                                                           <  1:291726‑M2/144‑1 (MQ=255)                                                                                                                                                                                                           gaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAg                               <  2:61203‑M2/204‑1 (MQ=255)                                                                                                                                                                                                           gaatcagctgagactgtgcagccataatcgaggacgcgcccgataatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGAGTTCCAGTAACTTTGGCGACTGGCACTGACAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGCCCACGCTTACTGTAGAAc                                 >  2:525834‑M2/48‑249 (MQ=255)                                                                                                                                                                                                                 gatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTaa                          >  1:587471‑M2/42‑250 (MQ=255)                                                                                                                                                                                                                          tggcagccataatcgagtgcgcgcccgctaatAGTTTGCCATATATTATTGGAATAGCGGTAAATCAGCTAAGTCTTAATACCACCAAACGGCACATTTAATCCAAATGAGATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTTAAAAATGTCGTACCGTACTTTCCAAGCTGCTTAGAGATGCTTATTTGTTCCCGACCACGCTTACTGTATAACAGATTAAAATAATCaa                 <  2:837062‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                             cagccataatcgaggacgcgcccggtaatGTTTTGCCATATATTATTGGAATAGCTGTAAATCAGCGAAGTCGTAATATCTCCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTTTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTCCTGTAGAACAGATTAAAATAATCAATaaa             <  1:607662‑M2/222‑1 (MQ=255)                                                                                                                                                                                                                                ccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATtg          <  1:362168‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                ccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGATCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATtg          >  1:54102‑M2/27‑251 (MQ=255)                                                                                                                                                                                                                                    aatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTgtgt       <  1:528535‑M2/228‑1 (MQ=255)                                                                                                                                                                                                                                       cgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTg    <  1:215089‑M2/231‑1 (MQ=255)                                                                                                                                                                                                                                        gaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCt  <  1:289003‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                                      taatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTc                                                       >  2:278242‑M2/5‑184 (MQ=255)                                                                                                                                                                                                                                                      taatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTc                                                       <  1:278242‑M2/180‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         AAGAGAACTGTAACCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATGCCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTAAGCGTAAAAGCGAGTAAATGATCCCGATCGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCT  >  NC_000913/1589477‑1589958

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