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%	157.0 / NA	45	M133R (ATG→AGG)	ydiK	UPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (19/26);  total (19/26)

TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTG  >  NC_000913/1769339‑1769596                                                                                                                                      |                                                                                                                               tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGACTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gggggg                                                                                                                         >  1:851828/1‑139 (MQ=255)      gACGGCAGCGGCCCGCTAGTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gggggctgtct                                                                                                                    >  1:707757/1‑133 (MQ=255)      gACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggggggggacg                                                                                                                    >  1:384043/1‑135 (MQ=255)           cAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGAt                                                                                                                >  2:1036768/1‑139 (MQ=255)                   cGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGa                                                                                                        <  2:667006/139‑1 (MQ=255)                        aTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  1:1055300/139‑1 (MQ=255)                        aTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  2:985021/139‑1 (MQ=255)                          tAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTcc                                                                                                  <  2:991857/139‑1 (MQ=255)                           aaaGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGc                                                                                                <  1:568909/139‑1 (MQ=255)                           aaaGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGc                                                                                                >  2:1011824/1‑139 (MQ=255)                            aaGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGcc                                                                                               <  2:851828/139‑1 (MQ=255)                              gCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCt                                                                                             <  2:384043/139‑1 (MQ=255)                                cATTTCCAGCGGTGACATGACGTTACCCGATCCGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTta                                                                                           >  1:979392/1‑139 (MQ=255)                                  tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   >  2:85639/1‑129 (MQ=255)                                  tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  1:85639/129‑1 (MQ=255)                                         cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCg                                                                                                         >  2:126807/1‑116 (MQ=255)                                         cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCg                                                                                                         <  1:126807/116‑1 (MQ=255)                                         cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcac                                                                                  <  1:901548/139‑1 (MQ=255)                                              aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gggg                                                                                                                           >  2:407792/1‑92 (MQ=255)                                              aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gggg                                                                                                                           <  1:407792/92‑1 (MQ=255)                                                             aTCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTg                                                              >  2:908033/1‑139 (MQ=255)                                                                 ggCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCa                                                          >  1:721669/1‑139 (MQ=255)                                                                  gCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAg                                                         <  2:244916/139‑1 (MQ=255)                                                                   cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        <  1:999445/139‑1 (MQ=255)                                                                   cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        <  1:908033/139‑1 (MQ=255)                                                                        cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  1:515338/91‑1 (MQ=255)                                                                        cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   >  2:515338/1‑91 (MQ=255)                                                                            aTACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCata                                               >  1:443335/1‑139 (MQ=255)                                                                                cATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg                                           <  2:443335/139‑1 (MQ=255)                                                                                       gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCg                                                                                                         >  1:1014848/1‑70 (MQ=255)                                                                                       gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCg                                                                                                         <  2:1014848/70‑1 (MQ=255)                                                                                       gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   >  2:96205/1‑76 (MQ=255)                                                                                       gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  1:96205/76‑1 (MQ=255)                                                                                            tGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAg                                                         <  2:3280/113‑1 (MQ=255)                                                                                            tGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAg                                                         >  1:3280/1‑113 (MQ=255)                                                                                             ggCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   >  2:24250/1‑70 (MQ=255)                                                                                             ggCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  1:24250/70‑1 (MQ=255)                                                                                                 cgAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCAtt                          >  1:691531/1‑139 (MQ=255)                                                                                                        gTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGccc                                                                                              >  1:887933/1‑64 (MQ=255)                                                                                                        gTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGccc                                                                                              <  2:887933/64‑1 (MQ=255)                                                                                                        gTATGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCt                   <  1:1036768/139‑1 (MQ=255)                                                                                                          aTGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 <  2:283820/139‑1 (MQ=255)                                                                                                          aTGCAGGCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 <  1:843944/139‑1 (MQ=255)                                                                                                                gCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGccc                                                                                              <  1:979831/56‑1 (MQ=255)                                                                                                                gCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGccc                                                                                              >  2:979831/1‑56 (MQ=255)                                                                                                                gCTGGCACAACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtc           >  2:820309/1‑139 (MQ=255)                                                                                                                     cacaACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 >  2:244707/1‑128 (MQ=255)                                                                                                                     cacaACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 <  1:244707/128‑1 (MQ=255)                                                                                                                     cacaACTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTc      >  1:54374/1‑139 (MQ=255)                                                                                                                         aCTTGCTGGATATG‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg  <  2:721669/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                               TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTG  >  NC_000913/1769339‑1769596

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