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%	50.1 / NA	17	Q679P (CAA→CCA)	elfC	putative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/8);  total (9/8)

CCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAAC  >  NC_000913/1001029‑1001475                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 ccTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAACCGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTa                                                                                                                                                                                                      <  1:252687/251‑1 (MQ=255)    ggATGACCTCTAACACCAGTAGCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGGAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAg                                                                                                                                                                                                     <  1:232949/249‑1 (MQ=255)    ggATGACCACTAACACAAGTAAAGACAACGAAGGGCACAATACACAAAACCTGGGTTTAACGGAGACATTACTCGCTGACTGTCACCTGAGCTACAGCGTGCGACAGACATATAACAGCGAGGCGCAAACGGCTAATGGTAGCGCCAGCATGGATTACAACGGGGCGTCTGCGGAAGCCAGAGTGGGACACAACTACAGCGATAACGACAGTCAACAACCACTGAACTACGCTCTTTCACGCAGATTAGtt                                                                                                                                                                                                   >  2:257105/1‑251 (MQ=255)              tAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCAAGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTcac                                                                                                                                                                                         >  1:146644/1‑251 (MQ=255)                                              cacaAAGCCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCGACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACTGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTTCCCATTCACAGGGCATTTCCCCGGGGCAATCGCTGGGGGaa                                                                                                                                                         <  2:146644/251‑1 (MQ=255)                                                    aCCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAAc                                                                                                                                                   >  1:232856/1‑251 (MQ=255)                                                                                                   aGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAAc                                                                                                      >  1:42348/1‑249 (MQ=255)                                                                                                      gTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGAGTTTGCGGATGACCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAACCTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGAGTGGCGAACAGCAc                                                                                                 >  1:34732/1‑251 (MQ=255)                                                                                                                       cAGCGCGGGGAAAACGGCAAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTCGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTAATTGCAGCACCAGGCGCCGAGAATACTCGTGTGGCGAACCGCACCGGGATGCAAactgact                                                                                >  2:177458/1‑251 (MQ=255)                                                                                                                           gAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACCAAGGGGCGTTTGCGGATGCCCGAGTGGGATACAACTACAGAGATAACGGCAGCCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTCGTTGCCCATTCCCAGGGCATTACCCTGGGGCAATCACTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGTACCGAGATGAAAAGTGACTgg                                                                              >  2:144526/1‑249 (MQ=255)                                                                                                                           gAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATGCTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGgcg                                                                            <  2:131568/251‑1 (MQ=255)                                                                                                                                 aaaaCGGCTAATGGTAGCGCTAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGATACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCTCCGGGCTGAAAACTGACTGGCGCGGata                                                                      <  2:92352/251‑1 (MQ=255)                                                                                                                                     cGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATAcc                                                                  <  1:1157/251‑1 (MQ=255)                                                                                                                                                   gcCAGCATGGATTACAAAGTGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTGCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGTGGATATACCGTTGTTCCTTATGc                                                    <  2:11245/251‑1 (MQ=255)                                                                                                                                                                 cAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGTTATACCGTTGGTCCTTATGCCACTTCTTATCggg                                      >  2:138839/1‑251 (MQ=255)                                                                                                                                                                         cGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAaatcg                               >  1:47491/1‑250 (MQ=255)                                                                                                                                                                                                       aCAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTAGGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAaac  <  2:121459/249‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 CCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAAC  >  NC_000913/1001029‑1001475

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