Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	750,061	C→A	M703I (ATG→ATT)	ybgQ ←	putative fimbrial usher protein YbgQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	750,061	0	C	A	100.0%	156.7 / NA	45	M703I (ATG→ATT)	ybgQ	putative fimbrial usher protein YbgQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/23);  total (22/23)

CTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCAGGCG  >  NC_000913/749930‑750192                                                                                                                                    |                                                                                                                                    cTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCAc                                                                                                                             <  1:199401/140‑1 (MQ=255)   ggCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCg                                                                                                                           <  1:328391/140‑1 (MQ=255)      gTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCa                                                                                                                        <  1:604487/140‑1 (MQ=255)         gCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACgg                                                                                                                     <  2:542408/140‑1 (MQ=255)            tATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTgg                                                                                                                  <  1:18953/140‑1 (MQ=255)                   aCGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCa                                                                                                                        >  2:495283/1‑127 (MQ=255)                   aCGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCa                                                                                                                        <  1:495283/127‑1 (MQ=255)                    cGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTgg                                                                                                          >  1:144450/1‑140 (MQ=255)                    cGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTgg                                                                                                          >  2:351678/1‑140 (MQ=255)                     gTTAACGTCTTTCCCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGGCATTAATGTTCACCGcc                                                                                                                         <  2:749924/124‑1 (MQ=255)                     gTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGcc                                                                                                                         >  1:749924/1‑124 (MQ=255)                               ttACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAAt                                                                                               >  1:358131/1‑140 (MQ=255)                                         gggAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGcacc                                                                                     <  2:358131/140‑1 (MQ=255)                                                            aCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAAtggt                                                                  <  1:748689/140‑1 (MQ=255)                                                                 tCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAAtggttggt                                                              <  1:231689/139‑1 (MQ=255)                                                                 tCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAAtggttggt                                                              >  2:231689/1‑139 (MQ=255)                                                                 tCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGgtg                                                             <  2:144450/140‑1 (MQ=255)                                                                 tCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGgtg                                                             >  2:783100/1‑140 (MQ=255)                                                                 tCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGgtg                                                             >  1:29550/1‑140 (MQ=255)                                                                 tCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGgtg                                                             >  1:679575/1‑140 (MQ=255)                                                                    cgcCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCaaa                                                                       <  1:569190/127‑1 (MQ=255)                                                                    cgcCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCaaa                                                                       >  2:569190/1‑127 (MQ=255)                                                                    cgcCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAg                                                          >  2:47503/1‑140 (MQ=255)                                                                    cgcCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATACCaaa                                                                       <  2:475450/127‑1 (MQ=255)                                                                    cgcCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATACCaaa                                                                       >  1:475450/1‑127 (MQ=255)                                                                     gcCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGt                                                         <  2:350644/140‑1 (MQ=255)                                                                           aaTCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGAt                                                   >  2:829375/1‑140 (MQ=255)                                                                                     tCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAAc                                         >  1:726895/1‑140 (MQ=255)                                                                                      ccTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGcc                                               >  1:340965/1‑133 (MQ=255)                                                                                      ccTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGcc                                               <  2:340965/133‑1 (MQ=255)                                                                                           atgaCGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAat                                   <  1:829375/140‑1 (MQ=255)                                                                                           atgaCGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAat                                   >  2:128568/1‑140 (MQ=255)                                                                                           atgaCGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAat                                   >  1:667032/1‑140 (MQ=255)                                                                                               cGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCt                               >  2:740074/1‑140 (MQ=255)                                                                                                 ttttCTGCAACTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGc                             >  2:144773/1‑140 (MQ=255)                                                                                                         aaCTGTTACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCAt                     <  1:128568/140‑1 (MQ=255)                                                                                                              ttACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTg                >  1:838095/1‑140 (MQ=255)                                                                                                               tACGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGc               <  2:29550/140‑1 (MQ=255)                                                                                                                aCGCCGTCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCt              >  2:218898/1‑140 (MQ=255)                                                                                                                      tCGGGTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCAt        <  1:509450/140‑1 (MQ=255)                                                                                                                        gggTAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCa      >  1:432868/1‑140 (MQ=255)                                                                                                                           tAAGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCAGGc   <  2:838095/140‑1 (MQ=255)                                                                                                                            aaGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCAGGCg  <  2:432868/140‑1 (MQ=255)                                                                                                                            aaGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCAGGCg  <  1:47503/140‑1 (MQ=255)                                                                                                                            aaGTCATTAATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCAGGCg  <  1:218898/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    CTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGCCACGGTGGAAGGCTGGTAACTGGAAATCAACGGCACCACCGCAATGCCAAAATGGTTGGTGTAGTCGAGATTGCCCTGAACCGGAATATCTGCCACGCCATCGGTGCTGACCATCAGGCG  >  NC_000913/749930‑750192

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