Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	309,710	C→T	R106H (CGC→CAC)	ecpB ←	putative fimbrial chaperone EcpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,710	0	C	T	100.0%	47.4 / NA	15	R106H (CGC→CAC)	ecpB	putative fimbrial chaperone EcpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/7);  total (8/7)

GCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTC  >  NC_000913/309596‑309837                                                                                                                   |                                                                                                                                gCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCgg                                                                                                               <  1:449694/133‑1 (MQ=255) gCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCgg                                                                                                               <  2:503829/133‑1 (MQ=255)  cgacgaCCAGTATGGCGCCGATGCCGGCGGAATCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCggg                                                                                                              <  2:175991/133‑1 (MQ=255)  cgacgaCCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCggg                                                                                                              <  1:439936/133‑1 (MQ=255)                  gCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                                                              >  1:310352/1‑133 (MQ=255)                      aTGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                                                          >  2:10681/1‑133 (MQ=255)                                                                        cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGc                                        >  2:505008/1‑133 (MQ=255)                                                                          cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCaa                                      >  1:256422/1‑133 (MQ=255)                                                                                     gCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCTGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCt                           <  2:310352/133‑1 (MQ=255)                                                                                     gCCTGATCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGCTAGCAGCATGCTCGCGTCAGTGAGTTGCAACTCATCCgggt                           >  2:241147/1‑131 (MQ=255)                                                                                          aTCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTcc                      >  1:445375/1‑133 (MQ=255)                                                                                           tCAAACCAGACAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCa                     <  2:445375/133‑1 (MQ=255)                                                                                                    aCAATGCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGTGGGAGAGAATAGGACCGCATCCGGCTTGTCGCATGAGAtga            >  2:231475/1‑131 (MQ=255)                                                                                                         gCGGTAGTAGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCGACTGATCAGGCTTGTCCATTGGGATAACCTGc       >  2:225374/1‑133 (MQ=255)                                                                                                            gtagtaGTGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCcg    <  1:505008/133‑1 (MQ=255)                                                                                                              agtagtGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCcgtc  >  2:481529/1‑133 (MQ=255)                                                                                                                   |                                                                                                                                GCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTC  >  NC_000913/309596‑309837

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