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%	136.8 / NA	44	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/20);  total (24/20)

CGCGCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819356‑3819621                                                                                                                                     |                                                                                                                                      cgcgCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTCCTACAGCCACGAAGGACATTAAGGt                                                                                                                                 >  2:1009505/1‑139 (MQ=255) cgcgCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGt                                                                                                                                 >  1:575465/1‑139 (MQ=255) cgcgCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGt                                                                                                                                 >  1:833452/1‑139 (MQ=255)   cgcTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                               <  1:120537/139‑1 (MQ=255)   cgcTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                               >  1:282390/1‑139 (MQ=255)        cAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                          <  1:1009505/139‑1 (MQ=255)                        aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          >  2:94001/1‑139 (MQ=255)                        aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          <  2:760389/139‑1 (MQ=255)                        aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                          >  2:714303/1‑139 (MQ=255)                          cATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTa                                                                                                                                     <  1:858728/110‑1 (MQ=255)                          cATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcag                                                                                                        <  2:658608/139‑1 (MQ=255)                          cATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcag                                                                                                        >  1:658608/1‑139 (MQ=255)                          cATGATGTGGTAATTATCTCCACCCTAGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTa                                                                                                                                     >  2:858728/1‑110 (MQ=255)                            tgatgTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcc                                                                                                      >  1:902986/1‑139 (MQ=255)                             gatgTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcca                                                                                                     >  2:65191/1‑139 (MQ=255)                                    tAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACt                                                                                              >  2:114097/1‑139 (MQ=255)                                       ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                           <  2:282390/139‑1 (MQ=255)                                        tATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGAGCAGCCAGCCAACTGccc                                                                                          >  1:723344/1‑139 (MQ=255)                                         aTCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                         <  2:723344/139‑1 (MQ=255)                                                   tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                               >  2:155288/1‑139 (MQ=255)                                                     gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAttt                                                                             >  2:62071/1‑139 (MQ=255)                                                        ttCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAGACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                          <  1:62071/139‑1 (MQ=255)                                                          cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCa                                                                        >  1:824153/1‑139 (MQ=255)                                                                  ttGCCCTCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcctt                                                                <  2:902986/139‑1 (MQ=255)                                                                    gCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  2:575465/139‑1 (MQ=255)                                                                     cccGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                             >  1:1056143/1‑139 (MQ=255)                                                                           tACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAAt                                                       <  1:94001/139‑1 (MQ=255)                                                                            aCTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              >  2:785078/1‑131 (MQ=255)                                                                            aCTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  1:785078/131‑1 (MQ=255)                                                                             cTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                     <  1:114097/139‑1 (MQ=255)                                                                                    tACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              <  1:134148/123‑1 (MQ=255)                                                                                    tACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                              >  2:134148/1‑123 (MQ=255)                                                                                           tAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGt                                       >  2:60909/1‑139 (MQ=255)                                                                                                            tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCgg                      <  1:65191/139‑1 (MQ=255)                                                                                                             cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                     >  1:960903/1‑139 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc            <  1:714303/139‑1 (MQ=255)                                                                                                                      cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc            >  2:1031675/1‑139 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTTTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAg         <  2:422461/139‑1 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAg         >  1:372810/1‑139 (MQ=255)                                                                                                                         ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGACAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAg         >  1:422461/1‑139 (MQ=255)                                                                                                                            cGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc            >  1:751919/1‑133 (MQ=255)                                                                                                                            cGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc            <  2:751919/133‑1 (MQ=255)                                                                                                                                ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCt  <  1:60909/139‑1 (MQ=255)                                                                                                                                ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCt  <  2:824153/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      CGCGCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819356‑3819621

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