Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,549,192	Δ39 bp	coding (385‑423/903 nt)	fimH →	type 1 fimbriae D‑mannose specific adhesin

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4549192	4549230	39	56 [0]	[3] 58	fimH	type 1 fimbriae D‑mannose specific adhesin

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4549191	0 (0.000)	55 (0.870)	52/480	0.3	100%	coding (384/903 nt)	fimH	type 1 fimbriae D‑mannose specific adhesin
	?	NC_000913	4549231 =	0 (0.000)					coding (424/903 nt)	fimH	type 1 fimbriae D‑mannose specific adhesin

CTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCA  >  NC_000913/4548982‑4549476                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        gcgcccgtcgtgaatgtggggcaaaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgctta                                                                                                                                                                                                                                                         <  2:420810‑M2/1‑1 (MQ=255)       gtcgtgaatgtggggcaaaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTgg                                                                                                                                                                                                                                                  <  1:399810‑M2/8‑1 (MQ=255)              atgtggggcaaaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTa                                                                                                                                                                                                                                           <  2:199302‑M2/15‑1 (MQ=255)                gtggggcaaaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTaa                                                                                                                                                                                                                                          <  2:124855‑M2/16‑1 (MQ=255)                    ggcaaaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTTGCTCATTAATTGcc                                                                                                                                                                                                                                     >  2:136356‑M2/231‑251 (MQ=255)                     gcaaaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCg                                                                                                                                                                                                                                    <  1:646799‑M2/22‑1 (MQ=255)                       taaacctggtcgtggatctttcgacgcaaatcttttgccataacgattatccggaaacaattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggaagtagctatccatttcctaccaccagcgcaacgccgcgcgttgtttttaattccataacggataagccgtggccggcggcgcttTAAGCTGGCTCATTAattgcagtt                                                                                                                                                                                                                                  >  2:69930‑M2/228‑247 (MQ=255)                                    ggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGAc                                                                                                                                                                                                                     <  1:629416‑M2/37‑1 (MQ=255)                                     gatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTTTTTTGCGAc                                                                                                                                                                                                                     >  2:314948‑M2/214‑250 (MQ=255)                                         tttcgacggaactcttttgccataacgcttatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgcaaaatatagtggcagtagctatccatttcctaccacccgcgcaacgccccgcgttgtttataattcgagaaaggataagccctggccgggggcgcttAAAGCTGGCTCATTTATTTCCGTGATTATTTTGCGACAGAcc                                                                                                                                                                                                                >  2:652760‑M2/210‑251 (MQ=255)                                             gacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACcaac                                                                                                                                                                                                             >  1:507252‑M2/206‑250 (MQ=255)                                                gcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggcccgaaaactatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgcggacggataagccggggccggtgacgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACt                                                                                                                                                                                                         >  1:36028‑M2/203‑251 (MQ=255)                                                  aaatctttttccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgttttatctaatttttccgtgaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataagtcgagaacggataagccgtgtccggtggcgcttAAAGCTGGCTCATTAATTGCCGGGCTTATTTTGCTACAGACCAACAACtat                                                                                                                                                                                                       <  2:559771‑M2/51‑1 (MQ=255)                                                                            gaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTtgtgtg                                                                                                                                                                             >  1:10784‑M2/175‑251 (MQ=255)                                                                              aaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGa                                                                                                                                                                           <  1:3525‑M2/79‑1 (MQ=255)                                                                                 cattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAata                                                                                                                                                                        <  2:47778‑M2/82‑1 (MQ=255)                                                                                            atgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCaat                                                                                                                                                             <  2:581748‑M2/93‑1 (MQ=255)                                                                                            atgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCa                                                                                                                                                               <  1:428457‑M2/91‑1 (MQ=255)                                                                                                cacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaacaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccgggggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAatg                                                                                                                                                         >  2:435488‑M2/155‑251 (MQ=255)                                                                                                 acactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtcaaatatagtggcagtagctatccatttcctaccaccagcgaaaacccgcgcgttgtttataattcgagaacggataagccgtggccgggggcgcttAAAGCTGGCTCATTACTTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATTTTTACGAcaataata                                                                                                                                                         >  1:307820‑M2/154‑249 (MQ=255)                                                                                                   actgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTCCGCCAATAatgatg                                                                                                                                                      >  2:130349‑M2/152‑251 (MQ=255)                                                                                                       caacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGAtgtggg                                                                                                                                                  >  1:507341‑M2/148‑250 (MQ=255)                                                                                                           gaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATgtggtggtg                                                                                                                                              >  1:630893‑M2/144‑251 (MQ=255)                                                                                                                tcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTAc                                                                                                                                         >  1:234418‑M2/139‑251 (MQ=255)                                                                                                                 cgtcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACATCGATGATGTCCAGTTTGTGTGGAATATTTGCGCCAATAAGGATGTGGTGGTTCCTAc                                                                                                                                         <  2:406897‑M2/113‑1 (MQ=255)                                                                                                                           cggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaaagccgcgcgttgtttataattcgagaacggataagccgtggccgggggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGCATATTTACGCAAATAATGATGTGGTGGTGCCTACTGGCGGCTGCg                                                                                                                              >  1:187905‑M2/128‑251 (MQ=255)                                                                                                                             gcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGAt                                                                                                                            <  2:630893‑M2/126‑1 (MQ=255)                                                                                                                              cgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttggttataattcgtgaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTTCGACAGACCAACAACTATAACAGCGGTGATTTCCAGTTTGTGTGACCTATTTACGCCAATAATGCTGAGGTGGTGCCTATTGGcggctggggtt                                                                                                                           >  1:9926‑M2/125‑246 (MQ=255)                                                                                                                                                cgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCa                                                                                                         <  1:140977‑M2/145‑1 (MQ=255)                                                                                                                                                         aaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTActc                                                                                                >  1:449296‑M2/98‑251 (MQ=255)                                                                                                                                                          aaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTActct                                                                                               <  1:130349‑M2/155‑1 (MQ=255)                                                                                                                                                              atagtggcagtagctatcaatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccgggggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCg                                                                                          >  1:519110‑M2/93‑251 (MQ=255)                                                                                                                                                                agtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGGCACCGGTAATCTG‑CCgg                                                                                         >  1:198737‑M2/91‑250 (MQ=255)                                                                                                                                                                   ggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGActc                                                                                     >  2:263017‑M2/88‑250 (MQ=255)                                                                                                                                                                   ggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACt                                                                                      >  1:472793‑M2/88‑250 (MQ=255)                                                                                                                                                                   ggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACt                                                                                      <  1:214721‑M2/163‑1 (MQ=255)                                                                                                                                                                    gcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTAc                                                                                    <  2:618495‑M2/165‑1 (MQ=255)                                                                                                                                                                                      ctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGAGATGTTTCTGCCCGTGATGTCACCGGTACTCTG‑CCGGCCTACCCT‑GGTTCAGTGCCAAtt                                                                 >  2:115780‑M2/69‑251 (MQ=255)                                                                                                                                                                                         ccaccagcgaaacgccgcgcgttgtttataattcgcgaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCct                                                              >  1:220412‑M2/66‑251 (MQ=255)                                                                                                                                                                                          caccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCctc                                                             <  1:652760‑M2/187‑1 (MQ=255)                                                                                                                                                                                                   aacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGttt                                                    <  2:372782‑M2/196‑1 (MQ=255)                                                                                                                                                                                                   aacgccgagcgttgtttataattcgagaacggatacgccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGATTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGATGCGATGTTTCTGCTCGTGGTGTCACCGTTAATCTGCCCGGA‑TACCCTCGGTT‑ATTGCCCATTCCCTTTACCGttt                                                    >  2:81916‑M2/56‑251 (MQ=255)                                                                                                                                                                                                      gccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTCCCCT‑GGTTCAGTGCCCATTCCTCTTACCGTTTAtt                                                 >  1:322096‑M2/53‑251 (MQ=255)                                                                                                                                                                                                       ccgctcgttgtttataattcgagaacgtataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTGTGTGTGGAATATTTACGCCGATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCCCCGTTTCTCTG‑CCGGACTCCCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATtg                                                <  1:94059‑M2/200‑1 (MQ=255)                                                                                                                                                                                                       ccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGTTGTCACCGATACTCTG‑CCGGACTACCCT‑GGGTCAGTGCCAATTCCTCTTACCGTTTattt                                                >  2:73994‑M2/52‑250 (MQ=255)                                                                                                                                                                                                              ttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGaa                                          >  1:235959‑M2/45‑250 (MQ=255)                                                                                                                                                                                                                              aacngataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTT‑CCGGACTACCCG‑AGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGt                         <  2:605601‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                                cggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTAt                       <  2:10784‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                   ataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAATCCAAAACCTGGGGTATTAc                    <  2:22059‑M2/228‑1 (MQ=255)                                                                                                                                                                                                                                     aagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAcc                   <  1:263017‑M2/229‑1 (MQ=255)                                                                                                                                                                                                                                        acgtggccggtggcgcttAAAGCTGGCTCATTAATTGAAGTGCTTATTTTGCGACAGGCCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGCATATTTACGCCAATAATGATGTGGTGGTGCCTCCTGGAGGCATCGATGTTTCTGCTTGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTTACGTTTATTGTGCGACAAGCCAAAACCCGGggaattaacnccc               >  2:595387‑M2/19‑245 (MQ=255)                                                                                                                                                                                                                                             gccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGca          <  2:472793‑M2/238‑1 (MQ=255)                                                                                                                                                                                                                                                  tggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCAACAGACCAACAACTATAACAGCGATGCTTTCCAGTTTGTGTGGAATATTTACGCCAATACTGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGGTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGCAAAGCCCAAACCTGGGGTATTACCTCTCCGGCacaaac     >  2:208407‑M2/9‑249 (MQ=255)                                                                                                                                                                                                                                                   tgcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAAcc     <  2:465782‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                                       cttAAAGCTGGTTCATTAATTGCCGTGCTTATCTTTCGTCAGACCAACTACAATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTATGCCAAGAATGATGTGCTGGTGCCTACTGGCGGCTGCGCTGTTTCTGCACGTGATGTAACCGTTCCTCTG‑CCGTACCACCCT‑GGTTCAGTGCTAATTCCTCCTACCGTTTATAATGCGAAAAGCCAAAACGTGGGGTATTTCCTCTCCGGtacagctgaa  <  2:78885/248‑7 (MQ=255)                                                                                                                                                                                                                                                       cttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGCAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGc   >  1:392573‑M2/4‑248 (MQ=255)                                                                                                                                                                                                                                                        ttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCctct                                                            >  2:538330/1‑190 (MQ=255)                                                                                                                                                                                                                                                        ttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATGCctct                                                            <  1:538330/190‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        CTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTG‑CCGGACTACCCT‑GGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCA  >  NC_000913/4548982‑4549476

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