Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,769,478	(G)7→8	coding (405/1113 nt)	ydiK →	UPF0118 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,769,471	1	.	G	100.0%	127.3 / NA	36	M133R (ATG→AGG)	ydiK	UPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (17/19);  total (17/19)

CGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTG  >  NC_000913/1769340‑1769601                                                                                                                                     |                                                                                                                                    cgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggggggg                                                                                                                             >  1:1201978/1‑139 (MQ=255)     gACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggggggggacg                                                                                                                         >  1:953649/1‑135 (MQ=255)                       aTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                        <  2:1232040/139‑1 (MQ=255)                       aTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                        <  2:1201978/139‑1 (MQ=255)                        ttAAAACCATTTTCCGCGGTGACATGGCGTTTCCCGATCTTGGGTGGCTTAAAACCATTCCGGTGATTGGCGCGAAGCTTTTTTCCGGCTTGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                        <  1:594360/139‑1 (MQ=255)                         tAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGTCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCg                                                                                                      <  1:769634/139‑1 (MQ=255)                                        cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcac                                                                                       <  2:953649/139‑1 (MQ=255)                                        cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcac                                                                                       <  1:27565/139‑1 (MQ=255)                                                       aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcac                                                                                       >  1:1208914/1‑124 (MQ=255)                                                       aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcac                                                                                       <  2:1208914/124‑1 (MQ=255)                                                       aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGtt                                                                        >  2:383621/1‑139 (MQ=255)                                                       aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggggg                                                                                                                               <  2:449009/83‑1 (MQ=255)                                                       aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggggg                                                                                                                               >  1:449009/1‑83 (MQ=255)                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTGCGTTGGGCAgg                                                             >  1:873054/1‑139 (MQ=255)                                                                       cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCt                                                                                                  >  1:10349/1‑97 (MQ=255)                                                                       cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCt                                                                                                  <  2:10349/97‑1 (MQ=255)                                                                           aTACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCata                                                    >  1:1005632/1‑139 (MQ=255)                                                                            tACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCatat                                                   <  2:1005632/139‑1 (MQ=255)                                                                                      gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGCTATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCttt                                         <  2:873054/139‑1 (MQ=255)                                                                                      gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCttt                                         >  2:265303/1‑139 (MQ=255)                                                                                           tGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt                                    >  2:688004/1‑139 (MQ=255)                                                                                              cgcgAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCa                                 <  1:265303/139‑1 (MQ=255)                                                                                                             aGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg                                                >  2:741446/1‑109 (MQ=255)                                                                                                             aGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg                                                <  1:741446/109‑1 (MQ=255)                                                                                                               gCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca                                                                                  >  1:568278/1‑73 (MQ=255)                                                                                                               gCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca                                                                                  <  2:568278/73‑1 (MQ=255)                                                                                                                  ggCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccac                                                                                 >  1:141430/1‑71 (MQ=255)                                                                                                                  ggCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccac                                                                                 <  2:141430/71‑1 (MQ=255)                                                                                                                  ggCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCtt                 >  1:426628/1‑135 (MQ=255)                                                                                                                  ggCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCtt                 <  2:426628/135‑1 (MQ=255)                                                                                                                  ggCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTctct             >  2:1113918/1‑139 (MQ=255)                                                                                                                  ggCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTctct             >  1:583860/1‑139 (MQ=255)                                                                                                                    cacaACTTGCTGGATATG‑GGGGGGGTCGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACAACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTc           >  1:958454/1‑139 (MQ=255)                                                                                                                    cacaACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGc                                                                                                               <  2:80138/39‑1 (MQ=255)                                                                                                                    cacaACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGc                                                                                                               >  1:80138/1‑39 (MQ=255)                                                                                                                    cacaACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTc           >  2:1163810/1‑139 (MQ=255)                                                                                                                        aCTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg       <  1:1163810/139‑1 (MQ=255)                                                                                                                             cTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTcc                                                                                                       <  2:1145098/38‑1 (MQ=255)                                                                                                                             cTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTcc                                                                                                       >  1:1145098/1‑38 (MQ=255)                                                                                                                             cTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCctg  <  1:688004/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    CGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTG  >  NC_000913/1769340‑1769601

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