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%	137.9 / NA	44	*118* (TAG→TAA)	ygbA	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/18);  total (26/18)

AACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856340‑2856583                                                                                                                       |                                                                                                                                   taCGTCTGCGGCTGAAACGCACGACCCTCCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTTTCCGGCATATTATTCATGAGGTTTTTTCGGACAat                                                                                                           <  1:1221468/138‑1 (MQ=255)               aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCAt                                                                                                                             >  2:443045/1‑107 (MQ=255)               aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCAt                                                                                                                             <  1:443045/107‑1 (MQ=255)                                   aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGc                                                                              >  2:1279894/1‑139 (MQ=255)                                   aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  1:1143866/1‑139 (MQ=255)                                        tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTc                                                                                                    >  1:1329902/1‑107 (MQ=255)                                        tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTc                                                                                                    <  2:1329902/107‑1 (MQ=255)                                        tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         <  2:70134/134‑1 (MQ=255)                                        tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  1:70134/1‑134 (MQ=255)                                        tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    >  1:393710/1‑139 (MQ=255)                                                   ggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  2:1261015/139‑1 (MQ=255)                                                     aaTTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                                                                  <  2:589314/128‑1 (MQ=255)                                                     aaTTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                                                                  >  1:589314/1‑128 (MQ=255)                                                          aTAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                                <  1:1005917/93‑1 (MQ=255)                                                          aTAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                                >  2:1005917/1‑93 (MQ=255)                                                              aTCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCAt                                                                                                                             <  1:498519/60‑1 (MQ=255)                                                              aTCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCAt                                                                                                                             >  2:498519/1‑60 (MQ=255)                                                               tCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                                 >  2:785378/1‑87 (MQ=255)                                                               tCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                                 <  1:785378/87‑1 (MQ=255)                                                                  gCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGa                                          >  1:1478783/1‑139 (MQ=255)                                                                  gCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGa                                          >  2:1151761/1‑139 (MQ=255)                                                                        aTTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                                <  1:1488703/79‑1 (MQ=255)                                                                        aTTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTGTCGGACGATATTTTTCCGGc                                                                                                >  2:1488703/1‑79 (MQ=255)                                                                             aCTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCg                                                                                                  >  1:566007/1‑72 (MQ=255)                                                                             aCTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCg                                                                                                  <  2:566007/72‑1 (MQ=255)                                                                             aCTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCg                                                                     >  2:974295/1‑101 (MQ=255)                                                                             aCTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCg                                                                     <  1:974295/101‑1 (MQ=255)                                                                                 cAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                           >  1:1218016/1‑139 (MQ=255)                                                                                  aaaTTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtc                          >  2:1291666/1‑139 (MQ=255)                                                                                        cc‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          >  1:1057203/1‑101 (MQ=255)                                                                                        cc‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          <  2:1057203/101‑1 (MQ=255)                                                                                                   aGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc              <  1:1279894/139‑1 (MQ=255)                                                                                                       cTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACg          >  1:1537238/1‑139 (MQ=255)                                                                                                        tGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGc         >  2:448982/1‑139 (MQ=255)                                                                                                             gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcccggcc               <  2:992072/128‑1 (MQ=255)                                                                                                             gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcccggcc               >  1:992072/1‑128 (MQ=255)                                                                                                               ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc              <  1:906168/127‑1 (MQ=255)                                                                                                               ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc              >  2:906168/1‑127 (MQ=255)                                                                                                               ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCt  >  2:1098818/1‑139 (MQ=255)                                                                                                               ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCt  >  1:1094271/1‑139 (MQ=255)                                                                                                                cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          >  2:465570/1‑82 (MQ=255)                                                                                                                cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          <  1:465570/82‑1 (MQ=255)                                                                                                                cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc              >  2:44220/1‑126 (MQ=255)                                                                                                                cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc              <  1:44220/126‑1 (MQ=255)                                                                                                                       |                                                                                                                                   AACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCC‑‑‑‑‑GGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856340‑2856583

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