Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,451,787	C→T	A184T (GCC→ACC)	yfcS ←	putative fimbrial chaperone YfcS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,451,787	0	C	T	100.0%	95.9 / NA	27	A184T (GCC→ACC)	yfcS	putative fimbrial chaperone YfcS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/14);  total (13/14)

AGCGGTACCGTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGCGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGTGAGCGTCA  >  NC_000913/2451711‑2451858                                                                             |                                                                        aGCGGTACCGTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTg                                                                     <  2:387982/81‑1 (MQ=255)    ggTACCGTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCt                                                                   >  1:388014/1‑80 (MQ=255)     gTACCGTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGa                                                                 >  2:388019/1‑81 (MQ=255)       aCCGTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGAt                                                                >  1:388011/1‑80 (MQ=255)          gTGGTTTTCGGCTCAATGACCAGAGGCGAGAGCCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATc                                                            <  2:387983/81‑1 (MQ=255)          gTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAAt                                                             >  1:388015/1‑80 (MQ=255)          gTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATc                                                            <  2:387984/81‑1 (MQ=255)             gTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCAc                                                          <  1:387985/80‑1 (MQ=255)                 tCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCAAgtag                                                      >  2:388016/1‑80 (MQ=255)                    gCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACg                                                  <  1:387986/81‑1 (MQ=255)                     cTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACgg                                                 >  2:388017/1‑81 (MQ=255)                           gACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTgg                                           >  1:388028/1‑81 (MQ=255)                               agagGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTa                                       >  2:388018/1‑81 (MQ=255)                                 agGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTAtt                                     <  2:387987/81‑1 (MQ=255)                                      agaACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACgg                                <  1:387988/81‑1 (MQ=255)                                          cccTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTaaa                            >  2:388021/1‑81 (MQ=255)                                           ccTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAAt                           >  2:388023/1‑81 (MQ=255)                                            cTGCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATg                          <  2:387989/81‑1 (MQ=255)                                              gCAGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGcc                        >  1:388013/1‑81 (MQ=255)                                                aGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTg                      <  1:387991/81‑1 (MQ=255)                                                aGCTGGATTACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTg                      <  1:387990/81‑1 (MQ=255)                                                       ttACCGTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTg               <  2:387992/81‑1 (MQ=255)                                                            gTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGt          >  1:388027/1‑81 (MQ=255)                                                            gTTTTTCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGt          <  2:387993/81‑1 (MQ=255)                                                             tttttCTGCGCACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGTg         <  2:387994/81‑1 (MQ=255)                                                                    gcgcACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGTGAGCGTCa  >  2:388020/1‑81 (MQ=255)                                                                     cgcACTGGTGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGTGAGCGTCa  <  2:387995/80‑1 (MQ=255)                                                                             |                                                                        AGCGGTACCGTGGTTTTCGGCTCAATGACCAGAGGCGAGAACCCTGCAGCTGGATTACCGTTTTTCTGCGCACTGGCGTTGCTGATAATCACGTAGTACGGCGTTGGGTTATTGACGGTAAATGCCTGACCGCTGCGGGTGAGCGTCA  >  NC_000913/2451711‑2451858

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: ‑