Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,591	G→A	G499D (GGC→GAC)	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,591	0	G	A	100.0%	87.5 / NA	26	G499D (GGC→GAC)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCT  >  NC_000913/4546348‑4546825                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           aTCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCGCTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGCGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGAGGGTTACCGTTATTCGACCAGCGGATATTTTACTTTCGCTGATACAACATACAGTAGAATGAATGACTacaac                                                                                                                                                                                                                                     >  1:238319/1‑251 (MQ=255)     ttATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCg                                                                                                                                                                                                                                 >  1:236144/1‑251 (MQ=255)          gTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAac                                                                                                                                                                                                                             >  2:79374/1‑250 (MQ=255)                        ggTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTAt                                                                                                                                                                                                              >  2:273984/1‑251 (MQ=255)                                    aaCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                                                                                                                                                                                  <  1:202160/251‑1 (MQ=255)                                              cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                                                                                                                                                        >  2:104794/1‑251 (MQ=255)                                                    gcgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCACGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGAGACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACTGAGTTATTCAGGTTAAGCATAAATTCACCGACTAt                                                                                                                                                                                  <  2:56204/251‑1 (MQ=255)                                                     cgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                 <  1:329146/251‑1 (MQ=255)                                                         ctgtctgtGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                                                                                                                                                               <  1:348831/249‑1 (MQ=255)                                                                  gATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                                                                                                                                                                    <  1:156333/251‑1 (MQ=255)                                                                          gCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                                                                                                                                                            <  2:215418/251‑1 (MQ=255)                                                                           cAGGCTAATTCCACACTTCCCGATGCCAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCACAATGAAACAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAATTTCCCCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                           >  2:309696/1‑251 (MQ=255)                                                                                   ttCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt                                                                                                                                                   >  1:222417/1‑251 (MQ=255)                                                                                                              tGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCg                                                                                                                        >  1:245670/1‑251 (MQ=255)                                                                                                                                   tctcTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAGTCACCGACTATTACAACCTCGCTTAGAACAAACGCGTGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAt                                                                                                    <  2:20705/250‑1 (MQ=255)                                                                                                                                            cAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGt                                                                                           >  2:281668/1‑250 (MQ=255)                                                                                                                                               aacgcgcAATGAATCAGGTACGCATCTTCAGTTAGTGGGATACCGTTATTCGACCAGCGGATATTTGAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGCCTATTACAACCTCGCTTATAACAACAGCGGGAAATTACAACTCACCGTTACTCATCAACTCGGGCGCAACTCAACACTGTATTTGAGTGGTAGcc                                                                                       >  2:161026/7‑251 (MQ=255)                                                                                                                                                              aGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                                                        >  1:186993/1‑251 (MQ=255)                                                                                                                                                                                                   aCCAACGGATATTTTAATTTCGCTGATACAACATACAGGCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAtt                                   <  2:199931/251‑1 (MQ=255)                                                                                                                                                                                                        cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGCAATGTCGATGAGCAATTCCAGGCTGGATTAAAt                               >  1:85116/1‑250 (MQ=255)                                                                                                                                                                                                           atatTTTAATTTCGCTGATGCAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCATCTCGGGCGCACATCAACACTGTATTTGAGTGTTAGCCATCAAACTTATTGGGGAACGGGTAATGTCGATGAGCAATTCCAGCCTGGATTAAATACTg                           <  2:238319/251‑1 (MQ=255)                                                                                                                                                                                                             atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGc                          <  1:193903/250‑1 (MQ=255)                                                                                                                                                                                                             atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGc                          <  1:10348/250‑1 (MQ=255)                                                                                                                                                                                                                      tCGCTGATACAACATACAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGa                 <  1:163009/250‑1 (MQ=255)                                                                                                                                                                                                                               caacaTACAGTCGAATGAATGAAAACAACATCGAAACAAAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTgg       >  2:244997/1‑251 (MQ=255)                                                                                                                                                                                                                                     aCAGTCGAATGAATGACTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt  <  1:299/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           ATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCT  >  NC_000913/4546348‑4546825

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