Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	309,218	Δ1 bp	coding (115/2526 nt)	ecpC ←	putative fimbrial usher protein EcpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,215	0	C	.	100.0%	170.5 / NA	37	coding (118/2526 nt)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (19/18);  total (19/18)

CTGGCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGC  >  NC_000913/309135‑309290                                                                                 |                                                                            cTGGCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCAc                                                                             <  2:67386/81‑1 (MQ=255)      tACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCAcccacc                                                                        >  2:67428/1‑78 (MQ=255)          ggCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTAtt                                                                   >  2:67433/1‑81 (MQ=255)           gCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTAttt                                                                  >  1:67431/1‑81 (MQ=255)           gCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTAttt                                                                  <  1:67387/81‑1 (MQ=255)           gCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTAttt                                                                  <  2:67388/81‑1 (MQ=255)            cGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTtg                                                                 >  1:67435/1‑81 (MQ=255)              agaTGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTtgct                                                               >  2:67443/1‑81 (MQ=255)               gaTGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTtgctg                                                              >  2:67446/1‑81 (MQ=255)                   aaTATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTgcgc                                                           <  1:67389/80‑1 (MQ=255)                   aaTATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCa                                                          <  2:67390/81‑1 (MQ=255)                        agagCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTAt                                                     >  2:67441/1‑81 (MQ=255)                        agagCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTAt                                                     >  2:67438/1‑81 (MQ=255)                        agagCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTAt                                                     >  2:67437/1‑81 (MQ=255)                            cGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTca                                                 <  2:67391/81‑1 (MQ=255)                                   cTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAgcgc                                          <  2:67393/81‑1 (MQ=255)                                   cTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAgcgc                                          <  1:67392/81‑1 (MQ=255)                                      aTGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTgg                                       >  2:67436/1‑81 (MQ=255)                                      aTGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTg                                        <  1:67394/80‑1 (MQ=255)                                       tGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGc                                      <  2:67395/81‑1 (MQ=255)                                       tGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGc                                      <  2:67396/81‑1 (MQ=255)                                         ccGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCAt                                    >  1:67434/1‑81 (MQ=255)                                         ccGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCAt                                    <  1:67397/81‑1 (MQ=255)                                          cGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATc                                   <  1:67398/81‑1 (MQ=255)                                           gTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCg                                  <  2:67400/81‑1 (MQ=255)                                           gTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCg                                  <  1:67399/81‑1 (MQ=255)                                              cTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCggg                                <  1:67401/80‑1 (MQ=255)                                                gAAGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCt                              >  2:67432/1‑80 (MQ=255)                                                  aGCGCCTGACTGAAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTgaa                           >  2:67452/1‑81 (MQ=255)                                                             gaAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACaaa                   >  2:67442/1‑78 (MQ=255)                                                             gaAGGCCTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACaaa                   >  1:67440/1‑78 (MQ=255)                                                                  ccTGCGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCAt           >  1:67451/1‑81 (MQ=255)                                                                      cGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGcc        <  1:67403/80‑1 (MQ=255)                                                                      cGGAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGcc        <  2:67402/80‑1 (MQ=255)                                                                       ggAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCg       >  2:67445/1‑80 (MQ=255)                                                                       ggAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGa      <  2:67404/81‑1 (MQ=255)                                                                        gAATAATCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGTCATCGGGCTGAACGAACAAAAAGACCATGCCg       >  1:67439/1‑79 (MQ=255)                                                                         aataatCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCTAAg    >  2:67448/1‑81 (MQ=255)                                                                            aatCA‑CCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGc  >  2:67444/1‑80 (MQ=255)                                                                                 |                                                                            CTGGCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGC  >  NC_000913/309135‑309290

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