Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,362,565	G→A	D152N (GAC→AAC)	yhcA →	putative fimbrial chaperone YhcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,362,565	0	G	A	100.0%	116.1 / NA	33	D152N (GAC→AAC)	yhcA	putative fimbrial chaperone YhcA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/20);  total (13/20)

TACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAA  >  NC_000913/3362488‑3362641                                                                              |                                                                             tACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAAc                                                                            >  2:635248/1‑80 (MQ=255)   ccTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAAcc                                                                           <  1:635204/79‑1 (MQ=255)   ccTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCa                                                                          <  1:635203/80‑1 (MQ=255)    cTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCaaa                                                                        >  1:635226/1‑81 (MQ=255)         cTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTgag                                                                   <  2:635205/81‑1 (MQ=255)           gATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGAt                                                                 >  2:635247/1‑81 (MQ=255)            aTCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCTTGGAATGGCGTAAAAATGGCAACCAAATTGAGAt                                                                 >  2:635243/1‑80 (MQ=255)              ccACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAg                                                              >  2:635260/1‑81 (MQ=255)                aCCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCa                                                            <  1:635206/81‑1 (MQ=255)                   cGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAAcc                                                         <  2:635207/81‑1 (MQ=255)                          ccTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTaa                                                  >  1:635250/1‑81 (MQ=255)                           cTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAAt                                                 <  2:635209/81‑1 (MQ=255)                           cTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAAt                                                 <  1:635208/81‑1 (MQ=255)                              ccGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAAtt                                                <  1:635210/79‑1 (MQ=255)                               cGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATg                                             >  2:635246/1‑81 (MQ=255)                                ggAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGt                                            <  2:635211/81‑1 (MQ=255)                                  aaGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGtcg                                          <  1:635212/81‑1 (MQ=255)                                   agagCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGtcgt                                         <  1:635213/81‑1 (MQ=255)                                        gagagaCCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTa                                    >  1:635252/1‑81 (MQ=255)                                          gagaCCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTAt                                   >  2:635251/1‑80 (MQ=255)                                           agaCCTGGAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTTAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGa                                 >  1:635254/1‑81 (MQ=255)                                                 ggAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTa                            >  2:635241/1‑80 (MQ=255)                                                  gAAATTCCTGGAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTa                            <  2:635214/79‑1 (MQ=255)                                                      ttCCTGGAATGGCGTAGAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAAttt                      <  2:635215/81‑1 (MQ=255)                                                            gAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTaaaa                  <  2:635216/79‑1 (MQ=255)                                                            gAATGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACg                >  2:635262/1‑81 (MQ=255)                                                              aTGGCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCt              <  1:635217/81‑1 (MQ=255)                                                                 gCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTccc           <  1:635219/81‑1 (MQ=255)                                                                 gCGTAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTcc            <  1:635218/80‑1 (MQ=255)                                                                    tAAAAATGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTc        <  2:635220/81‑1 (MQ=255)                                                                        aaTGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGa    <  1:635221/81‑1 (MQ=255)                                                                          tGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGaaa  <  2:635222/81‑1 (MQ=255)                                                                          tGGCAACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGaa   >  1:635258/1‑80 (MQ=255)                                                                              |                                                                             TACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAA  >  NC_000913/3362488‑3362641

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