Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,001,251	A→C	Q679P (CAA→CCA)	elfC →	putative fimbrial usher protein ElfC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,001,251	0	A	C	100.0%	72.1 / NA	20	Q679P (CAA→CCA)	elfC	putative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/9);  total (11/9)

TTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGAT  >  NC_000913/1001024‑1001486                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            ttACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGc                                                                                                                                                                                                                       >  1:73237/1‑250 (MQ=255)                 tctAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTc                                                                                                                                                                                                      >  2:164899/1‑251 (MQ=255)                                      gAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCCTTCACAGGGCATTACCCTGGGGc                                                                                                                                                                                  >  2:56168/1‑250 (MQ=255)                                                      aaaaCCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACt                                                                                                                                                                 >  2:18209/1‑251 (MQ=255)                                                                            aTTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAg                                                                                                                                           <  2:73237/251‑1 (MQ=255)                                                                                     tGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCag                                                                                                                                     >  1:110581/1‑248 (MQ=255)                                                                                            aaCCTGAGCTATAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTc                                                                                                                            <  1:44334/250‑1 (MQ=255)                                                                                                                 cAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCt                                                                                                      <  1:188692/251‑1 (MQ=255)                                                                                                                                  ggggACAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCCCAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGcgcg                                                                                     >  1:156737/1‑251 (MQ=255)                                                                                                                                                   gTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGaaa                                                                                                  <  1:172456/221‑1 (MQ=255)                                                                                                                                                   gTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGaaa                                                                                                  >  2:172456/1‑221 (MQ=255)                                                                                                                                                         ccAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGcc                                                              >  1:71055/1‑251 (MQ=255)                                                                                                                                                         ccAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGcc                                                              <  2:46217/251‑1 (MQ=255)                                                                                                                                                              aTGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCActtc                                                         >  1:101604/1‑251 (MQ=255)                                                                                                                                                                                 ttGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAaatcgaat                                       <  1:222333/250‑1 (MQ=255)                                                                                                                                                                                              aGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATgcgg                         <  1:167018/251‑1 (MQ=255)                                                                                                                                                                                                 gggCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGcgt                      <  1:164899/251‑1 (MQ=255)                                                                                                                                                                                                            aCAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAaac             <  1:196148/249‑1 (MQ=255)                                                                                                                                                                                                                 gATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAaacgtaacgt      >  1:1173/1‑251 (MQ=255)                                                                                                                                                                                                                     aCGGCAGTCAACAACCACTGAACTACGCTCTTTCCGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGAt  >  2:153851/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            TTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGAT  >  NC_000913/1001024‑1001486

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