Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	675,601	Δ1 bp	coding (947/978 nt)	ybeQ ←	Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	675,601	0	T	.	100.0%	196.2 / NA	47	coding (947/978 nt)	ybeQ	Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/27);  total (20/27)

GCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCG  >  NC_000913/675484‑675723                                                                                                                      |                                                                                                                           gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                      >  2:123255/1‑139 (MQ=255) gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                      >  2:209933/1‑139 (MQ=255) gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                      >  2:509636/1‑139 (MQ=255)    cGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGa                                                                                                   >  1:70793/1‑139 (MQ=255)         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGTCGAATGAAATTCCCCGGTCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              <  2:495432/139‑1 (MQ=255)         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              <  2:49798/139‑1 (MQ=255)         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              <  2:85018/139‑1 (MQ=255)         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              <  2:243107/139‑1 (MQ=255)               ttaTAGTCTGATTTCTCACAGCAATAAGGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAGATCTTCAGTCTCACGAATTTCC‑GATTTACATCTCTGGTTCCGGACGGCTGACCTTTAc                                                                                        >  2:675714/3‑139 (MQ=255)                 aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGACGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGa                                                                                      <  2:759696/139‑1 (MQ=255)                   aGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATa                                                                                    <  1:209933/139‑1 (MQ=255)                    gTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATAc                                                                                   <  1:639096/139‑1 (MQ=255)                         aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTc                                                                              >  2:598185/1‑139 (MQ=255)                              gCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGagag                                                                         <  2:438054/139‑1 (MQ=255)                               cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                        <  1:654017/139‑1 (MQ=255)                               cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                        <  2:634239/139‑1 (MQ=255)                                            agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCg                                                           <  2:102820/139‑1 (MQ=255)                                            agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCg                                                           <  1:759706/139‑1 (MQ=255)                                              aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                         >  2:13719/1‑139 (MQ=255)                                              aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                         >  2:331756/1‑139 (MQ=255)                                                   gggCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTc                                                    <  1:123255/139‑1 (MQ=255)                                                   gggCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTc                                                    >  2:204094/1‑139 (MQ=255)                                                         aaTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTc                                                    <  1:268068/133‑1 (MQ=255)                                                         aaTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTc                                                    >  2:268068/1‑133 (MQ=255)                                                                   cccGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACa                                                        <  2:670897/119‑1 (MQ=255)                                                                   cccGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACa                                                        >  1:670897/1‑119 (MQ=255)                                                                           ttGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGG‑CGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAAt                           >  1:197278/1‑139 (MQ=255)                                                                               ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        <  1:509636/139‑1 (MQ=255)                                                                               ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        <  1:307700/139‑1 (MQ=255)                                                                               ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        <  1:96997/139‑1 (MQ=255)                                                                                   aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAtct                                                                                                                 >  2:299838/1‑46 (MQ=255)                                                                                   aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAtct                                                                                                                 <  1:299838/46‑1 (MQ=255)                                                                                     ctttACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGAt                                  >  1:583752/4‑123 (MQ=255)                                                                                     ctttACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGAt                                  <  2:583752/120‑1 (MQ=255)                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                        >  2:242633/1‑82 (MQ=255)                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                        <  1:242633/82‑1 (MQ=255)                                                                                               aaTCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAg                                      >  1:56123/1‑109 (MQ=255)                                                                                               aaTCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAg                                      <  2:56123/109‑1 (MQ=255)                                                                                                aTCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCgaga                                                                          <  2:76950/72‑1 (MQ=255)                                                                                                aTCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCgaga                                                                          >  1:76950/1‑72 (MQ=255)                                                                                                    tCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGG‑CGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCg  <  2:197278/139‑1 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTct                                                   >  1:95941/1‑90 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTct                                                   <  2:95941/90‑1 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        <  2:531362/115‑1 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        >  1:531362/1‑115 (MQ=255)                                                                                                                aTTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                                              >  2:599743/1‑84 (MQ=255)                                                                                                                aTTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                                              <  1:599743/84‑1 (MQ=255)                                                                                                                      |                                                                                                                           GCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCG  >  NC_000913/675484‑675723

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