Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	860,638	A→G	Y657H (TAT→CAT)	ybiW ←	putative pyruvate formate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	860,638	0	A	G	100.0%	188.7 / NA	52	Y657H (TAT→CAT)	ybiW	putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/27);  total (25/27)

GGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCA  >  NC_000913/860512‑860765                                                                                                                               |                                                                                                                                ggTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTCACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                     <  2:1311965/139‑1 (MQ=255)    cAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCgg                                                                                                                  >  2:756762/1‑139 (MQ=255)      gTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAc                                                                                                                <  1:514497/139‑1 (MQ=255)      gTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAc                                                                                                                <  2:788085/139‑1 (MQ=255)        aCCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCa                                                                                                              <  2:1522212/139‑1 (MQ=255)           ggAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGa                                                                                                           <  1:740756/139‑1 (MQ=255)                 cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAg                                                                                                     >  2:190579/1‑139 (MQ=255)                 cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAg                                                                                                     >  2:1547566/1‑139 (MQ=255)                  gggCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgc                                                                                                    <  2:1560495/139‑1 (MQ=255)                  gggCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgc                                                                                                    <  2:1244922/139‑1 (MQ=255)                               cTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTa                                                                                       <  1:1547566/139‑1 (MQ=255)                                        gggTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAg                                                                              >  1:95628/1‑139 (MQ=255)                                          gtgtGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACTCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGtt                                                                            >  1:1156740/1‑139 (MQ=255)                                          gtgtGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGtt                                                                            >  2:1598402/1‑139 (MQ=255)                                                  tttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    <  2:305175/139‑1 (MQ=255)                                                  tttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    <  1:845657/139‑1 (MQ=255)                                                   ttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                                                                   >  2:760005/1‑139 (MQ=255)                                                    tACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGAtata                                                                  <  1:760005/139‑1 (MQ=255)                                                            tCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGt                                                                                                                               >  1:1153102/1‑70 (MQ=255)                                                                          aGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                            >  2:354821/1‑139 (MQ=255)                                                                           gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                     <  2:185317/65‑1 (MQ=255)                                                                           gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                     >  1:185317/1‑65 (MQ=255)                                                                           gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCa                                                                                                              <  1:1325459/72‑1 (MQ=255)                                                                           gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCa                                                                                                              >  2:1325459/1‑72 (MQ=255)                                                                              tGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    <  2:670621/111‑1 (MQ=255)                                                                              tGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    >  1:670621/1‑111 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    <  2:892956/109‑1 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    >  1:892956/1‑109 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                      >  1:1150231/1‑139 (MQ=255)                                                                                     cAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGa                                                                                               <  1:153500/77‑1 (MQ=255)                                                                                     cAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGa                                                                                               >  2:153500/1‑77 (MQ=255)                                                                                      aaaCGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                                <  2:1156740/139‑1 (MQ=255)                                                                                                 tAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                   <  2:1381265/61‑1 (MQ=255)                                                                                                 tAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                   >  1:1381265/1‑61 (MQ=255)                                                                                                 tAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtc                     >  1:1376059/1‑139 (MQ=255)                                                                                                      gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGc                >  2:1609714/1‑139 (MQ=255)                                                                                                         aTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                    >  2:845861/1‑132 (MQ=255)                                                                                                         aTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                    <  1:845861/132‑1 (MQ=255)                                                                                                         aTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGt             >  2:454280/1‑139 (MQ=255)                                                                                                         aTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGt             >  2:52374/1‑139 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                                 >  2:781085/1‑51 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                                 <  1:781085/51‑1 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACtt         <  1:756762/139‑1 (MQ=255)                                                                                                                 cGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgc     <  2:1150231/139‑1 (MQ=255)                                                                                                                  gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg    <  1:1598402/139‑1 (MQ=255)                                                                                                                  gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg    <  1:190579/139‑1 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                   >  1:1152161/1‑43 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                   <  2:1152161/43‑1 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    <  1:950879/74‑1 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                    >  2:950879/1‑74 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcgc   >  1:49963/1‑139 (MQ=255)                                                                                                                    aCCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCa  <  2:1376059/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                GGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCA  >  NC_000913/860512‑860765

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