Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,272	A→C	T393P (ACC→CCC)	fimD →	type I fimbriae usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,546,272	0	A	C	100.0%	121.3 / NA	43	T393P (ACC→CCC)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/24);  total (19/24)

TAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGA  >  NC_000913/4546034‑4546501                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       tAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATATATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAACGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTAAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGAAGCAGGAAAAACCCCGCTTTTTcc                                                                                                                                                                                                                            >  1:592940/1‑251 (MQ=255)  aaaCAAAATGGGTCTTACATTTATAATAGTACGGTGCCACCGGTGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCa                                                                                                                                                                                                                           <  2:238974/251‑1 (MQ=255)  aaaCAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTcc                                                                                                                                                                                                                            >  2:191324/1‑250 (MQ=255)       aaaTGGGTATGACATTTATAATAGTACGGTGACACCGGGGCCTTTTACCAGCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTa                                                                                                                                                                                                                      <  2:515463/251‑1 (MQ=255)           gggTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGAAGGTAATAGGGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTAcctt                                                                                                                                                                                                                  >  2:490710/1‑248 (MQ=255)               aTGACATTTATAATAGTCCGGCGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTACTAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGACGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTc                                                                                                                                                                                                              >  2:289523/1‑251 (MQ=255)                        ataatcatACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTATTAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGGGGATACCGTAGTGGAAATGCGCAGCAGGATAAACCCCGCTTTTTCCAGAGTACATTACTCCCCGGCCtt                                                                                                                                                                                                     >  2:288936/8‑251 (MQ=255)                                                 ttttACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATAtggtg                                                                                                                                                                            >  1:478031/1‑251 (MQ=255)                                                                 atCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCgg                                                                                                                                                             >  1:119652/1‑250 (MQ=255)                                                                  tCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGGGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGATGGACAAAATATGGTGGAACCCAACTGGCgggg                                                                                                                                                           >  2:469044/1‑249 (MQ=255)                                                                         cGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCGGGCTGGACAATATATGGTGGAACGCAACGGGCGG‑ATCGTTATc                                                                                                                                                   >  2:35462/1‑251 (MQ=255)                                                                                              gCAGGTAACGATCAAAGAGGCGGACGGCATCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGTAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCg                                                                                                                              <  2:139325/251‑1 (MQ=255)                                                                                                                  cTGACGGCAGCACGAAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGACTACCGTAGTGGAAATCCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACGCGACGGCCGTCCGGCTGGATGGACCATATATGGAGGAACTAAACTGGCGGGAT‑GTTGGCGTGGTTTTAATTTCGGGATTGGGAAAAGCTTGGGGGGACTg                                                                                                          >  2:333275/1‑251 (MQ=255)                                                                                                                    gACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTg                                                                                                          <  1:5759/249‑1 (MQ=255)                                                                                                                      cGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGgcg                                                                                                      <  1:160390/251‑1 (MQ=255)                                                                                                                         aagcccGCAGATTTTTACCGTACCCGATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATCCCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTGCATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGAGctc                                                                                                   <  2:506028/246‑1 (MQ=255)                                                                                                                           gcaCGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCtctg                                                                                                 <  2:142432/251‑1 (MQ=255)                                                                                                                              cGCAGATTTTTTCCGTGCCCTATTCGTCATTCCCGCTTTTGCAACGTGAAGGGCATACTCGTGATTACATTACGTCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCTGGAAAAACATGTGGGCACTGGGCACgctgtct                                                                                              <  2:553615/251‑1 (MQ=255)                                                                                                                                 aGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGAAACGTGACGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAAAGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCtgtg                                                                                           >  2:545233/1‑251 (MQ=255)                                                                                                                                 aGATTTTTACCGTACCCTATTCGTCAGGCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATGACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTGAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCtgtg                                                                                           <  1:61960/251‑1 (MQ=255)                                                                                                                                         aCCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCGTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGAc                                                                                   <  1:386646/251‑1 (MQ=255)                                                                                                                                         aCCGTACCCTATGCGTCAGTCCCGCATTTGCAACGTTAAGTGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATTCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCTG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCGCTGGGCTCTCTGTCTGTGGATATGAc                                                                                   <  2:251097/251‑1 (MQ=255)                                                                                                                                          ccGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGTTCTGTCTTTGGATATGACg                                                                                  <  2:412164/251‑1 (MQ=255)                                                                                                                                                ccTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCGCGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCt                                                                            <  2:335848/251‑1 (MQ=255)                                                                                                                                                       gTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                     <  1:25727/251‑1 (MQ=255)                                                                                                                                                         cAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACCATATATGGTGGAACGCAACTGGCGG‑GTCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTTGGCGCTCTGTCTGTGGATAGGACGCAGGCTAATTCcac                                                                    >  2:154430/1‑250 (MQ=255)                                                                                                                                                         cAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCTCTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTccgca                                                                   <  2:113304/251‑4 (MQ=255)                                                                                                                                                         cAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGatcgcccgcagcatggag                                                                   <  2:119652/251‑18 (MQ=255)                                                                                                                                                            tCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACtt                                                                >  1:188398/1‑251 (MQ=255)                                                                                                                                                             cccGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACtt                                                                <  1:469044/250‑1 (MQ=255)                                                                                                                                                                    ttGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATg                                                         <  1:591130/250‑1 (MQ=255)                                                                                                                                                                       cAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGCGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGGGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGGCGCAGGCTAATTCCACACTTCCCGATGACa                                                      >  2:432307/1‑250 (MQ=255)                                                                                                                                                                          cGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAA‑TGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTc                                                   >  1:130673/1‑249 (MQ=255)                                                                                                                                                                                 ggCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGAc                                           <  1:443095/251‑1 (MQ=255)                                                                                                                                                                                 ggCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGAc                                           >  2:520933/1‑251 (MQ=255)                                                                                                                                                                                  gCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGTGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACg                                          <  1:616736/251‑1 (MQ=255)                                                                                                                                                                                    aTACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACTCAACTGGCTG‑ATCGTTTTCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGa                                        <  1:555774/251‑1 (MQ=255)                                                                                                                                                                                    aTACTCGTTACTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGCAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAAGCTAATTCCACACTTCCCGATGACAGTCAGCATGACggc                                        >  2:184882/1‑250 (MQ=255)                                                                                                                                                                                         cGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATc                                   <  2:519138/251‑1 (MQ=255)                                                                                                                                                                                             aTTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGt                                <  1:373571/250‑1 (MQ=255)                                                                                                                                                                                                         cAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACAATACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCtat                   >  1:102168/1‑251 (MQ=255)                                                                                                                                                                                                                         tGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAACTCGCTCAATg   >  2:13861/1‑251 (MQ=255)                                                                                                                                                                                                                          ggAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGa  <  2:173286/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       TAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGG‑ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGA  >  NC_000913/4546034‑4546501

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