Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,819,488	C→T	H205H (CAC→CAT)  *561* (TAG→TAA)	yicG → ligB ←	UPF0126 family inner membrane protein DNA ligase, NAD(+)‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,819,488	0	C	T	100.0%	131.5 / NA	41	H205H (CAC→CAT)  *561* (TAG→TAA)	yicG ligB	UPF0126 family inner membrane protein DNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

CTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819360‑3819621                                                                                                                                 |                                                                                                                                      cTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaa                                                                                                                             <  2:867031/139‑1 (MQ=255)          tATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGa                                                                                                                    >  2:1132038/1‑139 (MQ=255)                        tgatgTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcc                                                                                                      >  2:450652/1‑139 (MQ=255)                          atgTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagccag                                                                                                    >  2:517286/1‑139 (MQ=255)                                   ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           <  2:1363868/139‑1 (MQ=255)                                    tATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGccc                                                                                          <  1:905352/139‑1 (MQ=255)                                     aTCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         <  1:450652/139‑1 (MQ=255)                                         ccACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          <  1:425544/118‑1 (MQ=255)                                         ccACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          >  2:425544/1‑118 (MQ=255)                                            cccTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttctt                                                                                  >  1:52884/1‑139 (MQ=255)                                               tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaa                                                                                                                             <  2:1105633/93‑1 (MQ=255)                                               tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaa                                                                                                                             >  1:1105633/1‑93 (MQ=255)                                               tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                               >  1:16331/1‑139 (MQ=255)                                                ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTCCCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCagccagcc                                                                                                  >  2:844918/1‑119 (MQ=255)                                                ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTCCCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCagccagcc                                                                                                  <  1:844918/119‑1 (MQ=255)                                                 gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAttt                                                                             >  2:1481424/1‑139 (MQ=255)                                                      cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           >  2:1347373/1‑120 (MQ=255)                                                      cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           <  1:1347373/120‑1 (MQ=255)                                                      cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCa                                                                        >  1:377503/1‑139 (MQ=255)                                                       ggCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTTAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         >  1:820796/1‑121 (MQ=255)                                                       ggCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTTAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         <  2:820796/121‑1 (MQ=255)                                                              ttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcctt                                                                >  1:371303/1‑139 (MQ=255)                                                                   cGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaaaa                                                                                                                           >  1:1020856/1‑75 (MQ=255)                                                                   cGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaaaa                                                                                                                           <  2:1020856/75‑1 (MQ=255)                                                                              gcTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                                <  1:517286/139‑1 (MQ=255)                                                                                  cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           <  2:111017/92‑1 (MQ=255)                                                                                  cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           >  1:111017/1‑92 (MQ=255)                                                                                       tAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                           <  2:103252/135‑1 (MQ=255)                                                                                       tAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                           >  1:103252/1‑135 (MQ=255)                                                                                         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTcc                                     <  2:16331/139‑1 (MQ=255)                                                                                                        tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTAGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCgg                      >  1:503428/1‑139 (MQ=255)                                                                                                                    gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctg                >  1:5040/1‑133 (MQ=255)                                                                                                                    gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctg                <  2:5040/133‑1 (MQ=255)                                                                                                                     ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         >  2:1136607/1‑59 (MQ=255)                                                                                                                     ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         <  1:1136607/59‑1 (MQ=255)                                                                                                                       aCGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                               >  2:323813/1‑99 (MQ=255)                                                                                                                       aCGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                               <  1:323813/99‑1 (MQ=255)                                                                                                                       aCGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGaa       <  2:371303/139‑1 (MQ=255)                                                                                                                        cGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                          >  2:1036731/1‑103 (MQ=255)                                                                                                                        cGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                          <  1:1036731/103‑1 (MQ=255)                                                                                                                            ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCt  <  1:1481424/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      CTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819360‑3819621

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