Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,054	(C)7→8	coding (1279/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	308,047	1	.	C	100.0%	75.5 / NA	23	A429G (GCG→GGG)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (11/12);  total (11/12)

TCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCT  >  NC_000913/307908‑308172                                                                                                                                             |                                                                                                                               tCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCg                                                                                                                                 <  1:362535/140‑1 (MQ=255)     aaTTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccc                                                                                                                           <  1:1074988/140‑1 (MQ=255)     aaTTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑cccc                                                                                                                           <  1:517522/140‑1 (MQ=255)               cAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCgg                                                                                                                  <  2:582932/140‑1 (MQ=255)                          ttttCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGAc                                                                                                       <  2:602394/140‑1 (MQ=255)                           tttCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACg                                                                                                      >  1:270565/1‑140 (MQ=255)                                tGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGttt                                                                                                 <  2:808780/140‑1 (MQ=255)                                    tAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacc                                                                                             <  2:745207/140‑1 (MQ=255)                                              cGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGa                                                                                   >  1:186985/1‑140 (MQ=255)                                              cGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGa                                                                                   <  2:680103/140‑1 (MQ=255)                                                          tCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                                                       <  2:270565/140‑1 (MQ=255)                                                                   agTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                                              >  1:600946/1‑140 (MQ=255)                                                                         gcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgcc                                                        >  2:923972/1‑140 (MQ=255)                                                                                ctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                                                 <  2:600946/140‑1 (MQ=255)                                                                                                  gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg                               <  2:583898/140‑1 (MQ=255)                                                                                                        tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCg                         <  1:923972/140‑1 (MQ=255)                                                                                                         cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGa                        >  1:550616/1‑140 (MQ=255)                                                                                                         cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGa                        >  1:229184/1‑140 (MQ=255)                                                                                                              gCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc                                  <  2:231497/126‑1 (MQ=255)                                                                                                              gCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc                                  >  1:231497/1‑126 (MQ=255)                                                                                                                  gCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc               >  2:141682/1‑140 (MQ=255)                                                                                                                  gCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc               >  1:558384/1‑140 (MQ=255)                                                                                                                  gCATATTTTGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGc               >  1:257996/1‑140 (MQ=255)                                                                                                                         ttGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCaa        <  2:550616/140‑1 (MQ=255)                                                                                                                         ttGCAGGTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCaa        <  2:558384/140‑1 (MQ=255)                                                                                                                               gTTAACGTTGATCGC‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAActctct  >  1:279818/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                               TCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCT  >  NC_000913/307908‑308172

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