Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,856,453	C→T	*118* (TAG→TAA)	ygbA ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,856,453	0	C	T	100.0%	112.3 / NA	35	*118* (TAG→TAA)	ygbA	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/20);  total (15/20)

CTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856345‑2856583                                                                                                             |                                                                                                                                    cTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTATT‑CGGACGATAtttt                                                                                                      >  2:984605/1‑139 (MQ=255)       cTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGc                                                                                                <  2:812030/139‑1 (MQ=255)       cTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGc                                                                                                <  2:654192/139‑1 (MQ=255)                    cccTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCg                                                                                   >  2:1127298/1‑139 (MQ=255)                      ctgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGTCAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGGACTTCTGGCACCGGa                                                                                 <  1:984605/139‑1 (MQ=255)                             cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGt                                                                          <  2:336774/139‑1 (MQ=255)                              aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCgta                                                                                 >  1:869121/1‑129 (MQ=255)                              aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 <  2:869121/131‑1 (MQ=255)                              aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         <  1:605368/139‑1 (MQ=255)                                   tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    >  1:1461263/1‑139 (MQ=255)                                             gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          <  2:495351/139‑1 (MQ=255)                                             gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          <  2:678478/139‑1 (MQ=255)                                              ggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  2:1461263/139‑1 (MQ=255)                                                   tGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt                                                    <  1:1328890/139‑1 (MQ=255)                                                        aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    <  2:676042/118‑1 (MQ=255)                                                        aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    >  1:676042/1‑118 (MQ=255)                                                          tCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATc                                             >  2:542724/1‑139 (MQ=255)                                                             gCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGa                                          <  1:542724/139‑1 (MQ=255)                                                                  gATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACg                                     >  2:930478/1‑139 (MQ=255)                                                                        aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    <  2:1275935/102‑1 (MQ=255)                                                                        aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    >  1:1275935/1‑102 (MQ=255)                                                                           gCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAg                               <  2:1325720/136‑1 (MQ=255)                                                                           gCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAg                               >  1:1325720/1‑136 (MQ=255)                                                                                         aGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc              <  1:1127298/139‑1 (MQ=255)                                                                                             cTGCCTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACg          <  2:466612/139‑1 (MQ=255)                                                                                                ccTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGc                                                                                                >  2:1288772/1‑50 (MQ=255)                                                                                                ccTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGc                                                                                                <  1:1288772/50‑1 (MQ=255)                                                                                                ccTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCAt       >  2:1408171/1‑139 (MQ=255)                                                                                                 cTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                         <  2:143866/104‑1 (MQ=255)                                                                                                 cTGTCCGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                         >  1:143866/1‑104 (MQ=255)                                                                                                   gTCCGGCATATTATTCATGAGGTTTTTTCCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  1:1123663/87‑1 (MQ=255)                                                                                                   gTCCGGCATATTATTCATGAGGTTTTTTCCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         >  2:1123663/1‑87 (MQ=255)                                                                                                     ccGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         >  1:916254/1‑84 (MQ=255)                                                                                                     ccGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  2:916254/84‑1 (MQ=255)                                                                                                     ccGGCATATTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCt  >  2:1094175/1‑139 (MQ=255)                                                                                                             |                                                                                                                                    CTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTT‑CGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856345‑2856583

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: ‑