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

CCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCG  >  NC_000913/675470‑675723                                                                                                                                    |                                                                                                                           ccAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACa                                                                                                                    <  2:703833/139‑1 (MQ=255)               gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                      >  1:229973/1‑139 (MQ=255)                    tttGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACt                                                                                                 >  2:833112/1‑139 (MQ=255)                             tAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTAc                                                                                        >  1:744293/1‑139 (MQ=255)                               aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                                        <  2:45945/121‑1 (MQ=255)                               aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                                        >  1:364907/1‑121 (MQ=255)                               aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                                        <  2:364907/121‑1 (MQ=255)                               aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATATTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                                        >  1:45945/1‑121 (MQ=255)                                 aGTCTGATTTCGAACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATa                                                                                    <  1:110013/139‑1 (MQ=255)                                     tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTcc‑gactga                                                                                                                      <  1:502557/101‑1 (MQ=255)                                     tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTcc‑gactga                                                                                                                      >  2:502557/1‑101 (MQ=255)                                     tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              <  1:690382/125‑1 (MQ=255)                                     tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              >  2:690382/1‑125 (MQ=255)                                     tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAgg                                                                                >  2:408233/1‑139 (MQ=255)                                           cGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              >  1:479549/1‑119 (MQ=255)                                           cGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTCGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCGGGTTCCGGATGACTGaa                                                                                              <  2:479549/119‑1 (MQ=255)                                                          tgaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                      <  2:545393/95‑1 (MQ=255)                                                          agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                                      >  1:545393/1‑96 (MQ=255)                                                            aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                         >  2:426179/1‑139 (MQ=255)                                                            aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                         <  2:229973/139‑1 (MQ=255)                                                              tatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAg                                                       <  2:744293/139‑1 (MQ=255)                                                                    ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                 >  2:157245/1‑139 (MQ=255)                                                                     cGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcc                                                >  2:455706/1‑139 (MQ=255)                                                                      gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              <  2:610421/92‑1 (MQ=255)                                                                      gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                                              >  1:610421/1‑92 (MQ=255)                                                                              ttCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTa                                       <  1:220254/139‑1 (MQ=255)                                                                                             ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        <  1:408233/139‑1 (MQ=255)                                                                                             ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                        <  1:763510/139‑1 (MQ=255)                                                                                                 aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                    <  1:426179/139‑1 (MQ=255)                                                                                                 aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                    <  1:157245/139‑1 (MQ=255)                                                                                                                cTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATgg                         >  2:481963/1‑119 (MQ=255)                                                                                                                cTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATgg                         <  1:481963/119‑1 (MQ=255)                                                                                                                   cAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCg  >  2:230238/1‑139 (MQ=255)                                                                                                                     gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                    <  2:111240/119‑1 (MQ=255)                                                                                                                     gTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                    >  1:111240/1‑119 (MQ=255)                                                                                                                                    |                                                                                                                           CCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCG  >  NC_000913/675470‑675723

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