Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,200,454	(A)8→7	coding (176/3633 nt)	yehI →	DUF4132 domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,200,447	0	A	.	100.0%	149.0 / NA	36	coding (169/3633 nt)	yehI	DUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (20/16);  total (20/16)

AGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCG  >  NC_000913/2200354‑2200566                                                                                              |                                                                                                                        aGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCa                                                                           <  1:99184/139‑1 (MQ=255)             tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                 >  2:1395913/1‑88 (MQ=255)             tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                 <  1:1395913/90‑3 (MQ=255)                ccGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            <  2:1279882/139‑1 (MQ=255)                          cAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAAc                                                  <  1:975595/139‑1 (MQ=255)                                  gTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCgg                                          >  1:1225165/1‑139 (MQ=255)                                        aCTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAaca                                    >  1:1333490/1‑139 (MQ=255)                                           ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTg                                                                                   <  2:885669/89‑1 (MQ=255)                                           ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTg                                                                                   <  2:360355/89‑1 (MQ=255)                                           ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTg                                                                                   >  1:360355/1‑89 (MQ=255)                                           ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTg                                                                                   >  1:885669/1‑89 (MQ=255)                                           ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACaa                                 >  1:1021223/1‑139 (MQ=255)                                               gTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAAccc                             >  2:695181/1‑139 (MQ=255)                                                 ctcATCCAGACATTGGGTGATGAATTGCTTCAAAAAGCTGAGAAG‑TAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             >  1:848635/1‑89 (MQ=255)                                                 ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             <  2:848635/89‑1 (MQ=255)                                                    aTCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        >  1:705337/1‑107 (MQ=255)                                                    aTCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        <  2:705337/107‑1 (MQ=255)                                                     tCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATc                       <  2:1225165/139‑1 (MQ=255)                                                            aTTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggccc                >  1:1373193/1‑139 (MQ=255)                                                              tGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggcccgg              >  2:1208763/1‑139 (MQ=255)                                                              tGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggcccgg              >  1:1114878/1‑139 (MQ=255)                                                                ggTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCATTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGAc            >  2:1009840/1‑139 (MQ=255)                                                                ggTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGAc            >  1:313597/1‑139 (MQ=255)                                                                ggTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGAc            >  2:398868/1‑139 (MQ=255)                                                                     tgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATgg       >  2:1383524/1‑139 (MQ=255)                                                                     tgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATgg       >  2:765206/1‑139 (MQ=255)                                                                     tgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATgg       >  2:411336/1‑139 (MQ=255)                                                                     tgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATgg       >  2:1406975/1‑139 (MQ=255)                                                                     tgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATgg       <  2:1373193/139‑1 (MQ=255)                                                                         ttGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCATTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATggcggc   <  1:1009840/139‑1 (MQ=255)                                                                          tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  1:1208763/139‑1 (MQ=255)                                                                          tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  1:398868/139‑1 (MQ=255)                                                                          tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:1021223/139‑1 (MQ=255)                                                                          tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  1:765206/139‑1 (MQ=255)                                                                            cTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  >  2:1444984/1‑137 (MQ=255)                                                                            cTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  1:1444984/137‑1 (MQ=255)                                                                                       tGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACg                                                              >  2:371532/1‑66 (MQ=255)                                                                                       tGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACg                                                              <  1:371532/66‑1 (MQ=255)                                                                                              |                                                                                                                        AGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCG  >  NC_000913/2200354‑2200566

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