Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,223	T→G	V139G (GTG→GGG)	fimH →	type 1 fimbriae D‑mannose specific adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,223	0	T	G	100.0%	178.1 / NA	58	V139G (GTG→GGG)	fimH	type 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/30);  total (28/30)

CGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAC  >  NC_000913/4548986‑4549458                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            cGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCt                                                                                                                                                                                                                                >  1:177470/1‑251 (MQ=255)   aCGCAAATCTTTTGCCATAACTATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGTTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTgg                                                                                                                                                                                                                              <  2:448957/251‑1 (MQ=255)   aCGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTgg                                                                                                                                                                                                                              <  2:142484/251‑1 (MQ=255)             tttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAAtt                                                                                                                                                                                                                    >  1:354370/1‑251 (MQ=255)              ttGCCATAACGATTATCCGGAAACCATGACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGATCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAAtt                                                                                                                                                                                                                    <  2:354370/250‑1 (MQ=255)               tGCCATAGCGATTATCCGTAAACCATTACAGACTATGGCACACTGCAACGAGGCTCGGCTGATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTGGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGCGGGCGATTACAGCTGGCTCATTAATTg                                                                                                                                                                                                                   <  2:329838/250‑1 (MQ=255)                 ccATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCg                                                                                                                                                                                                                <  1:595328/251‑1 (MQ=255)                      aCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTTCTACCACCAGCGAAACGCCGCGCGTTGTTTTTAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGGTTAAAGCTGGCTCCTTAATTGCCGGGCtt                                                                                                                                                                                                           >  2:151106/1‑251 (MQ=255)                        gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGCAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAt                                                                                                                                                                                                         >  1:485493/1‑251 (MQ=255)                              ccGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCg                                                                                                                                                                                                   <  1:153970/251‑1 (MQ=255)                                    aCCATTACAGACTATGTCACTCTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCTTGCTTATTTTGCGACAGAc                                                                                                                                                                                             <  1:564225/251‑1 (MQ=255)                                    aCCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAGATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGAc                                                                                                                                                                                             <  2:171324/251‑1 (MQ=255)                                      cATTACAGACTATGTCACACTGCACCGAGGCTTGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTCAAAAATAGTGGCAGGAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATCAGCCGTGGCCGGTGGCGCTTTATTTGACCCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCCTGCTCCTTAATTGCCGGGCTTTGTTTGCGaacgaccc                                                                                                                                                                                           >  2:129458/1‑244 (MQ=255)                                        ttACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATTGTTGCAGTATCTATCCATTTCCTACCACCAGCGAAACACCGCGCGTTGGTTATAATTCCAGAACGGATAAGCCGTGGCCGGTGGCGATTTATTTTACGCCTGTGAGCAATGCGGGCGGGGGGGCGGTTAAAGATGGCTCATTAATTGCCGTGGTtttttttggccaggacaac                                                                                                                                                                                         >  2:419693/1‑238 (MQ=255)                                           cAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGGGATTAAAGCTGGCTCATTTATTTCCGTGCTTATTTTGCGGCAGACcaacca                                                                                                                                                                                       >  1:293790/1‑248 (MQ=255)                                              aCTAGGTCACACTGAAACGAGGCTCGGCTTCAGGCGGCGTGTTATCGAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTAGAGAACGGATAAGCCGTTGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACta                                                                                                                                                                                    <  2:341109/250‑1 (MQ=255)                                               cTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAAAGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACtat                                                                                                                                                                                   >  1:17461/1‑250 (MQ=255)                                                tATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTTTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGGTTCAAGCTGGCTCCTTAATTGCCGTGGTTTTTTTGCGACAGACCAACAACtata                                                                                                                                                                                  >  2:451335/1‑250 (MQ=255)                                                             aaCGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTc                                                                                                                                                                    >  1:396256/1‑251 (MQ=255)                                                                         cTTATGGCGGCGGGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGGGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAAAAGCGATGATTTCCAGTTTGTGTgg                                                                                                                                                        >  1:369264/1‑251 (MQ=255)                                                                          ttATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTggt                                                                                                                                                       >  1:179759/1‑250 (MQ=255)                                                                             tgtcggcgTGCTATCTAATTTTTCCGGGACCGTAAAATATAGCGGCAGTAGCTATCCATTTCCTACCACCCGCGAAACGACGCGCGTTGTTTATAATTCGAGAACGGCTAAGCCGTGGCAGGTGGCGCTTTATTTGACGACTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCACCTTAATTGCCGTGCTTATTTTGCGACAGAACCACAACTCTAACCGCGATGATTTCCAGTTTGTGTGGAat                                                                                                                                                     >  2:251955/4‑250 (MQ=255)                                                                              ggcggcGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAatat                                                                                                                                                   <  1:627697/251‑1 (MQ=255)                                                                              ggcggcGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACGATAACAGCGATGATTTCCAGTTTGTGTGGAatat                                                                                                                                                   <  1:202092/251‑1 (MQ=255)                                                                                                ttttCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGCGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCCGTTTGTGTGGAATATTTACGCCAATAATGAtgt                                                                                                                                 >  1:617884/1‑251 (MQ=255)                                                                                                     cGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCaataa                                                                                                                                       >  2:415258/1‑240 (MQ=255)                                                                                                     cGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCaataa                                                                                                                                       <  1:415258/240‑1 (MQ=255)                                                                                                        gACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGc                                                                                                                         <  2:534063/251‑1 (MQ=255)                                                                                                          ccGTAAAATATAGTGGCATTAGCTATCCATTTCCGACCACCAGCGAAACTCCGCTCGTTGTTTATAATTCGAGCTCGGATAAGCCGTGGCCGGTGGTGCTTTATTTGGCGCCTGTGAGCAGTGCGTGCGGGGGGGCGATTAAATCTGGCTCATTAAGTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGTTGATTTCCAGTTTGTGTGGAATAGTTACGCCAATAATGATGTGGTGTTGCCt                                                                                                                       <  2:354197/251‑1 (MQ=255)                                                                                                            gTAAAGTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGGGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGACTGTGAGCAGTGCGGGCGGGGGGGCGATTCAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAAAAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCCATAATGTTGTGGTGGTGCCTAc                                                                                                                     >  2:541179/1‑251 (MQ=255)                                                                                                                        ggCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGCGCAGTGCGGGCGGGGGGGCGATTACAGCTGGCTAATTAATTGCCGTGCTTATTTTGCGACAGCCCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTggcggacgcg                                                                                                          >  2:146403/1‑245 (MQ=255)                                                                                                                                   tCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGt                                                                                                                                                                                                               >  1:528410/1‑138 (MQ=255)                                                                                                                                   tCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGt                                                                                                                                                                                                               <  2:528410/138‑1 (MQ=255)                                                                                                                                   tCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGGTGTTTCTGCTc                                                                                              >  2:457002/1‑251 (MQ=255)                                                                                                                                     cATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGt                                                                                            >  1:60048/1‑251 (MQ=255)                                                                                                                                           cTACCACCAGCGTAAAGGCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGCGGCCGGAGGCGCTTTATTTGACGCCTGTGAGTAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCAAATAATGATGTGGTGGTGCATACTGGCGGCTCAGATGTTTATTCTCGTGATGTc                                                                                      >  2:627887/1‑251 (MQ=255)                                                                                                                                              ccacccacgactcgccgcTCGTTGTTTATAATTCGAGAATGGATAAGCCGTGGCCGATGGCGCTTGTTATGACGCCTGTGAGCAGTGCGGGCGGGGGGGCTATTAAATCTGGCTCATTAATTGCCGTGCTTATTTTGCGACACACTAACATCTATAACAGCGATGATTTCCAGTTTGTTTGGAATATTTACGCCAAAAATGATGTGGTGGTGCCTACTAGCTGCTGCTTTGTTTATGCTCGTCATGTctcc                                                                                   <  2:382368/239‑4 (MQ=255)                                                                                                                                               caccaGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGGTTCTGCTCGTGATGTCACCg                                                                                  >  1:252142/1‑251 (MQ=255)                                                                                                                                                      gAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGGTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTGCTGCTCGTGATTTCACCGTTActct                                                                           <  2:42936/251‑1 (MQ=255)                                                                                                                                                        aaCGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGCGGCCGGTGGCGCTTTCTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGAGGAATCTTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGGGATGTCACCGTTATTCTGc                                                                         >  2:52586/1‑251 (MQ=255)                                                                                                                                                           gccgcGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCgg                                                                      <  2:177470/251‑1 (MQ=255)                                                                                                                                                                    tgttTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCt                                                              <  1:241913/250‑1 (MQ=255)                                                                                                                                                                    tgttTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCt                                                              <  1:253906/250‑1 (MQ=255)                                                                                                                                                                     gttTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTg                                                             <  1:240828/250‑1 (MQ=255)                                                                                                                                                                        tataATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGGCACCGTTACTCTGCCGGACTACCCTGGtt                                                          <  2:252142/250‑1 (MQ=255)                                                                                                                                                                           aaTTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtttt                                                                                                                                                                                                      >  1:509693/1‑107 (MQ=255)                                                                                                                                                                           aaTTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTAGTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtttt                                                                                                                                                                                                      <  2:509693/107‑1 (MQ=255)                                                                                                                                                                                gagaACGGATAAGCCGTGGCCAGTGGCGCTTTTTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTCATTGCCGTGATTATTTTGCGACAGACAAACAACTATAACAGCGATGATTTCCAGTTTGCGTGGACTATTTACGCCACTAATGATGTGGTGGTGCCTACTGGCGGCTCCGCTGTTTCTGCTCGGTATGTCACCGTTACTCTTCCGGAGTATCCTGGttatgggcca                                                  >  2:607994/1‑242 (MQ=255)                                                                                                                                                                                 agaACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCGACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAAt                                                <  2:620282/251‑1 (MQ=255)                                                                                                                                                                                  gaACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAAtt                                               <  1:457002/251‑1 (MQ=255)                                                                                                                                                                                      ggATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGGTGTTTCTGCTCGTGATGTCAACTGTACTCTGCCGGCCTACCCTGGTTCAGTGCCAATTCctc                                           >  2:531151/1‑251 (MQ=255)                                                                                                                                                                                       gATAAGCCGTGGCCGGTGGCGCTGTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCctct                                          <  2:573525/251‑1 (MQ=255)                                                                                                                                                                                             ccGTGGCCGGTGGCGCTTTATGTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGt                                    <  1:103043/251‑1 (MQ=255)                                                                                                                                                                                                tGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCGGTGCTTATTTTGCGCAAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGACTATTTAAGCCAATACTGATGTGGTGGTGCCTACTGGCGGCTGCGCTGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTCACCTGGTTCCGTGCCGATTCCTCTTACCGtttt                                 >  1:281828/1‑250 (MQ=255)                                                                                                                                                                                                           gcTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTTATGGCACCGTTACTCTGCCCGCCTACCCTGGTTCAGTGCCAATTTCTCTTACCGTTTAtttgtccaaa                       >  2:306884/1‑242 (MQ=255)                                                                                                                                                                                                                    ggcGCCTGTGAGCAGGGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAc              <  2:17461/248‑1 (MQ=255)                                                                                                                                                                                                                        ccTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTggg         <  2:485493/251‑1 (MQ=255)                                                                                                                                                                                                                               gCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAc  <  1:151106/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            CGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAC  >  NC_000913/4548986‑4549458

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