Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,159,005	(T)7→8	coding (2142/2220 nt)	ygiQ ←	radical SAM superfamily protein YgiQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,158,998	1	.	T	100.0%	158.3 / NA	45	A717T (GCG→ACG)	ygiQ	radical SAM superfamily protein YgiQ
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (19/26);  total (19/26)

CTTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGC‑TTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATC  >  NC_000913/3158863‑3159130                                                                                                                                         |                                                                                                                                     ctTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGC‑tttt                                                                                                                                  >  1:692307/1‑140 (MQ=255) ctTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGC‑tttt                                                                                                                                  >  2:619972/1‑140 (MQ=255)  ttGCGTATGTACGAAAAATAAAAAAAAAACCCCCGGTGGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGC‑ttttt                                                                                                                                 <  2:350835/140‑1 (MQ=255)       tATGTACGCAGAATAAAAAAAAAACCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGC‑ttttttttgc                                                                                                                            <  2:549656/140‑4 (MQ=255)           tACGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTa                                                                                                                         >  1:386164/1‑140 (MQ=255)           tACGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTa                                                                                                                         >  2:219824/1‑140 (MQ=255)             cGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGc                                                                                                                       >  2:758964/1‑140 (MQ=255)             cGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGc                                                                                                                       >  1:350835/1‑140 (MQ=255)             cGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGc                                                                                                                       >  1:470578/1‑140 (MQ=255)              gCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCa                                                                                                                      >  2:709577/1‑140 (MQ=255)               cAGAATAAAAAAAAAACCCCGGGTAGCGTGTCAGTTGCCGGGAATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAg                                                                                                                     <  2:658006/140‑1 (MQ=255)                        aaaaaaaTCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGAc                                                                                                            >  1:126567/1‑140 (MQ=255)                         aaaaaaCCCCCGGTAGCGTGTCAGTTGCGGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACg                                                                                                           <  1:589451/140‑1 (MQ=255)                          aaaaaTCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGc                                                                                                          >  2:495094/1‑140 (MQ=255)                            aaaTCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTg                                                                                                        >  2:1090763/1‑140 (MQ=255)                                    ggTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCa                                                                                                <  1:735015/140‑1 (MQ=255)                                    ggTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCa                                                                                                <  1:628239/140‑1 (MQ=255)                                     gTAGCGTGTCGGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCAt                                                                                               <  2:24386/140‑1 (MQ=255)                                     gTAGCGTGTCGGTGGCCGGGGATATTTTTTACCGTCCATCCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCAt                                                                                               <  2:676116/140‑1 (MQ=255)                                     gTAGCGTGTCAGTTGCGGGGGATTTTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCAt                                                                                               <  1:1122441/140‑1 (MQ=255)                                          gtgtCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCg                                                                                          <  1:219824/140‑1 (MQ=255)                                            gtCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTa                                                                                        >  1:676116/1‑140 (MQ=255)                                            gtCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTa                                                                                        >  2:440519/1‑140 (MQ=255)                                                      ggggATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCa                                                                              >  1:1125342/1‑140 (MQ=255)                                                          atatTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGc                                                                          <  2:126567/140‑1 (MQ=255)                                                             ttttttAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCgg                                                                       <  2:470578/140‑1 (MQ=255)                                                               ttttACCGCCAAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGAc                                                                     <  1:228196/140‑1 (MQ=255)                                                                         cAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCg                                                           <  1:758964/140‑1 (MQ=255)                                                                                cgcTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTgg                                                    >  2:78713/1‑140 (MQ=255)                                                                                 gcTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTTCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGc                                                   >  2:58072/1‑140 (MQ=255)                                                                                   tttAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGcga                                                 <  1:619972/140‑1 (MQ=255)                                                                                    ttAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGcgac                                                >  1:1079321/1‑140 (MQ=255)                                                                                     tAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGcgacg                                               <  2:403248/140‑1 (MQ=255)                                                                                       ggctgTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAg                                             >  2:941786/5‑140 (MQ=255)                                                                                         ccGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCt                                           <  1:78713/140‑1 (MQ=255)                                                                                         ccGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCt                                           <  1:440519/140‑1 (MQ=255)                                                                                            tttCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCAGGCGACGAGCTGCa                                        >  2:707744/1‑140 (MQ=255)                                                                                             ttcttcGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTc                                                       <  2:139757/124‑1 (MQ=255)                                                                                             ttcttcGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTc                                                       >  1:139757/1‑124 (MQ=255)                                                                                              tcttcGCACCAGCATTCACCGGACGAGATTTCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACg                                      <  1:58072/140‑1 (MQ=255)                                                                                              tcttcGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACg                                      <  1:559847/140‑1 (MQ=255)                                                                                                 tcGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCAt                                   <  2:386164/140‑1 (MQ=255)                                                                                                     aCCAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATctct                               <  2:692307/140‑1 (MQ=255)                                                                                                      ccAGCATTCACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCtt                              <  1:709577/140‑1 (MQ=255)                                                                                                              cACCGGACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGtt                      <  1:1107684/140‑1 (MQ=255)                                                                                                                   gACGAGATTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGgcgc                 <  1:707744/140‑1 (MQ=255)                                                                                                                      gagaTTGCGTAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAgg              <  2:1125342/140‑1 (MQ=255)                                                                                                                               tAGACGACGCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCa     <  1:941786/140‑1 (MQ=255)                                                                                                                                  acgacgCTTTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATc  >  1:937720/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     CTTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGTAGCGTGTCAGTTGCCGGGGATATTTTTTAACGTCCAACCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTCACCGGACGAGATTGCGTAGACGACGC‑TTTTTTTGCCGTAGCAGGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATC  >  NC_000913/3158863‑3159130

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