Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,191,135	G→A	R171C (CGC→TGC)	yehC ←	putative fimbrial chaperone YehC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,191,135	0	G	A	100.0%	103.8 / NA	30	R171C (CGC→TGC)	yehC	putative fimbrial chaperone YehC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/12);  total (18/12)

GACATTAACACTCTGACTTTCTAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATA  >  NC_000913/2191001‑2191271                                                                                                                                       |                                                                                                                                         gACATTAACACTCTGACTTTCTAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTa                                                                                                                                     >  1:281729/1‑140 (MQ=255) gACATTAACACTCTGACTTTCTAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTa                                                                                                                                     >  1:336026/1‑140 (MQ=255)       aaCACTCTGACTTTCTAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCa                                                                                                                               >  2:250855/1‑140 (MQ=255)                  tttCTAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAAt                                                                                                                    <  2:336026/140‑1 (MQ=255)                      tAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCg                                                                                                                >  1:152468/1‑140 (MQ=255)                        aGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCt                                                                                                              >  1:203349/1‑140 (MQ=255)                          ggggCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCttt                                                                                                            >  1:361072/1‑140 (MQ=255)                              cAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCtttattt                                                                                                        >  2:260452/1‑140 (MQ=255)                              cAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCtttattt                                                                                                        >  1:402760/1‑140 (MQ=255)                                          tttCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATa                                                                                            >  1:494134/1‑140 (MQ=255)                                                       tttGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGt                                                                               >  2:172461/1‑140 (MQ=255)                                                           aCATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTaaaa                                                                           <  2:203349/140‑1 (MQ=255)                                                                 ttaGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACtt                                                                     >  1:443208/1‑140 (MQ=255)                                                                  taGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTa                                                                    <  1:172461/140‑1 (MQ=255)                                                                     cTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAAtt                                                                 <  1:64001/140‑1 (MQ=255)                                                                        tGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTg                                                              >  1:5369/1‑140 (MQ=255)                                                                         gACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGt                                                             >  1:92554/1‑140 (MQ=255)                                                                          aCATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGtt                                                            >  2:522533/1‑140 (MQ=255)                                                                               cGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCa                                                       >  2:490297/1‑140 (MQ=255)                                                                                    tCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCa                                                  <  2:494134/140‑1 (MQ=255)                                                                                        cACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCa                                                  <  1:366172/136‑1 (MQ=255)                                                                                        cACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCa                                                  >  2:366172/1‑136 (MQ=255)                                                                                                     gAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGt                                         <  1:54380/132‑1 (MQ=255)                                                                                                     gAGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGt                                         >  2:54380/1‑132 (MQ=255)                                                                                                      aGTCATTTTTTATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTa                                <  2:361072/140‑1 (MQ=255)                                                                                                           ttttttATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACcttc                           <  1:490297/140‑1 (MQ=255)                                                                                                                tATCACCAAACCATTGCCACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGct                      <  1:181450/140‑1 (MQ=255)                                                                                                                                ccACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACC‑TCTTGCTCTGGACTATTTGGTGGa     <  2:81226/140‑1 (MQ=255)                                                                                                                                 cACTGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGa     <  2:5369/140‑1 (MQ=255)                                                                                                                                    tGCAATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAata  >  1:90457/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         GACATTAACACTCTGACTTTCTAAGGGGGCAATCATAATAGTTTCATAATTGACTTTGACATTATTAGCTTTGACATCCGAAATCGTCACCCAATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATA  >  NC_000913/2191001‑2191271

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