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

GGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTT  >  NC_000913/2856321‑2856586                                                                                                                                     |                                                                                                                                      gggCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCa                                                                                                                                 >  2:1072509/1‑139 (MQ=255) gggCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCa                                                                                                                                 <  2:396418/139‑1 (MQ=255)      tAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAgg                                                                                                                            >  2:658351/1‑139 (MQ=255)                gcgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                                  <  1:172427/139‑1 (MQ=255)                           ggggCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTc                                                                                                       >  2:438635/1‑139 (MQ=255)                             ggCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCg                                                                                                     <  2:1188406/139‑1 (MQ=255)                             ggCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCg                                                                                                     <  1:1551795/139‑1 (MQ=255)                                tGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                                                                  >  1:352205/1‑139 (MQ=255)                                  aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATAtttt                                                                                                         <  2:583627/130‑1 (MQ=255)                                  aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATAtttt                                                                                                         >  1:583627/1‑130 (MQ=255)                                  aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATGCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGt                                                                                                >  2:1241115/1‑139 (MQ=255)                                       cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTgg                                                                                           >  2:1039022/1‑139 (MQ=255)                                            cccTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCg                                                                                      <  1:1072509/139‑1 (MQ=255)                                              ctgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                    <  2:1773022/139‑1 (MQ=255)                                               tgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGtt                                                                                               >  2:136409/1‑127 (MQ=255)                                               tgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGtt                                                                                               <  1:136409/127‑1 (MQ=255)                                               tgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                                                                                   >  1:537872/1‑139 (MQ=255)                                                  tgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCt                                                                                <  1:438635/139‑1 (MQ=255)                                                      aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                            >  2:1977590/1‑139 (MQ=255)                                                            gAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgat                                                                      <  1:976855/139‑1 (MQ=255)                                                               gCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgag                                                                   <  2:537872/139‑1 (MQ=255)                                                                  gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTT‑‑CGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGAc                                                              <  2:692722/139‑1 (MQ=255)                                                                        aaTTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                                  <  1:1339310/83‑1 (MQ=255)                                                                        aaTTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                                  >  2:1339310/1‑83 (MQ=255)                                                                                      cATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                            >  2:539543/1‑139 (MQ=255)                                                                                           aTTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGc                                       <  1:1977590/139‑1 (MQ=255)                                                                                            tttAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGcc                                      <  1:539543/139‑1 (MQ=255)                                                                                             ttAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCca                                     >  1:1483588/1‑139 (MQ=255)                                                                                                     aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGc                                       <  2:1998807/129‑1 (MQ=255)                                                                                                     aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGc                                       >  1:1998807/1‑129 (MQ=255)                                                                                                     aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtc                             >  1:1265168/1‑139 (MQ=255)                                                                                                     aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtc                             >  1:722793/1‑139 (MQ=255)                                                                                                         tGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcgg                         >  2:931547/1‑139 (MQ=255)                                                                                                              gACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                                                                     >  2:1847108/1‑90 (MQ=255)                                                                                                              gACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                                                                     <  1:1847108/90‑1 (MQ=255)                                                                                                                 aGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc                 <  1:1241115/139‑1 (MQ=255)                                                                                                                  gATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGTCCa                <  1:931547/139‑1 (MQ=255)                                                                                                                  gATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCa                <  1:658351/139‑1 (MQ=255)                                                                                                                   aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                           >  2:1682977/1‑95 (MQ=255)                                                                                                                   aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                           <  1:1682977/95‑1 (MQ=255)                                                                                                                          tGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATaa        <  2:352205/139‑1 (MQ=255)                                                                                                                           gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                                    <  1:921280/46‑1 (MQ=255)                                                                                                                           gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                                    >  2:921280/1‑46 (MQ=255)                                                                                                                            tCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATc      <  1:1818539/139‑1 (MQ=255)                                                                                                                             ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                            >  2:1849585/1‑84 (MQ=255)                                                                                                                             ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                            <  1:1849585/84‑1 (MQ=255)                                                                                                                                gCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGc            <  2:1540746/129‑1 (MQ=255)                                                                                                                                gCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGc            >  1:1540746/1‑129 (MQ=255)                                                                                                                                gCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGtt  >  2:1988226/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      GGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTT  >  NC_000913/2856321‑2856586

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