Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,194,362	(G)6→7	coding (1499/1662 nt)	yqiK →	PHB family membrane protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,194,356	1	.	G	100.0%	86.2 / NA	23	V498G (GTG→GGG)	yqiK	PHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (11/12);  total (11/12)

GCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCAT  >  NC_000913/3194223‑3194475                                                                                                                                       |                                                                                                                        gCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑ggggg                                                                                                                    >  1:696282/1‑139 (MQ=255)            ggCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCt                                                                                                          >  1:422044/1‑139 (MQ=255)                  gcctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGa                                                                                                    <  2:422044/139‑1 (MQ=255)                   cctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGaa                                                                                                   <  2:696282/139‑1 (MQ=255)                   cctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGaa                                                                                                   <  2:452257/139‑1 (MQ=255)                               gagaAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCa                                                                                       <  1:193142/139‑1 (MQ=255)                                   aaTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCg                                                                                   <  1:556358/139‑1 (MQ=255)                                   aaTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCg                                                                                   <  1:62706/139‑1 (MQ=255)                                    aTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGc                                                                                  >  2:660186/1‑139 (MQ=255)                                    aTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGc                                                                                  >  2:412137/1‑139 (MQ=255)                                        gTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGctctc                                                                              >  2:121380/1‑139 (MQ=255)                                            aGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑ggg                                                                                                                      >  1:437479/1‑94 (MQ=255)                                            aGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑ggg                                                                                                                      <  2:437479/94‑1 (MQ=255)                                                           tCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCAc                                                           >  2:450840/1‑139 (MQ=255)                                                                ggTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTg                                                      >  2:633568/1‑139 (MQ=255)                                                                            attCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGc                                                                                  >  2:748162/1‑99 (MQ=255)                                                                            attCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGc                                                                                  <  1:748162/99‑1 (MQ=255)                                                                            attCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTg                                          <  1:633568/139‑1 (MQ=255)                                                                            attCAGGTCGATGGTCTGAATCATGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTg                                          >  2:434883/1‑139 (MQ=255)                                                                                  gTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAAt                                    >  1:638757/1‑139 (MQ=255)                                                                                         gTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTg                             <  2:638757/139‑1 (MQ=255)                                                                                                 cGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCa                     <  1:681127/139‑1 (MQ=255)                                                                                                   tGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAgg                   <  1:434883/139‑1 (MQ=255)                                                                                                   tGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAgg                   >  2:698374/1‑139 (MQ=255)                                                                                                   tGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAgg                   >  1:103868/1‑139 (MQ=255)                                                                                                                    aTGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCAt  <  2:370081/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                        GCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCAT  >  NC_000913/3194223‑3194475

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