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%	117.9 / NA	33	L102P (CTT→CCT)	yqiH	putative fimbrial chaperone YqiH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/16);  total (17/16)

TGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATG  >  NC_000913/3190108‑3190261                                                                              |                                                                             tGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTcc                                                                           >  2:588677/1‑81 (MQ=255)   ttCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTcc                                                                           <  2:588645/79‑1 (MQ=255)       tccAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGa                                                                     >  1:588682/1‑81 (MQ=255)          aaTGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACag                                                                  <  2:588647/81‑1 (MQ=255)          aaTGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACag                                                                  <  2:588646/81‑1 (MQ=255)                 cGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACg                                                           >  2:588683/1‑81 (MQ=255)                         tCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTa                                                   >  1:588686/1‑81 (MQ=255)                            cAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTAttt                                                <  1:588648/81‑1 (MQ=255)                             aGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTAtttt                                               <  2:588649/81‑1 (MQ=255)                                gAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAAc                                            >  1:588673/1‑81 (MQ=255)                                  gCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGt                                          >  2:588680/1‑81 (MQ=255)                                  gCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGt                                          <  2:588650/81‑1 (MQ=255)                                   cAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTg                                         >  2:588688/1‑81 (MQ=255)                                    aaaTTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGc                                        <  2:588651/81‑1 (MQ=255)                                     aaTTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGc                                        >  2:588697/1‑80 (MQ=255)                                     aaTTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCg                                       <  2:588652/81‑1 (MQ=255)                                        tCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCgag                                     >  1:588681/1‑80 (MQ=255)                                         cAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCgaga                                    >  2:588678/1‑80 (MQ=255)                                           aaTAAGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAAt                                 <  2:588653/81‑1 (MQ=255)                                               aGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATAcc                              <  1:588654/80‑1 (MQ=255)                                               aGAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATAcc                              <  1:588655/80‑1 (MQ=255)                                                gAATGATAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGc                            >  1:588679/1‑81 (MQ=255)                                                     aTAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaa                        <  1:588657/80‑1 (MQ=255)                                                     aTAACTCAGGAAAAAATTGCTCAACCTACTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaa                       <  2:588656/81‑1 (MQ=255)                                                      tAACTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaaa                      >  1:588695/1‑81 (MQ=255)                                                       aaCTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaaa                      >  1:588685/1‑80 (MQ=255)                                                       aaCTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaaa                      >  1:588658/1‑80 (MQ=255)                                                       aaCTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaaa                      >  1:588687/1‑80 (MQ=255)                                                         cTCAGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCaaaa                      >  1:588690/1‑78 (MQ=255)                                                            aGGAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATa                <  1:588659/81‑1 (MQ=255)                                                              gAAAAAATTGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATaaa              <  2:588660/81‑1 (MQ=255)                                                                     ttGCTCAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAAtg       <  2:588661/81‑1 (MQ=255)                                                                          cAACCTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATg  >  1:588691/1‑81 (MQ=255)                                                                              |                                                                             TGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCAAAACAGATAAAAAAAATGTGATG  >  NC_000913/3190108‑3190261

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