Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	557,194	A→G	Q681R (CAG→CGG)	sfmD →	fimbrial biogenesis usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	557,194	0	A	G	100.0%	61.5 / NA	18	Q681R (CAG→CGG)	sfmD	fimbrial biogenesis usher protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/5);  total (13/5)

AGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGA  >  NZ_CP009273/557061‑557326                                                                                                                                      |                                                                                                                                     aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCt                                                                                                                                   <  2:234453/137‑1 (MQ=255)     aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTgg                                                                                                                             >  2:351497/1‑139 (MQ=255)                             cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg                                                                                                     >  2:27631/1‑139 (MQ=255)                                          cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc                                                                                        >  1:525285/1‑139 (MQ=255)                                           acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcc                                                                                       <  2:525285/139‑1 (MQ=255)                                                         ggCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGa                                                                          >  1:27562/1‑138 (MQ=255)                                                                          cTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag                                                                                         <  2:102659/106‑1 (MQ=255)                                                                          cTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag                                                                                         >  1:102659/1‑106 (MQ=255)                                                                                        gggggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAg                                          <  2:208283/139‑1 (MQ=255)                                                                                           ggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAACCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCg                                       >  2:397620/1‑139 (MQ=255)                                                                                               ggTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgtg                                    >  1:384945/1‑138 (MQ=255)                                                                                                          gACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTaa                                                        >  1:307822/1‑107 (MQ=255)                                                                                                          gACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTaa                                                        <  2:307822/107‑1 (MQ=255)                                                                                                          gACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaa                         >  1:143217/1‑138 (MQ=255)                                                                                                           aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa                       >  1:308959/1‑139 (MQ=255)                                                                                                           aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa                       >  1:190018/1‑139 (MQ=255)                                                                                                               cgATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa                    >  1:35414/1‑138 (MQ=255)                                                                                                                                aCTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGa  >  2:239974/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     AGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGA  >  NZ_CP009273/557061‑557326

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