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

TGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAA  >  NC_000913/4548942‑4549436                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCtt                                                                                                                                                                                                                                                       <  2:559074/250‑1 (MQ=255)  gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTa                                                                                                                                                                                                                                                      <  2:420810‑M1/250‑2 (MQ=255)        gtcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctgg                                                                                                                                                                                                                                               <  1:399810‑M1/251‑9 (MQ=255)               aTGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcatta                                                                                                                                                                                                                                        <  2:199302‑M1/251‑16 (MQ=255)                 gtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaa                                                                                                                                                                                                                                       <  2:124855‑M1/250‑17 (MQ=255)                     ggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagcttgctcattaattgcc                                                                                                                                                                                                                                  >  2:136356‑M1/1‑230 (MQ=255)                      gCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccg                                                                                                                                                                                                                                 <  1:646799‑M1/251‑23 (MQ=255)                        taaacCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACAATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGAAGTAGCTATCCATTTCCTACCACCAGCGCAACGCCGCGCGTTGTTTTTAATTCCATAACGGATAAGCCGTGGCCGGCGGCGCTTtaagctggctcattaattgcagtt                                                                                                                                                                                                                               >  2:69930‑M1/2‑227 (MQ=255)                                     ggATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgac                                                                                                                                                                                                                  <  1:629416‑M1/251‑38 (MQ=255)                                      gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgctttttttgcgac                                                                                                                                                                                                                  >  2:314948‑M1/1‑213 (MQ=255)                                          tttCGACGGAACTCTTTTGCCATAACGCTTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGCAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGCAACGCCCCGCGTTGTTTATAATTCGAGAAAGGATAAGCCCTGGCCGGGGGCGCTTaaagctggctcatttatttccgtgattattttgcgacagacc                                                                                                                                                                                                             >  2:652760‑M1/1‑209 (MQ=255)                                              gACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaac                                                                                                                                                                                                          >  1:507252‑M1/1‑205 (MQ=255)                                                 gCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGCCCGAAAACTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGCGGACGGATAAGCCGGGGCCGGTGACGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaact                                                                                                                                                                                                      >  1:36028‑M1/1‑202 (MQ=255)                                                   aaaTCTTTTTCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTTTTATCTAATTTTTCCGTGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAGTCGAGAACGGATAAGCCGTGTCCGGTGGCGCTTaaagctggctcattaattgccgggcttattttgctacagaccaacaactat                                                                                                                                                                                                    <  2:559771‑M1/251‑52 (MQ=255)                                                                             gAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtg                                                                                                                                                                          >  1:10784‑M1/1‑174 (MQ=255)                                                                               aaCCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtgga                                                                                                                                                                        <  1:3525‑M1/251‑80 (MQ=255)                                                                                  cATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaata                                                                                                                                                                     <  2:47778‑M1/251‑83 (MQ=255)                                                                                             aTGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaat                                                                                                                                                          <  2:581748‑M1/251‑94 (MQ=255)                                                                                             aTGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgcca                                                                                                                                                            <  1:428457‑M1/249‑92 (MQ=255)                                                                                                 cacaCTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTAACACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatg                                                                                                                                                      >  2:435488‑M1/1‑154 (MQ=255)                                                                                                  acacTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTCAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAACCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTaaagctggctcattacttgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatttttacgacaataata                                                                                                                                                      >  1:307820‑M1/1‑153 (MQ=255)                                                                                                    acTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttccgccaataatgatg                                                                                                                                                   >  2:130349‑M1/1‑151 (MQ=255)                                                                                                        cAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggg                                                                                                                                               >  1:507341‑M1/1‑147 (MQ=255)                                                                                                            gAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtg                                                                                                                                           >  1:630893‑M1/1‑143 (MQ=255)                                                                                                                 tCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctac                                                                                                                                      >  1:234418‑M1/1‑138 (MQ=255)                                                                                                                  cgtcttaTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacatcgatgatgtccagtttgtgtggaatatttgcgccaataaggatgtggtggttcctac                                                                                                                                      <  2:406897‑M1/247‑114 (MQ=255)                                                                                                                            cggcgTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAAGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggcatatttacgcaaataatgatgtggtggtgcctactggcggctgcg                                                                                                                           >  1:187905‑M1/1‑127 (MQ=255)                                                                                                                              gcgTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgat                                                                                                                         <  2:630893‑M1/251‑127 (MQ=255)                                                                                                                               cgTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGGTTATAATTCGTGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttatttttcgacagaccaacaactataacagcggtgatttccagtttgtgtgacctatttacgccaataatgctgaggtggtgcctattggcggctggggtt                                                                                                                        >  1:9926‑M1/1‑124 (MQ=255)                                                                                                                                                 cGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtca                                                                                                      <  1:140977‑M1/251‑146 (MQ=255)                                                                                                                                                          aaaaTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactc                                                                                             >  1:449296‑M1/1‑97 (MQ=255)                                                                                                                                                           aaaTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactct                                                                                            <  1:130349‑M1/251‑156 (MQ=255)                                                                                                                                                               ataGTGGCAGTAGCTATCAATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccg                                                                                        >  1:519110‑M1/1‑92 (MQ=255)                                                                                                                                                                 aGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatggcaccggtaatctgccgg                                                                                       >  1:198737‑M1/1‑90 (MQ=255)                                                                                                                                                                    ggCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactc                                                                                   >  2:263017‑M1/1‑87 (MQ=255)                                                                                                                                                                    ggCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggact                                                                                    >  1:472793‑M1/1‑87 (MQ=255)                                                                                                                                                                    ggCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggact                                                                                    <  1:214721‑M1/250‑164 (MQ=255)                                                                                                                                                                     gCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactac                                                                                  <  2:618495‑M1/251‑166 (MQ=255)                                                                                                                                                                                       cTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgagatgtttctgcccgtgatgtcaccggtactctgccggcctaccctggttcagtgccaatt                                                                >  2:115780‑M1/1‑68 (MQ=255)                                                                                                                                                                                          ccaccaGCGAAACGCCGCGCGTTGTTTATAATTCGCGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcct                                                             >  1:220412‑M1/1‑65 (MQ=255)                                                                                                                                                                                           caccaGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctc                                                            <  1:652760‑M1/251‑188 (MQ=255)                                                                                                                                                                                                    aaCGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgttt                                                   <  2:372782‑M1/251‑197 (MQ=255)                                                                                                                                                                                                    aaCGCCGAGCGTTGTTTATAATTCGAGAACGGATACGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgattattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggatgcgatgtttctgctcgtggtgtcaccgttaatctgcccggataccctcggttattgcccattccctttaccgttt                                                   >  2:81916‑M1/1‑55 (MQ=255)                                                                                                                                                                                                       gccgcGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactcccctggttcagtgcccattcctcttaccgtttatt                                                >  1:322096‑M1/1‑52 (MQ=255)                                                                                                                                                                                                        ccgcTCGTTGTTTATAATTCGAGAACGTATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtgtgtgtggaatatttacgccgataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtccccgtttctctgccggactcccctggttcagtgccaattcctcttaccgtttattg                                               <  1:94059‑M1/251‑201 (MQ=255)                                                                                                                                                                                                        ccgcGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgttgtcaccgatactctgccggactaccctgggtcagtgccaattcctcttaccgtttattt                                               >  2:73994‑M1/1‑51 (MQ=255)                                                                                                                                                                                                               ttgttTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaa                                         >  1:235959‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                               aaCNGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactcttccggactacccgagttcagtgccaattcctcttaccgtttattgtgcgaaaagccaaaacctggggt                        <  2:605601‑M1/251‑224 (MQ=255)                                                                                                                                                                                                                                 cGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaaaagccaaaacctggggtat                      <  2:10784‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                                    aTAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaaaatccaaaacctggggtattac                   <  2:22059‑M1/251‑229 (MQ=255)                                                                                                                                                                                                                                      aaGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaaaagccaaaacctggggtattacc                  <  1:263017‑M1/250‑230 (MQ=255)                                                                                                                                                                                                                                         acgTGGCCGGTGGCGCTTaaagctggctcattaattgaagtgcttattttgcgacaggccaacaactataacagcgatgatttccagtttgtgtggcatatttacgccaataatgatgtggtggtgcctcctggaggcatcgatgtttctgcttgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctctttacgtttattgtgcgacaagccaaaacccggggaattaacnccc              >  2:595387‑M1/2‑18 (MQ=255)                                                                                                                                                                                                                                              gCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaaaagccaaaacctggggtattacctctccggca         <  2:472793‑M1/251‑239 (MQ=255)                                                                                                                                                                                                                                                   tGGCGCTTaaagctggctcattaattgccgtgcttattttgcaacagaccaacaactataacagcgatgctttccagtttgtgtggaatatttacgccaatactgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctgggtcagtgccaattcctcttaccgtttattgtgcgcaaagcccaaacctggggtattacctctccggcacaaac    >  2:208407‑M1/1‑8 (MQ=255)                                                                                                                                                                                                                                                    tgcgctTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaaaagccaaaacctggggtattacctctccggcacaacc    <  2:465782‑M1/249‑244 (MQ=255)                                                                                                                                                                                                                                                        cTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgcaaagccaaaacctggggtattacctctccggcacaaccgc  >  1:392573‑M1/1‑3 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      TGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAA  >  NC_000913/4548942‑4549436

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