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%	186.2 / NA	51	Y657H (TAT→CAT)	ybiW	putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/29);  total (22/29)

GGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCG  >  NC_000913/860512‑860763                                                                                                                               |                                                                                                                              ggTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                   <  2:651845/139‑1 (MQ=255)   ccAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTCTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGAGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCg                                                                                                                 >  1:809435/2‑139 (MQ=255)      gTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAc                                                                                                              <  2:305404/139‑1 (MQ=255)            gAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGAc                                                                                                        >  1:558240/1‑139 (MQ=255)                 cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAg                                                                                                   >  1:1285741/1‑139 (MQ=255)                     cTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAGCGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               <  2:295694/139‑1 (MQ=255)                       tGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGa                                                                                             <  1:298045/139‑1 (MQ=255)                         cgcCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgc                                                                                                  >  2:89901/1‑132 (MQ=255)                         cgcCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgc                                                                                                  <  1:89901/132‑1 (MQ=255)                           ccTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTAt                                                                                         <  2:690225/139‑1 (MQ=255)                              tCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGt                                                                                      >  1:968986/1‑139 (MQ=255)                               cTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTa                                                                                     <  1:1138712/139‑1 (MQ=255)                                        gggTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAg                                                                            <  1:1136919/139‑1 (MQ=255)                                                gCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTgccgcc                                                                                                                      <  2:509401/89‑1 (MQ=255)                                                gCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTgccgcc                                                                                                                      >  1:509401/1‑89 (MQ=255)                                                  tttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                  <  2:1285741/139‑1 (MQ=255)                                                   ttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                                                                 >  2:165176/1‑139 (MQ=255)                                                    tACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATCATGGTACTGTTTCAGTTCGTCGAtata                                                                >  2:628304/1‑139 (MQ=255)                                                     aCGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                               >  2:1299803/1‑139 (MQ=255)                                                               ggTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAg                                                     <  1:492704/139‑1 (MQ=255)                                                                 tgtTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcg                                                   <  1:635746/139‑1 (MQ=255)                                                                 tgtTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcg                                                   <  2:968986/139‑1 (MQ=255)                                                                   ttGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGTTGACGTACCCGCGTAATGGTTGCCGCCAACCgg                                                                                                                >  2:249642/1‑76 (MQ=255)                                                                   ttGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCgg                                                                                                                <  1:249642/76‑1 (MQ=255)                                                                               gcgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                     <  1:1299803/139‑1 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGAtata                                                                <  2:912721/111‑1 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGAtata                                                                >  1:912721/1‑111 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                    >  2:828323/1‑139 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                    >  2:1103000/1‑139 (MQ=255)                                                                                cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                    >  1:673705/1‑139 (MQ=255)                                                                                 gcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                                   <  1:165176/139‑1 (MQ=255)                                                                                     cAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCTTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGAc                               <  2:809435/139‑1 (MQ=255)                                                                                                  aGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                  <  2:1191751/139‑1 (MQ=255)                                                                                                  aGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                  <  2:845194/139‑1 (MQ=255)                                                                                                  aGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                  <  1:1103000/139‑1 (MQ=255)                                                                                                   gCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATcgtcgt                 >  2:167657/1‑139 (MQ=255)                                                                                                    cGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                   <  1:1201029/40‑1 (MQ=39)                                                                                                    cGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                   >  2:1201029/1‑40 (MQ=255)                                                                                                      gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGc              >  2:1274575/1‑139 (MQ=255)                                                                                                            gatgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACt        <  1:167657/139‑1 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 <  2:195339/49‑1 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 >  1:195339/1‑49 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                                                          <  1:935467/56‑1 (MQ=255)                                                                                                             atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                                                          >  2:935467/1‑56 (MQ=255)                                                                                                                  gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATCATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg  <  1:628304/139‑1 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACTGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                          <  1:45842/98‑1 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACTGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                          >  2:45842/1‑98 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 <  2:104022/43‑1 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 >  1:104022/1‑43 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                                   >  1:998199/1‑105 (MQ=255)                                                                                                                   tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                                   <  2:998199/105‑1 (MQ=255)                                                                                                                               |                                                                                                                              GGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCG  >  NC_000913/860512‑860763

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