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%	147.5 / NA	46	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 (24/22);  total (24/22)

CGCGCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819356‑3819621                                                                                                                                     |                                                                                                                                      cgcgCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGt                                                                                                                                 >  1:247198/1‑139 (MQ=255)                            tgatgTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcc                                                                                                      >  1:1105367/1‑139 (MQ=255)                            tgatgTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcc                                                                                                      >  2:463817/1‑139 (MQ=255)                                    tAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACt                                                                                              >  2:489353/1‑139 (MQ=255)                                         aTCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         <  1:489353/139‑1 (MQ=255)                                                cccTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttctt                                                                                  >  1:379643/1‑139 (MQ=255)                                                 ccTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTg                                                                                 >  2:770771/1‑139 (MQ=255)                                                    ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                              <  2:677111/139‑1 (MQ=255)                                                     gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                               <  2:312781/89‑1 (MQ=255)                                                     gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                               >  1:312781/1‑89 (MQ=255)                                                          cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                          <  2:798092/89‑1 (MQ=255)                                                          cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                          >  1:798092/1‑89 (MQ=255)                                                             cTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCAttt                                                                     <  1:770771/139‑1 (MQ=255)                                                                ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcc                                                                  >  2:118666/1‑139 (MQ=255)                                                                                      cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           <  1:81470/92‑1 (MQ=255)                                                                                      cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           >  2:81470/1‑92 (MQ=255)                                                                                        ccTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                          <  2:379643/139‑1 (MQ=255)                                                                                            aaGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                          >  1:328379/1‑103 (MQ=255)                                                                                            aaGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                          <  2:328379/103‑1 (MQ=255)                                                                                                   gATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                          >  2:798930/1‑48 (MQ=255)                                                                                                   gATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                          <  1:798930/48‑1 (MQ=255)                                                                                                       aCCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              >  1:24237/1‑104 (MQ=255)                                                                                                       aCCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  2:24237/104‑1 (MQ=255)                                                                                                            tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                          >  1:253678/1‑87 (MQ=255)                                                                                                            tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                          <  2:253678/87‑1 (MQ=255)                                                                                                            tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              >  1:295188/1‑99 (MQ=255)                                                                                                            tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  2:295188/99‑1 (MQ=255)                                                                                                            tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCgg                      >  1:196719/1‑139 (MQ=255)                                                                                                                  actacAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctg                <  1:118666/139‑1 (MQ=255)                                                                                                                   ctacAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctgc               >  2:883741/1‑139 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  1:1036037/89‑1 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              >  1:674822/1‑89 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              >  2:1036037/1‑89 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  2:674822/89‑1 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc            >  2:509087/1‑139 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc            <  2:949473/139‑1 (MQ=255)                                                                                                                        gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                          <  2:860446/107‑1 (MQ=255)                                                                                                                        gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                          >  1:860446/1‑107 (MQ=255)                                                                                                                        gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCa          >  2:753920/1‑139 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          <  2:117171/42‑1 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          >  1:117171/1‑42 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTccc                                    >  2:366849/1‑112 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTccc                                    <  1:366849/112‑1 (MQ=255)                                                                                                                             gAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACt     <  2:452071/139‑1 (MQ=255)                                                                                                                                ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCt  <  2:247198/139‑1 (MQ=255)                                                                                                                                ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCt  <  2:1105367/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      CGCGCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819356‑3819621

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