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

GGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTT  >  NC_000913/2856322‑2856587                                                                                                                                    |                                                                                                                                       ggCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCAt                                                                                                                                 <  1:255308/139‑1 (MQ=255)                      gTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATAtt                                                                                                            >  1:18040/1‑139 (MQ=255)                             gCTGAAACGCACGACCCTGCTGTCACGGATGAAGAGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTTCAGGACAGAGCTGCCTGTGCGGAATATTATTCATGAGGTTTTTTAGGACGATATTTTTCCgg                                                                                                     >  2:761301/1‑139 (MQ=255)                                       aCGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGc                                                                                           >  2:649783/1‑139 (MQ=255)                                        cGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCa                                                                                          <  2:18040/139‑1 (MQ=255)                                           cccTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCg                                                                                       >  1:617139/1‑139 (MQ=255)                                               gctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCg                                                                                       <  2:415461/135‑1 (MQ=255)                                               gctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCg                                                                                       >  1:415461/1‑135 (MQ=255)                                                     aCGGATGAAGCGGCTGGGAATTGATAAATCCGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCTTGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                             >  2:536999/1‑139 (MQ=255)                                                         aTGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCg                                                                         >  2:538998/1‑139 (MQ=255)                                                          tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                        <  2:392158/139‑1 (MQ=255)                                                             aGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGt                                                                                                                            <  2:3180/84‑1 (MQ=255)                                                             aGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGt                                                                                                                            >  1:3180/1‑84 (MQ=255)                                                             aGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                                    >  2:603712/1‑108 (MQ=255)                                                             aGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                                    <  1:603712/108‑1 (MQ=255)                                                                   tGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGAc                                                               <  1:13033/139‑1 (MQ=255)                                                                    gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                              <  1:512387/139‑1 (MQ=255)                                                                    gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                              <  2:151962/139‑1 (MQ=255)                                                                    gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                              <  2:335657/139‑1 (MQ=255)                                                                     ggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                             <  1:316351/139‑1 (MQ=255)                                                                      gAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                            <  2:617139/139‑1 (MQ=255)                                                                      gAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                            <  1:169237/139‑1 (MQ=255)                                                                         ttGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACgg                                                         <  1:761301/139‑1 (MQ=255)                                                                                    gCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGa                                              >  1:202437/1‑139 (MQ=255)                                                                                     cATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCTTGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                             <  1:536999/139‑1 (MQ=255)                                                                                      aTTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                             <  1:574963/138‑1 (MQ=255)                                                                                      aTTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                             >  2:574963/1‑138 (MQ=255)                                                                                               aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAg                                   <  1:538998/139‑1 (MQ=255)                                                                                                   cAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGcc                                       <  1:157293/131‑1 (MQ=255)                                                                                                   cAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGcc                                       >  2:157293/1‑131 (MQ=255)                                                                                                   cAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                               >  1:158289/1‑139 (MQ=255)                                                                                                    aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtc                              >  2:750745/1‑139 (MQ=255)                                                                                                          ccGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcc                        <  1:750745/139‑1 (MQ=255)                                                                                                            ggACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcccg                      <  2:71874/139‑1 (MQ=255)                                                                                                                  aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCaa                >  2:161701/1‑139 (MQ=255)                                                                                                                     tGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACTATATTTTTCTGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                              <  2:475297/90‑1 (MQ=255)                                                                                                                     tGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATAGATGAGATGACg                                                              >  1:475297/1‑90 (MQ=255)                                                                                                                          gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt        >  1:636733/1‑139 (MQ=255)                                                                                                                             cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                            >  2:366424/1‑84 (MQ=255)                                                                                                                             cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                            <  1:366424/84‑1 (MQ=255)                                                                                                                                cATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGttt  <  1:161701/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       GGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTT  >  NC_000913/2856322‑2856587

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