Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,231,366	A→G	S116S (TCA→TCG)	yohK →	PF04712 family membrane protein YohK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,231,366	0	A	G	100.0%	156.4 / NA	44	S116S (TCA→TCG)	yohK	PF04712 family membrane protein YohK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/17);  total (27/17)

GCCTATCCTTTATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATT  >  NC_000913/2231238‑2231500                                                                                                                                 |                                                                                                                                       gCCTATCCTTTATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATcgc                                                                                                                             <  2:402122/140‑1 (MQ=255)           tataTGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATc                                                                                                                   >  1:136883/1‑140 (MQ=255)           tataTGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATc                                                                                                                   >  1:546092/1‑140 (MQ=255)             tatGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCt                                                                                                                 >  2:438651/1‑140 (MQ=255)                gAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGcc                                                                                                              >  1:708916/1‑140 (MQ=255)                   cagcTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGaa                                                                                                           >  1:555078/1‑140 (MQ=255)                      cTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATc                                                                                                        >  1:224414/1‑140 (MQ=255)                      cTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATc                                                                                                        >  1:303402/1‑140 (MQ=255)                       tACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCa                                                                                                       <  2:617979/140‑1 (MQ=255)                            cAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTcac                                                                                                  <  2:669795/140‑1 (MQ=255)                            cAGATCCGCGCGCGCTGGAAATCAATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTcac                                                                                                  <  2:793021/140‑1 (MQ=255)                               aTCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTcaccac                                                                                               <  2:555078/140‑1 (MQ=255)                                tCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACg                                                                                              >  2:459246/1‑140 (MQ=255)                                tCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACg                                                                                              >  2:525629/1‑140 (MQ=255)                                tCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACg                                                                                              >  1:369016/1‑140 (MQ=255)                                tCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACg                                                                                              >  1:62203/1‑140 (MQ=255)                                    cgcgcgcTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTa                                                                                          >  1:517459/1‑140 (MQ=255)                                                    atcatcACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTgg                                                                          <  2:136883/140‑1 (MQ=255)                                                       atcaCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTggcgg                                                                       >  1:542367/1‑140 (MQ=255)                                                       atcaCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTggcgg                                                                       >  1:35940/1‑140 (MQ=255)                                                         caCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCa                                                                     >  2:703013/1‑140 (MQ=255)                                                             aTCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTAt                                                                 <  1:317078/140‑1 (MQ=255)                                                               cTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTggcggc                                                                      >  1:18451/1‑133 (MQ=255)                                                               cTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTggcggc                                                                      <  2:18451/133‑1 (MQ=255)                                                                 gTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGTCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGc                                                             >  1:578353/1‑140 (MQ=255)                                                                 gTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGc                                                             <  2:787288/140‑1 (MQ=255)                                                                 gTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGc                                                             >  1:764678/1‑140 (MQ=255)                                                                   tttATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGGATTGGCgg                                                           >  1:113489/1‑140 (MQ=255)                                                                      aTCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTTGCGGGAt                                                        >  1:386702/1‑140 (MQ=255)                                                                      aTCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATCGGCGGTAt                                                        >  1:722290/1‑140 (MQ=255)                                                                         ggCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCggggtttt                                                     >  2:140771/1‑134 (MQ=255)                                                                          gCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCg                                                    >  1:446022/1‑140 (MQ=255)                                                                            aGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGc                                                  <  2:764678/140‑1 (MQ=255)                                                                              cGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCaa                                                <  2:446022/140‑1 (MQ=255)                                                                              cGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCaa                                                >  1:645814/1‑140 (MQ=255)                                                                               gTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAAt                                               <  2:791477/140‑1 (MQ=255)                                                                                   ttGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTa                                                         >  2:296331/1‑126 (MQ=255)                                                                                   ttGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTa                                                         <  1:296331/126‑1 (MQ=255)                                                                                    tGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGAAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAgcg                                          <  2:303402/140‑1 (MQ=255)                                                                                        aTGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGt                                      <  1:438651/140‑1 (MQ=255)                                                                                                                   ttgATGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGGATTTCCGCCAACAGCGCCGGTTTGGGGGTTTTTGTCGGCGTCCTccg           >  1:63918/1‑138 (MQ=255)                                                                                                                       tGGGCGCTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATCGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGcgcg       <  2:708916/140‑1 (MQ=255)                                                                                                                           cgcTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGGATTGGCGGGATTTCCACAAACAACCCAGGTTTGGTGGTTTTTGTTGGCGTCCTCCGCggggggt   >  2:38003/1‑137 (MQ=255)                                                                                                                            gcTTCGCCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATCGGCGGTATTCCGGCAATCTGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTAtt  <  2:722290/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                       GCCTATCCTTTATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGTTGATGGGCGCTTCACCGGAAATCGCCGCGTCAATCCTGCCGAAATCAGTCACCACGCCTATTGCAATGGCGGTTGGCGGCAGTATTGGCGGTATTCCGGCAATCAGCGCAGTTTGCGTGATTTTCGTCGGCATCCTCGGCGCGGTATT  >  NC_000913/2231238‑2231500

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