Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,001,722	Δ1 bp	coding (2507/2601 nt)	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,719	0	G	.	100.0%	81.5 / NA	20	coding (2504/2601 nt)	elfC	putative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (8/12);  total (8/12)

CTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGATTACGC  >  NC_000913/1001641‑1001793                                                                               |                                                                           cTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTggg                                                                          <  2:205612/81‑1 (MQ=255)    aTAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTgggc                                                                         <  2:205613/79‑2 (MQ=255)     tAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTgggcgaa                                                                      <  1:205614/81‑5 (MQ=255)           cATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTccc                                                               <  2:205615/81‑1 (MQ=255)               aTTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATc                                                           <  1:205617/81‑1 (MQ=255)               aTTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATc                                                           <  1:205618/81‑1 (MQ=255)               aTTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATc                                                           <  2:205616/81‑1 (MQ=255)                  gatgatgatGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCaa                                                         <  1:205619/80‑1 (MQ=255)                  gatgatgatGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCaa                                                         <  1:205620/80‑1 (MQ=255)                      atgatgGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAAttt                                                     <  1:205621/80‑1 (MQ=255)                           ggCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCata                                               >  2:205651/1‑81 (MQ=255)                            gctcgctcTATATGGCGTGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCatat                                              <  1:920230/81‑1 (MQ=255)                              tcgctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATc                                             >  2:205654/1‑80 (MQ=255)                              tcgctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATc                                             >  1:205644/1‑80 (MQ=255)                                   ctATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTAc                                       >  2:205657/1‑81 (MQ=255)                                   ctATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTAc                                       >  1:205656/1‑81 (MQ=255)                                                  aCCGGCGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGaacaa                        >  1:205655/1‑81 (MQ=255)                                                       cGAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAAt                   >  2:205646/1‑81 (MQ=255)                                                        gAAGGGAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAAtt                  <  1:205622/81‑1 (MQ=255)                                                            ggAACAATAAGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAAct              <  2:205623/81‑1 (MQ=255)                                                                     aGCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGAtt      <  1:205624/80‑1 (MQ=255)                                                                      gCGTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGATTAc    <  1:205625/81‑1 (MQ=255)                                                                        gTGCGCT‑GGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGATTAcgc  >  1:205665/1‑81 (MQ=255)                                                                               |                                                                           CTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGATTACGC  >  NC_000913/1001641‑1001793

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