Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,190,185	T→C	L102P (CTT→CCT)	yqiH →	putative fimbrial chaperone YqiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,190,185	0	T	C	100.0%	86.3 / NA	26	L102P (CTT→CCT)	yqiH	putative fimbrial chaperone YqiH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/10);  total (16/10)

GTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATG  >  NC_000913/3190109‑3190261                                                                             |                                                                             gTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTcc                                                                           >  2:592625/1‑80 (MQ=255)         aaTGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGa                                                                     >  2:592620/1‑78 (MQ=255)                cGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAAc                                                            <  2:592594/80‑1 (MQ=255)                     aaaTCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCtctt                                                      >  2:592635/1‑81 (MQ=255)                     aaaTCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCtctt                                                      >  2:592615/1‑81 (MQ=255)                      aaTCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCtctt                                                      <  1:592596/80‑1 (MQ=255)                      aaTCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCtctt                                                      <  2:592595/80‑1 (MQ=255)                             gTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTAtttt                                               >  2:592643/1‑80 (MQ=255)                                  cAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTg                                         >  2:592630/1‑81 (MQ=255)                                  cAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTg                                         <  2:592597/81‑1 (MQ=255)                                   aaaTTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGc                                        <  1:592599/81‑1 (MQ=255)                                   aaaTTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGc                                        <  2:592598/81‑1 (MQ=255)                                    aaTTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGc                                        >  1:592652/1‑80 (MQ=255)                                      ttCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCgag                                     <  1:592600/81‑1 (MQ=255)                                        cAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCgaga                                    >  2:592648/1‑80 (MQ=255)                                         aaaTAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGaaa                                  <  1:592601/81‑1 (MQ=255)                                          aaTAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAAt                                 <  1:592602/81‑1 (MQ=255)                                           aTAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGTAACGCTCTTCTATTTTAACGTGCGAGAAATa                                >  2:592640/1‑81 (MQ=255)                                            tAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATAc                               >  2:592634/1‑81 (MQ=255)                                             aaGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATAcc                              >  2:592619/1‑81 (MQ=255)                                                  tGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCa                         <  1:592603/81‑1 (MQ=255)                                                   gATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCa                         >  2:592624/1‑80 (MQ=255)                                                       aCTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaaa                      >  2:592631/1‑79 (MQ=255)                                                             gAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATaa               >  1:592623/1‑80 (MQ=255)                                                                         cAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATg  >  1:592638/1‑81 (MQ=255)                                                                         cAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATg  >  2:592616/1‑81 (MQ=255)                                                                             |                                                                             GTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATG  >  NC_000913/3190109‑3190261

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