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%	263.1 / NA	63	coding (947/978 nt)	ybeQ	Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (29/34);  total (29/34)

GCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGCA  >  NC_000913/675484‑675731                                                                                                                      |                                                                                                                                   gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                              >  1:601630/1‑139 (MQ=255)      tttGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACt                                                                                                         <  1:356076/139‑1 (MQ=255)         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                      <  2:1149451/139‑1 (MQ=255)               tAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTAc                                                                                                >  2:183041/1‑139 (MQ=255)                     tCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAtct                                                                                                                         >  1:979758/1‑108 (MQ=255)                     tCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAtct                                                                                                                         <  2:979758/108‑1 (MQ=255)                             cGCACAGCAATAAAGGGAATATGGGCCGAATGAAATTCCCCGGGCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGGCATCGCTGGTTCCGGATAATTGAACTGTACTATACAGGACgaaa                                                                                  >  2:1260166/1‑137 (MQ=255)                              gCACAGCAATAA‑‑AGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGcc                                                                               >  1:905019/1‑139 (MQ=255)                               cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                <  2:601630/139‑1 (MQ=255)                               cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                >  2:438281/1‑139 (MQ=255)                                   gCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGt                                                                            >  1:738073/1‑139 (MQ=255)                                           gagaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATc                                                                    <  2:265139/139‑1 (MQ=255)                                            agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                      <  2:591471/104‑1 (MQ=255)                                            agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                      >  1:591471/1‑104 (MQ=255)                                             gaATATGCGCCGAATGACATTCCCCGGTCCTTGGTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTGTACGATACAGGTCGAGAGCCTGTCGGAGATCGa                                                                  <  1:1260166/139‑1 (MQ=255)                                              aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                                 >  2:1162122/1‑139 (MQ=255)                                              aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                                 >  1:91179/1‑139 (MQ=255)                                               atatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACa                                                                <  1:438281/139‑1 (MQ=255)                                               atatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACa                                                                <  2:1185331/139‑1 (MQ=255)                                                 atGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAgg                                                              >  2:237648/1‑139 (MQ=255)                                                  tGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGt                                                             <  2:905019/139‑1 (MQ=255)                                                  tGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGt                                                             <  2:950382/139‑1 (MQ=255)                                                   gggCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                      >  2:1008924/1‑97 (MQ=255)                                                   gggCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                      <  1:1008924/97‑1 (MQ=255)                                                   gggCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTc                                                            >  1:952441/1‑139 (MQ=255)                                                   gggCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTc                                                            >  1:1019049/1‑139 (MQ=255)                                                      ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGt                                                                            >  2:4236/1‑120 (MQ=255)                                                      ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGt                                                                            <  1:4236/120‑1 (MQ=255)                                                      ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         <  1:183041/139‑1 (MQ=255)                                                      ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         >  2:875327/1‑139 (MQ=255)                                                        gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGTCAGGTcttcct                                                       >  1:704335/1‑139 (MQ=255)                                                          aTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAtct                                                                                                                         <  1:1011681/71‑1 (MQ=255)                                                          aTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAtct                                                                                                                         >  2:1011681/1‑71 (MQ=255)                                                            gAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTtctc                                                   <  2:952441/139‑1 (MQ=255)                                                              atttCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                 <  2:738073/137‑1 (MQ=255)                                                                ttCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTa                                               <  2:1019049/139‑1 (MQ=255)                                                                        cccTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         <  2:638288/121‑1 (MQ=255)                                                                        cccTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         >  1:638288/1‑121 (MQ=255)                                                                          cTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATa                                     <  1:317254/139‑1 (MQ=255)                                                                             gTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACAGCTCTGGCTCCGGATGACTGAACTTTACGATACAGGGCGAGAGCCTGTCGGAGATCGACAGGGCTTCCTTCTCCGTAGTGATACATAATg                                  >  2:589942/1‑139 (MQ=255)                                                                               ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                                <  1:746073/139‑1 (MQ=255)                                                                               ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                                <  2:426361/139‑1 (MQ=255)                                                                               ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                                <  1:1162122/139‑1 (MQ=255)                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         <  2:162378/105‑1 (MQ=255)                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         >  1:162378/1‑105 (MQ=255)                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTtata                       <  1:875327/139‑1 (MQ=255)                                                                                               aaTCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         >  1:688311/1‑98 (MQ=255)                                                                                               aaTCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                         <  2:688311/98‑1 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                <  1:1100317/69‑1 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                <  1:700753/69‑1 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                >  2:1100317/1‑69 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                <  2:934097/69‑1 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                >  1:934097/1‑69 (MQ=255)                                                                                                     cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                >  2:700753/1‑69 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCgaga                                                                                  <  2:836495/65‑1 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGt                                                             <  2:990643/86‑1 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGt                                                             >  1:990643/1‑86 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                 >  1:188033/1‑98 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                 <  2:188033/98‑1 (MQ=255)                                                                                                       gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTATACGATACAGGTCgaga                                                                                  >  1:836495/1‑65 (MQ=255)                                                                                                        tctcACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCa       <  1:149653/139‑1 (MQ=255)                                                                                                             cGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                                 >  1:842032/1‑76 (MQ=255)                                                                                                             cGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGCa  >  2:123338/1‑139 (MQ=255)                                                                                                             cGAATTTCC‑GACTCACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                                 <  2:842032/76‑1 (MQ=255)                                                                                                                      |                                                                                                                                   GCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGCA  >  NC_000913/675484‑675731

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